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Day MONTH Year 
(use item 2 if no earlier priority) 
Measured from the earliest priority date in item 3, this PCT/USA National Phase Application Request is being filed 
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(a) [XX ] 20 months from above item 3 date (b) [ ] 30 months from above item 3 date, 
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Title of Invention METHOD OF SIMULATING INDOOR BEHAVIOR OF PESTICIDAL COMPOUND 



6^ Inventor(s) MATOBA. Yoshihide et al 

Applicant herewith submits the following under 35 U.S.C. 371 to effect filing: 

7^: [ X ] Please immediately start national examination procedures (35 U.S.C. 371(f)). 



8r [ ] A copy of the International Application as filed (35 U.S.C. 371(c)(2)) is transmitted herewith (file if in 
■■^ -English but, if in foreign language, file only if not transmitted to PTO by the International Bureau) including: 

a. [ ] Request; 
& [ ] Abstract; 

c. pgs. Spec, and Claims; 

d. sheet(s) Drawing which are [ ] informal [ ] formal of size [ ] A4 [ ] 13" [ ] 14" 

9. [ X ] A copy of the International Application has been transmitted by the International Bureau. 

10. A translation of the International Application into English (35 U.S.C. 371(c)(2)) 

a. [ XX] is transmitted herewith including: (1) [ ] Request; (2) [XX ] Abstract; 

(3) 265 pgs. Spec, and Claims; 

(4) 21 sheet(s) Drawing which are: 

[ ] informal [XX ] formal of size [XX ] A4 [ ] 11" 

b. [ ] is not required, as the application was filed in English. 

c. [ ] is not herewith, but will be filed when required by the forthcoming PTO Missing Requirements Notice 

per Rule 494(c) if box 4(a) is X'd or Rule 495(c) if box 4(b) is X'd. 

d. [ ] Translation verification attached (not required now). 

1 1 ■ I X ] PLEASE AMEND the specification before its first line by inserting as a separate paragraph: 

a. IX] -This application is the national phase of international application PCT / JP97 /03924 

filed October 29. 1997 which designated the U.S.- 

b. [ ] -This application also claims the benefit of U.S. Provisional Appliction No 

60/ , filed .-- 
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Re: USA National Filing of PCT/ JP97 / 03924 p age 2 of 3 

12. [ ] Amendments to the claims of the International Application under PCT Article 19 (35 U.S.C. 371(c)(3)), 

i e., before 18th month from first priority date above in item 3, are transmitted herewith (file if in 
English but, if in foreign language, file only if not transmitted by the International Bureau) 
including: 

13. [X] PCT Article 1 9 claim amendments (if any) have been transmitted by the International Bureau. 

14. [ ] Translation of the amendments to the claims under PCT Article 19 (35 U.S.C. 371(c)(3)), i.e., of claim 

amendments made before 18th month, is attached (required by 20th month from the date in item 3 
if box 4(a) above is X'd. or 30th month if box 4(b) is X'd. or else amendments will be considered 
cancelled) . 

15. A declaration of the inventor (35 U.S.C. 371(c)(4)) 

a. [XX ] is submitted herewith [ ] Original [XX] Facsimile/Copy 

b. [ ] is not herewith, but will be filed when required by the forthcoming PTO Missing Requirements Notice 

per Rule 494(c) if box 4(a) is X'd or Rule 495(c) if box 4(b) is X'd. 

16. An International Search Report (ISR): 

a. Was prepared by [ ] European Patent Office [XX ] Japanese Patent Office [ ] Other 

b. [ X ] has been transmitted by the International Bureau to PTO. 

c. [XX ] copy herewith ( 2 pg(s).) [ ] plus Annex of family members ( pg(s).). 

1U International Preliminary Examination Report (IPER): 

S=^ a. [ X ] has been transmitted (if this letter is filed after 28 months from date in item 3) in English by the 
a International Bureau with Annexes (if any) in original language. 

|U b. [ ] copy herewith in English 

hj c.1 [ ] IPER Annex(es) in original language ("Annexes" are amendments made to claims/spec/drawings 
during Examination) including attached amended: 

|2 c.2 [ ] Specification/claim pages # [ ] Drawing Sheets # 

-•== 5 c.3 [ ] Which resulted in cancellation of pages # claims # 

Dwg Sheets # 

Q d. [ ] Translation of Annex(es) to IPER ( required by 30th month due date, or else annexed 
HJ amendments will be considered cancelled) . 

18. ; Information Disclosure Statement including: 

? 1J a. [ ] Attached Form PTO-1449 listing documents 

b. [ ] Attached copies of documents listed on Form PTO-1449 
yj c. [ X ] A concise explanation of relevance of ISR references is given in the ISR. 

19. [XX] Assignment document and Cover Sheet for recording are attached. Please mail the recorded 

assignment document back to the person whose signature, name and address appear at the end of this 
letter. 

20. [ ] Copy of Power to IA agent. 

21. [ ] Drawings: sheet(s) per set: [ ] 1 set informal; [ ] Formal of size [ ] A4 [ ] 11" 

22. [ ] (No.) Verified Statement(s) establishing "small entity" status under Rules 9 & 27 

23. Priority is hereby claimed under 35 U.S.C. 119/365 based on the priority claim and the certified copy, 
both filed in the International Application during the international stage based on the filing 

in (country) JAPAN of: 

Application No. Filing Date Application No. Filing Date 

(1) 294203/1996 Nov. 6. 1996 (4) 299073/1996 Nov. 11. 1996 

(2) 295366/1996 Nov. 7. 1996 (5) 

(3) 299059/1996 Nov. 11. 1996 (6)_ 



[ X ] See Form PCT/IB/304 sent to US/DO with copy of priority documents. If copy has not been received, 
please proceed promptly to obtain same from the IB . 
b. [ ] Copy of Form PCT/IB/304 attached. 

24. Attached: VERIFICATION 

25. Preliminary Amendment: ATTACHED 
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Re: USA National Filing of PCT/ JP97 / 03924 Page 3 of 3 
25.5 Per item 17.c3, cancel original pages # claims # , Drawing Sheets # 

26. Calculation of the U.S. National Fee (35 U.S.C. 371 (c)(1)) and other fees is as follows : 

based on amended claim(s) per above item(s) [ ] 12, [ ] 14, [ ] 17, [ XX] 25 [ ] 25.5 (hilite) 

Large/Small Entity fe g»fc 
TOTAL EFFECTIVE CLAIMS 36 -20 = * 16 x $ 22/$ll = $ 352.00 pawm 

INDEPENDENT CLAIMS 2 - 3= * 0 x $ 82/$41 = $ 

*lf answer <0, enter "0" 

If any proper (ignore improper) MULTIPLE DEPENDENT CLAIM is present, add $270/$135 + mm 

BASIC NATIONAL FEE (37 CFR 1 .492(a)(1)-(4)): >B ASIC FEE REQUIRED, NOW ^ * 
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A. If country code letters in item 1 are not "US", "BR". "BB". "TT", "MX*. IL or "NZ' I 

See item 16 re: I 

1. Search Report was not prepared by EPO or JPO — - add $1070/$535 + <96Q9sd 

2. Search Report was prepared by EPO or JPO add $930/ $465 + 930.00 <9rom> 

SKIP B, C, D AND E UNLESS country code letters in item 1 are "US". "BR". "BB". "TT", "MX". "IL" or "NZ 

-> [ ] B. If neither international search fee nor international 

preliminary examination fee was paid to USPTO , add $1070/$535 + mm 

( 30 

( only) -> [ ] C. If international search fee was paid to USPTO 

( one) but not international preliminary examination fee, add $790/ $395 + 0ssm 

< of) 

(these)-> [ ] D. If international preliminary examination fee was paid to 

( 4) USPTO - add $720/$360 + earn 

(boxes) 

-> [ ] E. If international pre liminaiy examination fee was paid 

to USPTO and Rules 492(a)(4) and 496(b) satisfied. add $98/ $49 + ( 962«3) 

27 SUBTOTAL = $ 1,282.00 

28. If Assignment box 19 above is X'd, add Assignment Recording fee of $40.00 + 40.00 (58» 

29. Attached is a check to cover the TOTAL FEES $ 1,322.00 

CHARGE STATEMENT: The Commissioner is hereby authorized to charge any fee specifically authorized hereafter, or any missing or insufficient fee(s) 
filed, or asserted to be filed, or which should have been filed herewith or concerning any paper filed hereafter, and which may be required under Rules 
16-18 and 492 (missing or insufficient fee only) now or hereafter relative to this application and the resulting Official document under Rule 20, or credit 
any overpayment, to our Account/Order Nos. shown in the heading hereof for which purpose a duplicate copy of this sheet is attached. 
This CHARGE STATEMENT does not authorize charge of the issue fee until/unless an issue fee transmittal form 
is filed. 

Pillsbury Madison & Sutro LLP 
Intellectual Property Group 

1100 New York Avenue 

Ninth Floor, East Tower By Atty: Kendrew H. Colton Reg. No. 30,368 

Washington, D.C. 20005-3918 

Tel: (202) 861-3000 A Qf r I A^Q ft J^/ Fax: (202)822-0944 

Attv/Sec: KHC:mhn Sio: t&~X < C^/JjU Te , . r2rm afi , o finfi 



Atty/Sec: KHC:mhn sig: f ^ > ^S^*; Tel.: (202) 861- 3606 

NOTE: File in duplicate with 2 postcarcf/eceipts (PAT-103) & attachments. 
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In re Patent Application of : MATOBA, Yoshihide et al 
Atty.DKT. : M.# 255089 
Filed : Herewith 

Title: METHOD OF SIMULATING INDOOR BEHAVIOR OF PESTICIDAL 
COMPOUND 

*************** 

Date: July 6, 1998 

PRELIMINARY AMENDMENT 

Honorable Commissioner 
Patent and Trademark Office 
Washington, D.C. 20231 

Dear Sir: 

Prior to rendering an initial Official Action in this 
application, please amend it as follows: 

In the Specification 

Page 2 07 

The examples of values in Table 35 assume those of a 
six-mat room (9.72 [m 2 ]) of a typical apartment in Japan. 

The cross-sectional area (AJ of the space portion of 
the carpet in Table 36 is an area in which the space portion 
is in contact with the indoor air (see Figs. 13A, 13B, and 
13C) . Further, the examples of values in Table 36 are based 
on the values shown in Figs. 13A, 13B, and 13C. Here, one 
space portion is set to a size of 0.1 [mm] x 0.1 [mm] x 3 
[mm], and one square of carpet fiber is set to a size of 
2.5 [mm] x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 
[pieces] of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as: 



Ac = (0.1 [ram] x o.l [mm] * 16 [pieces] )/l [cm] 

= 1.6 x 10 -3 [m 2 /m 2 floor] 
A cS = (0.1 [mm] x 3 [mm] * 4 [faces] x 16 [pieces] )/l 

[cm] 

= 0.192 [mVm 2 floor] 
As = 1 - K 

= 1 - 1.6 x 10" 3 [mVm 2 floor] 
A a5 = A s + A^ 

= 1.19 [mVm 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2 ) . The 
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DESCRIPTION 

Method of Simulating Indoor Behavior of Pesticidal 
Compound 
Technical Field 

The present invention relates to a method of 
simulating, when a chemical containing a pesticidal 
compound is used indoors, an indoor behavior of the 
pesticidal compound including an estimation method of 
estimating the indoor behavior of the compound and a safety- 
evaluating method of evaluating its safety in human bodies 
by using the estimation method; and, in particular, to a 
method of simulating an indoor behavior of a pesticidal 
compound in the case where a chemical containing the 
pesticidal compound is residually sprayed, sprayed in an 
indoor space, heated to vaporize, or sprayed over the whole 
floor area. 
Ba ck g rou nd Art 

Conventionally known is a f ugaciousness (hereinafter 
referred to as Fugacity) model for simulating a behavior 
of a chemical material in global environment. The 
above-mentioned fugacity model utilizes Fugacity whose unit 
is an external force by which the chemical material escapes 
from one medium to another medium, i.e., pressure. For 
example, when the chemical material concentration differs 
between two media A and B, concentrations in the respective 
media are expressed by: 
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N A /V A = f A Z A 
N B /V B = f B Z B 

wherein N is chemical mass, V is volume of medium, f is 
Fugacity, and Z is Fugacity capacity of medium. 

Here, while the mass N changes over time according to 
transference and degradation of the chemical material 
between the media A and B; assuming that volume V and Fugacity 
capacity Z are constant, the above-mentioned two 
expressions are represented as: 

(df A /dt)V A Z A = dN A /dt 

= - ( Degradation ) A ± ( Trans ference) M 

(df B /dt)V B Z B = dN B /dt 

= - ( Degradation ) B ± (Transf erence) AB 

When Degradation and Transference in these two 
differential equations are given, unknown parameters f A and 
f B can be determined by calculation. When these parameters 
are respectively multiplied by Fugacity capacity Z A and Z B/ 
the chemical material concentrations in the respective 
media in a specific period of time can be simulated. 

As an apparatus for simulating a behavior of a chemical 
material, Japanese Patent Application Laid-Open No. 64- 
88811 discloses a configuration of closed-space simulator 
which can perform simulation in response to any capacity 
of closed space without actually constructing a closed space 
when evaluating temperature change of a specific gas 
component such as carbonic acid gas. 



2 



SCC 97-06 



As a configuration for evaluating influence of a 
harmful material on human bodies, Japanese Patent 
Application Laid-Open No. 3-89146 discloses a configuration 
of percutaneous absorption evaluation apparatus employing 
a vertical type diffusion cell which is a system closer to 
a clinical state than is a parallel type cell, thereby being 
capable of simultaneously measuring, in real time, the 
process of a chemical being emitted from its base on the 
skin by optoacoustic measurement and the process of the 
chemical infiltrating through the skin by absorptiometry. 

Also, Japanese Patent Application Laid-Open No. 7- 
218496 discloses a configuration of system which uses the 
fact that a dissolution parameter inherently existing in 
a chemical material and a logarithmic value of median lethal 
dose of the chemical material with respect to a mammal are 
in a specific correlation, estimating acute toxicity of the 
chemical material with respect to the mammal. 

Further, simulations of indoor behavior of a 
pesticidal compound when insecticides are sprayed in an 
indoor space, electrically heated to vaporize in a room, 
and sprayed over the whole floor surface are respectively 
disclosed in Y . Matoba et al . , "A SIMULATION OF INSECTICIDES 
IN INDOOR AEROSOL SPACE SPRAYING," Chemosphere r Vol.26, 
No. 6, pp. 1167-1186, 1993; Y. Matoba et al., "INDOOR 
SIMULATION OF INSECTICIDES SUPPLIED WITH AN ELECTRIC 
VAPORIZER BY THE FUGACITY MODEL , " Chemosphere , Vol.28, No . 4 , 
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pp. 767-786, 1994; and Y. Matoba et al., "INDOOR SIMULATION 
OF INSECTICIDES IN BROADCAST SPRAYING," Chemosphere f 
Vol.30, No. 2, pp. 345-365, 1995. 
Disclosure of Invention 

The above-mentioned simulation models, however, do 
not mention how to solve differential equations, and minute 
time units set when solving the differential equations are 
assumed to be constant. Theoretically, the smaller is the 
minute time unit, the longer becomes the calculation time; 
whereas the solution would not converge when the minute time 
unit is large. Accordingly, in the case where a 
differential equation containing a parameter which changes 
over time is to be solved, when the minute time unit is set 
to a constant value so that the solution does not diverge, 
there is a problem that the processing speed of a computer 
must be enhanced. 

Also, the above-mentioned simulation models fail to 
mention any security with respect to human bodies. 

In order to solve the above-mentioned conventional 
problems, it is an object of the present invention to provide 
a method of simulating an indoor behavior of a pesticidal 
compound, which can process simultaneous differential 
equations accurately in a short time by automatically 
setting a minute time unit. 

In order to achieve the above-mentioned object, the 
method of simulating an indoor behavior of a pesticidal 



SCC 97-06 



compound in accordance with the present invention comprises 
a step of dividing an indoor environment into predetermined 
media (constituents) and forming a differential equation 
concerning a fugacity of the compound in each of the media; 
5 a step of determining the fugacity of the compound in each 

of the media from the differential equation; a step of 
determining the indoor behavior of the compound from the 
fugacity of the compound in each of the media; and a step 
of changing, in response to a fluctuation in mass balance 

10 of the compound indoors , a minute time unit used when solving 

the differential equation. 

As the indoor environment is divided into 
predetermined media, and exchanges of the chemical compound 
between the media and the like are taken into account, 

15 simulation results close to the actual behavior of the 

compound can be obtained, while the minute time unit can 
be set automatically in response to fluctuation in mass 
balance when solving simultaneous differential equations 
including a parameter which changes over time . Accordingly , 

2 0 when a computer processes the above-mentioned differential 

equation, accurate solutions can be obtained in a short time. 

Preferably, the method of simulating an indoor 
behavior of a pesticidal compound in accordance with the 
present invention further comprises a step of evaluating 

25 safety of the compound with respect to a human body according 

to the indoor behavior of the compound. 
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As a consequence of this configuration, the safety of 
the pesticidal compound with respect to the human body can 
be evaluated accurately in a short time. Accordingly, when 
formulating a chemical such as insecticide including the 
above-mentioned compound, simulation can be easily repeated 
while changing conditions, thereby making it easier to 
formulate a chemical having a high safety conforming to the 
aimed object. 

Further, the method of simulating an indoor behavior 
of a pesticidal compound in accordance with the present 
invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is residually sprayed; whereas the above- 
mentioned media are a spraying site, suspended particles 
which are divided into at least one kind according to size, 
indoor air, a floor, a wall, and a ceiling. 

Preferably, in this case, the differential equation 
at the spraying site is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound at the spraying site, temporal change in volume 
of the spraying site, amount of attachment of the suspended 
particles to the spraying site, amount of transference of 
the compound between the spraying site and another medium, 
and change in amount of degradation of the compound at the 
spraying site; the differential equation in the suspended 
particles is a differential equation stating a relationship 
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among temporal change of fugacity of the compound in the 
suspended particles, temporal change in volume of the 
suspended particles, amount of transference of the compound 
between the suspended particles and another medium, and 
change in amount of degradation of the compound in the 
suspended particles; the differential equation in the 
indoor air is a differential equation stating a relationship 
among temporal change of fugacity of the compound in the 
indoor air, amount of discharge of the compound outdoors, 
amount of transference of the compound between the indoor 
air and another medium, and change in amount of degradation 
of the compound in the indoor air; the differential equation 
at the floor is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound at the floor, temporal change in volume of the floor, 
amount of attachment of the suspended particles to the floor, 
amount of transference of the compound between the floor 
and another medium, and change in amount of degradation of 
the compound at the floor; the differential equation at the 
wall is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the wall, 
temporal change in volume of the wall, amount of attachment 
of the suspended particles to the wall, amount of 
transference of the compound between the wall and another 
medium, and change in amount of degradation of the compound 
at the wall; and the differential equation at the ceiling 
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is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the ceiling, 
temporal change in volume of the ceiling, amount of 
attachment of the suspended particles to the ceiling, amount 
of transference of the compound between the ceiling and 
another medium, and change in amount of degradation of the 
compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is spatially sprayed; whereas the above- 
mentioned media are suspended particles which are divided 
into at least one kind according to size, indoor air, a floor, 
a wall, and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the suspended particles, temporal change 
in volume of the suspended particles, amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 
in the suspended particles; the differential equation in 
the indoor air is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
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compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 
amount of degradation of the compound in the indoor air; 
the differential equation at the floor is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
particles to the floor, amount of transference of the 
compound between the floor and another medium, and change 
in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound at the wall, temporal change in volume of 
the wall, amount of attachment of the suspended particles 
to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 
equation at the ceiling is a differential equation stating 
a relationship among temporal change of fugacity of the 
compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 
of degradation of the compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
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invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is heated to vaporize; whereas the above- 
mentioned media are condensed particles which are divided 
into at least one kind according to generation and extinction, 
high-concentration air, medium-concentration air, low- 
concentration air, a floor, a wall, and a ceiling which is 
divided into at least one kind according to compound 
concentration . 

Preferably, in this case, the differential equation 
in the condensed particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the condensed particles, temporal change 
in volume of the condensed particles, amount of transference 
of the compound between the condensed particles and another 
medium, and change in amount of degradation of the compound 
in the condensed particles; the differential equation in 
the high-concentration air is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the high-concentration air, amount of 
discharge of the compound, amount of transference of the 
compound between the high-concentration air and another 
medium, and change in amount of degradation of the compound 
in the high-concentration air; the differential equation 
in the medium-concentration air is a differential equation 
stating a relationship among temporal change of fugacity 
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of the compound in the medium-concentration air, amount of 
transference of the compound between the medium- 
concentration air and another medium, and change in amount 
of degradation of the compound in the medium-concentration 
air; the differential equation in the low-concentration air 
is a differential equation stating a relationship among 
temporal change of fugacity of the compound in the low- 
concentration air, amount of discharge of the compound 
outdoors, amount of transference of the compound between 
the low-concentration air and another medium, and change 
in amount of degradation of the compound in the low- 
concentration air; the differential equation at the floor 
is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the floor, 
temporal change in volume of the floor, amount of 
transference of the compound between the floor and another 
medium, and change in amount of degradation of the compound 
at the floor; the differential equation at the wall is a 
differential equation stating a relationship among temporal 
change of fugacity of the compound at the wall, temporal 
change in volume of the wall, amount of transference of the 
compound between the wall and another medium, and change 
in amount of degradation of the compound at the wall; and 
the differential equation at the ceiling is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the ceiling, temporal change 
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in volume of the ceiling, amount of transference of the 
compound between the ceiling and another medium, and change 
in amount of degradation of the compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is sprayed over the whole floor; whereas the 
above-mentioned media are suspended particles which are 
divided into at least one kind according to size, indoor 
air, a floor, a wall, and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the suspended particles, temporal change 
in volume of the suspended particles, amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 
in the suspended particles; the differential equation in 
the indoor air is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 
amount of degradation of the compound in the indoor air; 
the differential equation at the floor is a differential 
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equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
particles to the floor, amount of transference of the 
compound between the floor and another medium, and change 
in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound at the wall, temporal change in volume of 
the wall, amount of attachment of the suspended particles 
to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall; and the differential 
equation at the ceiling is a differential equation stating 
a relationship among temporal change of fugacity of the 
compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 
of degradation of the compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention may be such that the floor is constituted by a 
rug having ears of fiber, whereas a space between the ears 
is added to the above-mentioned media. 

Preferably, in this case, the differential equation 
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in the space between the ears is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the space between the ears, temporal 
change in volume of the solution containing the compound 
in the space between the ears, amount of attachment of the 
compound into the space portion between the ears by falling, 
amount of transference of the compound between the space 
portion between the ears and another medium, and change in 
amount of degradation of the compound in the space portion 
between the ears . 

Even in the case where a rug having ears of fiber is 
spread on the floor, when the space between the ears is 
further added to the media, the behavior of the compound 
can be simulated accurately, thus allowing various kinds 
of simulations to be performed. 

In order to achieve the above-mentioned object, the 
computer program product of the present invention is a 
computer program product to be used together with an 
information processing apparatus comprising input means for 
receiving a data input from outside, display means, and 
readout means for reading out information from a 
computer-usable storage medium; the computer program 
product comprising a computer-usable storage medium which 
has a program area for storing a program and has a 
computer-readable program materialized in the storage 
medium for causing, according to data input from the input 
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means, the display means to display a result of simulation 
of an indoor behavior of a pesticidal compound; the computer 
program product comprising, in the program area, a program 
for dividing an indoor environment into predetermined media 
and forming a differential equation concerning a fugacity 
of the compound, a program for determining the fugacity of 
the compound in each of the media from the differential 
equation; a program for determining the indoor behavior of 
the compound from the fugacity of the compound in each of 
the media, and a program for changing, in response to a 
fluctuation in mass balance of the compound indoors , a minute 
time unit used when solving the differential equation. 

As the indoor environment is divided into 
predetermined media, and exchanges of the chemical compound 
between the media and the like are taken into account, 
simulation results close to the actual behavior of the 
compound can be obtained, while the minute time unit can 
be set automatically in response to fluctuation in mass 
balance when solving simultaneous differential equations 
including a parameter which changes over time . Accordingly , 
when a computer processes the above-mentioned differential 
equation, accurate solutions can be obtained in a short time. 

Preferably, the computer program product further 
comprises, in the program area, a program for evaluating 
safety of the compound with respect to a human body according 
to the indoor behavior of the compound. 
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As a consequence of this configuration, the safety of 
the pesticidal compound with respect to the human body can 
be evaluated accurately in a short time. Accordingly, when 
formulating a chemical such as insecticide including the 
5 above-mentioned compound, simulation can be easily repeated 

while changing conditions, thereby making it easier to 
formulate a chemical having a high safety conforming to the 
aimed object. 

The computer program product of the present invention 

10 may be such that the above-mentioned compound is introduced 

into an indoor space as a solution containing the compound 
is residually sprayed; whereas the above-mentioned media 
are a spraying site, suspended particles which are divided 
into at least one kind according to size, indoor air, a floor, 

15 a wall, and a ceiling. 

Preferably, in this case, the differential equation 
at the spraying site is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound at the spraying site, temporal change in volume 

2 0 of the spraying site, amount of attachment of the suspended 

particles to the spraying site, amount of transference of 
the compound between the spraying site and another medium, 
and change in amount of degradation of the compound at the 
spraying site; the differential equation in the suspended 

25 particles is a differential equation stating a relationship 

among temporal change of fugacity of the compound in the 
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suspended particles, temporal change in volume of the 
suspended particles, amount of transference of the compound 
between the suspended particles and another medium, and 
change in amount of degradation of the compound in the 
suspended particles; the differential equation in the 
indoor air is a differential equation stating a relationship 
among temporal change of fugacity of the compound in the 
indoor air, amount of discharge of the compound outdoors, 
amount of transference of the compound between the indoor 
air and another medium, and change in amount of degradation 
of the compound in the indoor air; the differential equation 
at the floor is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound at the floor, temporal change in volume of the floor, 
amount of attachment of the suspended particles to the floor, 
amount of transference of the compound between the floor 
and another medium, and change in amount of degradation of 
the compound at the floor; the differential equation at the 
wall is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the wall, 
temporal change in volume of the wall, amount of attachment 
of the suspended particles to the wall, amount of 
transference of the compound between the wall and another 
medium, and change in amount of degradation of the compound 
at the wall; and the differential equation at the ceiling 
is a differential equation stating a relationship among 
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temporal change of fugacity of the compound at the ceiling, 
temporal change in volume of the ceiling, amount of 
attachment of the suspended particles to the ceiling, amount 
of transference of the compound between the ceiling and 
5 another medium, and change in amount of degradation of the 

compound at the ceiling. 

The computer program product of the present invention 
may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 

10 is spatially sprayed; whereas the above-mentioned media are 

suspended particles which are divided into at least one kind 
according to size, indoor air, a floor, a wall, and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 

15 stating a relationship among temporal change of fugacity 

of the compound in the suspended particles, temporal change 
in volume of the suspended particles, amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 

20 in the suspended particles; the differential equation in 

the indoor air is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
compound outdoors, amount of transference of the compound 

25 between the indoor air and another medium, and change in 

amount of degradation of the compound in the indoor air; 
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the differential equation at the floor is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
5 particles to the floor, amount of transference of the 

compound between the floor and another medium, and change 
in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 

10 of the compound at the wall, temporal change in volume of 

the wall, amount of attachment of the suspended particles 
to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 

15 equation at the ceiling is a differential equation stating 

a relationship among temporal change of fugacity of the 
compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 

2 0 between the ceiling and another medium, and change in amount 

of degradation of the compound at the ceiling. 

The computer program product of the present invention 
may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 

25 is heated to vaporize; whereas the above-mentioned media 

are condensed particles which are divided into at least one 



19 



SCC 97-06 



kind according to generation and extinction, high- 
concentration air, medium-concentration air, low- 
concentration air, a floor, a wall, and a ceiling which is 
divided into at least one kind according to compound 
5 concentration. 

Preferably, in this case, the differential equation 
in the condensed particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the condensed particles, temporal change 

10 in volume of the condensed particles , amount of transference 

of the compound between the condensed particles and another 
medium, and change in amount of degradation of the compound 
in the condensed particles; the differential equation in 
the high-concentration air is a differential equation 

15 stating a relationship among temporal change of fugacity 

of the compound in the high-concentration air, amount of 
discharge of the compound, amount of transference of the 
compound between the high-concentration air and another 
medium, and change in amount of degradation of the compound 

20 in the high-concentration air; the differential equation 

in the medium-concentration air is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the medium-concentration air, amount of 
transference of the compound between the medium- 

25 concentration air and another medium, and change in amount 

of degradation of the compound in the medium-concentration 
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air; the differential equation in the low-concentration air 
is a differential equation stating a relationship among 
temporal change of fugacity of the compound in the low- 
concentration air, amount of discharge of the compound 
5 outdoors, amount of transference of the compound between 

the low-concentration air and another medium, and change 
in amount of degradation of the compound in the low- 
concentration air; the differential equation at the floor 
is a differential equation stating a relationship among 

10 temporal change of fugacity of the compound at the floor, 

temporal change in volume of the floor, amount of 
transference of the compound between the floor and another 
medium, and change in amount of degradation of the compound 
at the floor; the differential equation at the wall is a 

15 differential equation stating a relationship among temporal 

change of fugacity of the compound at the wall, temporal 
change in volume of the wall, amount of transference of the 
compound between the wall and another medium, and change 
in amount of degradation of the compound at the wall; and 

20 the differential equation at the ceiling is a differential 

equation stating a relationship among temporal change of 
fugacity of the compound at the ceiling, temporal change 
in volume of the ceiling, amount of transference of the 
compound between the ceiling and another medium, and change 

25 in amount of degradation of the compound at the ceiling. 

The computer program product of the present invention 
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may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 
is sprayed over the whole floor; while the above-mentioned 
media are suspended particles which are divided into at least 
5 one kind according to size, indoor air, a floor, a wall, 

and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
stating a relationship among temporal change of fugacity 

10 of the compound in the suspended particles, temporal change 

in volume of the suspended particles, amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 
in the suspended particles; the differential equation in 

15 the indoor air is a differential equation stating a 

relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 

2 0 amount of degradation of the compound in the indoor air; 

the differential equation at the floor is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 

25 particles to the floor, amount of transference of the 

compound between the floor and another medium, and change 
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in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound at the wall, temporal change in volume of 
the wall, amount of attachment of the suspended particles 
to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 
equation at the ceiling is a differential equation stating 
a relationship among temporal change of fugacity of the 
compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 
of degradation of the compound at the ceiling. 

The computer program product of the present invention 
may be such that the floor is constituted by a rug having 
ears of fiber, whereas a space between the ears is added 
to the above-mentioned media. 

Preferably, in this case, the differential equation 
in the space between the ears is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the space between the ears, temporal 
change in volume of the solution containing the compound 
in the space between the ears , amount of attachment of the 
compound into the space portion between the ears by falling, 
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amount of transference of the compound between the space 
portion between the ears and another medium, and change in 
amount of degradation of the compound in the space portion 
between the ears . 
5 Even in the case where a rug having ears of fiber is 

spread on the floor, when the space between the ears is 
further added to the media, the behavior of the compound 
can be simulated accurately, thus allowing various kinds 
of simulations to be performed. 

10 Brief Description of Drawings 

Fig. 1 is a flowchart showing a method of simulating 
an indoor behavior of a pesticidal compound in accordance 
with an embodiment of the present invention; 

Fig. 2 is an explanatory view showing a model of indoor 

15 environment in the case where an indoor behavior of a 

compound is performed in the first embodiment of the present 
invention; 

Fig. 3A is a graph showing temporal change of 
concentration in indoor air; 
20 Fig. 3B is an explanatory view in the case where mean 

concentration in indoor air is determined by use of the graph 
in Fig. 3A; 

Fig. 4A is a graph showing temporal change in residual 
amount ; 

25 Fig. 4B is an explanatory view in the case where mean 

residual amount is determined by use of the graph in Fig. 
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4A; 

Fig. 5 is a flowchart showing a method of setting an 
estimation nick time width; 

Fig. 6A is an explanatory view showing a modeled state 
5 of a carpet in the case where an indoor behavior of a 

pesticidal compound is estimated in accordance with the 
second embodiment of the present invention; 

Fig. 6B is an enlarged view of a carpet fiber in the 
case where an indoor behavior of a pesticidal compound is 
10 estimated in accordance with the second embodiment of the 

present invention; 

Fig. 6C is an enlarged view of a space portion of the 
carpet in the case where an indoor behavior of a pesticidal 
compound is estimated in accordance with the second 
15 embodiment of the present invention; 

Fig. 7 is a graph showing temporal change of 
concentration in the air in Example 1 in the above-mentioned 
first embodiment; 

Fig. 8A is a graph showing temporal change in 
2 0 concentration in the air in the case where ventilation rate 

is changed in Example 2; 

Fig. 8B is a graph showing temporal change in 
concentration in the air in the case where a window is opened 
and closed in Example 3 ; 
25 Fig. 9A is a graph showing temporal change in residual 

amount on the floor in Example 1 ; 
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Fig. 9B is a graph showing temporal change in residual 
amount on the wall and ceiling in Example 1; 

Fig. 1 OA is a graph comparing actually measured values 
of concentration in the air and their simulation results 
5 in Example 1 ; 

Fig. 1 OB is a graph comparing actually measured values 
of residual amount at the spraying site, floor, wall, and 
ceiling and their simulation results in Example 1; 

Fig. 11 is an explanatory view showing a modeled state 
10 of indoor environment in the case where an indoor behavior 

of a compound is estimated in the third embodiment of the 
present invention; 

Fig. 12A is a graph comparing actually measured values 
of concentration in the air and their simulation results 
15 in Examples 6 to 10; 

Fig. 12B is a graph comparing actually measured values 
of concentration in the air and their simulation results 
in Examples 6 to 10; 

Fig. 12C is a graph comparing actually measured values 
20 of residual amount on the floor and ceiling and their 

simulation results in Example 6; 

Fig. 13A is an explanatory view showing a modeled state 
of a carpet in the case where an indoor behavior of a 
pesticidal compound is estimated in accordance with the 
25 fourth embodiment of the present invention; 

Fig. 13B is an enlarged view of a carpet fiber in the 
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case where an indoor behavior of a pesticidal compound is 
estimated in accordance with the fourth embodiment of the 
present invention; 

Fig. 13C is an enlarged view of a space portion of the 
5 carpet in the case where an indoor behavior of a pesticidal 

compound is estimated in accordance with the fourth 
embodiment of the present invention; 

Fig. 14 is an explanatory view showing an example of 
size of a typical room, together with an airflow; 
10 Fig. 15 is an explanatory view showing a modeled state 

of indoor environment in the case where an indoor behavior 
of a compound is estimated in the fifth embodiment of the 
present invention; 

Fig. 16A is an explanatory view showing a modeled state 
15 of indoor environment in the case where an indoor behavior 

of the compound is estimated in the fifth embodiment of the 
present invention; 

Fig. 16B is an enlarged view of an electric heating 
evaporator in the fifth embodiment of the present invention; 
2 0 Fig. 17 is an explanatory view showing a modeled state 

of indoor environment in the case where an indoor behavior 
of a compound is estimated in the seventh embodiment of the 
present invention; 

Fig. 18 is an explanatory view showing a falling state 
25 of a particle zone; 

Fig. 19 is a conf igurational view of a storage medium 
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20 employed in each of the above-mentioned embodiments of 
the present invention; 

Fig. 20 is a conf igurational view of an information 
processing apparatus for operating a program stored in the 
5 storage medium 20; and 

Fig. 21 is a conf igurational view of the information 
processing apparatus for operating a program stored in the 
storage medium 20. 

g est M odes fo r C arrying Out the invention 

10 (1) First Embodiment 

With reference to the drawings , the first embodiment 
of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a 
pesticidal compound (hereinafter simply referred to as 

15 compound) contained in a chemical (solution) such as 

insecticide in accordance with this embodiment estimates, 
for example, the behavior of the compound when the chemical 
is residually sprayed in order to exterminate indoor vermin 
such as cockroach. Here, residual spraying refers to a 

2 0 process in which a chemical is sprayed to a local site such 

as the boundary between a floor and a wall . The chemical 
encompasses all kinds of insecticides and the like, 
including pyrethroid insecticidal compounds, 
organophosphorus compounds, carbamate compounds, and 

25 insect growth restrainers (IGR). 

First, as shown in Fig. 2, an environment is classified 
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into seven kinds of media consisting of spraying site 1, 
large particles (suspended particles) 2, small particles 
(suspended particles) 3, whole indoor air 4, floor 5, wall 
6, and ceiling 7. 

When a chemical is sprayed to the boundary between the 
floor 5 and the wall 6, part of the chemical attaches to 
the spraying site 1, whereas the rest dissipates into the 
air as large and small particles 2 and 3. The chemical 
attached to the spraying site 1 not only infiltrates into 
this site but also evaporates. On the other hand, among the 
large and small particles 2 and 3 , the smaller particles 
3, due to their weight, scatter higher than the larger 
particles 2 . When these particles attach onto interior 
materials of the floor 5, wall 6, and ceiling 7, the compound 
in the particles moves into the air 4 after infiltrating 
into the interior materials or is degraded by light or 
oxidation. Such large and small particles 2 and 3 and 
evaporated chemical have decreased due to ventilation of 
the room. 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned seven kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is residually sprayed, i.e., compound 
distribution in the room, and temporal change in amount of 
attachment to the spraying site 1 , floor 5 , wall 6 , or ceiling 
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7 are estimated and analyzed. 

Here, since the floor 5 mainly comprises three kinds, 
i.e., tatami, flooring, or carpet (rug), the differential 
equation is formed so as to correspond to the kind of floor. 
5 Explained in this embodiment is the case where the floor 

5 is made of flooring. 

With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 

10 be mainly divided into a step of dividing an indoor 

environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2) ; a step of determining the fugacity 
of the compound in each medium from the differential equation 

15 (S3); a step of determining the indoor behavior of the 

compound from the fugacity of the compound in each medium 
(S4); a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 

20 see Fig. 5) ; and a step of evaluating safety of the compound 

with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step ( hereinafter abridged as S ) 1, a primary 
condition is inputted. The primary condition comprises 

2 5 physicochemical properties of the compound (see Table 1), 

indoor environmental behavior properties of the compound 
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(see Table 2 ) , indoor environment (see Table 3) , and product 
properties (see Table 4). 
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TABLE 1 



PHYSICOCHEMICAL 
PROPERTY OF COMPOUND 


SYMBOL 


EXAMPLE OF VALUE 


MOLECULAR WEIGHT 




350.46 Cg/mole) 


SPECIFIC GRAVITY 


P 


1.11 Cg/cm3] 


VAPOR PRESSURE 


ps 


1.60 XI 0" 4 [Pa] 


WATER SOLUBILITY 


CS 


6.28 XI 0" 3 [mole/m3] 


OCTANOL/ WATER 
DISTRIBUTION COEFFICIENT 


Kow 


7.59 X10 4 



TABLE 2 



INDOOR ENVIRONMENTAL 

BEHAVIOR 
PROPERTY OF COMPOUND 


SYMBOL 


EXAMPLE OF VALUE 


DEGRADATION HALF-LIFE 
IN SPRAYING SITE 
PARTICLE 
AIR 
FLOOR 

WALL /CEILING 


z l 

Ti <i=2,3) 

x 4 

^5 

Tj(j=6,7) 


13 [day] 
7.7 Ch] 
7.7 Ch] 

13 Cday] (IN CASE OF FLOORING) 
31 [day] (IN CASE OF WALLPAPER) 



EXAMPLE OF VALUE | 


5 s 
§ 1 s 
§ 3 1 

O ^ ft 
^ ^ O O O 

« W W W 

1 2 1 % I § _ 2 £ £ £ 

<r> o\ o ^ ^ °. ^ 
tN <N VD <— i — < O OOOOO 


& 

g § 


SYMBOL | 


j> H -9- JE O 8 | 


j 8 


INDOOR ENVIRONMENT 


ROOM SIZE 
TEMPERATURE 
RELATIVE HUMIDITY 
ABSOLUTE HUMIDITY 

VENTILATION RATE 

INDOOR VAPOR PRESSURE IN 
PARTICLE - CONSTITUTING SOLVENT 

AIR DIFFUSION COEFFICIENT 

OIL COMPONENT CONTENT 


WET - BULB TEMPERATURE 
VAPOR PRESSURE OF WATER 


OIL BASE 


WATER BASE 
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Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2) . The 
secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 5) determined by the 
5 primary condition alone, and a temporal change condition 

accompanying a temporal change (see Table 6). Their 
calculations will be explained later in detail. 
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TABLE 6 



TEMPORAL CHANGE CONDITION 


SYMBOL 


SUSPENDED PARTICLE 

DIAMETER 
FALLING SPEED 
FLOATING NUMBER 
FLOATING HEIGHT 


di 0=2,3) 
vi (i=2,3) 
ni Q=2,3) 
Hi 0=2,3) 


TRANSFERENCE SPEED OF COMPOUND 
IN SPRAYING SITE 

SUSPENDED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


kl 

ki 0=2,3) 
k 4 

kj 0=5,6,7) 


TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN SPRAYING SITE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN SUSPENDED PARTICLE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN AIR AND FLOOR / WALL / CEILING 


Dm 

Di4 0=2,3) 
D 4 j (3=5,6,7) 


VOLUME OF SPRAYING SITE 
VOLUME OF SUSPENDED PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


Vl 

Vi (i=2,3) 
Vj 0=5,6,7) 
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Using the secondary condition, seven kinds of Fugacity 
are calculated (S3). Namely, seven kinds of differential 
equations concerning the spraying site 1, large and small 
particles 2 and 3, air 4, floor 5, wall 6, and ceiling 7 
are simultaneously formed and are solved by Runge- 
Kutta-Gill method, whereby seven kinds of Fugacity are 
computed over time. Here, an estimation nick time width 
(minute time unit) set when solving the differential 
equations is automatically set so as to be varied in response 
to a fluctuation in mass balance of the compound. 

At S4, using thus computed Fugacity (f 3 , f 4 ) of the 
small particles 3 and air 4, together with Fugacity (f 2 ) of 
the large particles 2 when necessary, temporal 
concentration of the compound in the indoor air is computed; 
whereas temporal residual amounts of the compound to the 
spraying site 1 and floor 5 are computed by use of Fugacity 
(f x , f 5 ) of the spraying site 1 and floor 5. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 
an estimated exposure amount in inhalation indicating a 
degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 
estimated exposure amount in inhalation (S7). At S8, the 
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inhalation safety coefficient is compared with a reference 
value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to Si, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the spraying site 1 or floor 5 to which the 
chemical is attached (S9) . Thereafter, at S10, it is judged 
whether to perform a percutaneous safety evaluation or not. 
If the percutaneous safety evaluation is to be performed, 
a percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 
condition is considered. 

If the percutaneous safety evaluation is not selected 
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at S10, calculated according to the estimated percutaneous 
exposure amount is an estimated oral exposure amount 
indicating the degree of exposure in the case where the 
chemical attached to a hand or the like is taken orally (S13 ) . 
Subsequently, according to the estimated oral exposure 
amount, an oral safety coefficient is computed (S14) . This 
oral exposure may occur, in particular, when an infant puts 
a chemical-attached hand into the mouth. At S15, as with 
the safety evaluation in inhalation, the oral safety 
coefficient is compared with a reference value defined in 
each country, whereby the safety is evaluated. If there is 
a problem in safety at S15, the operation returns to SI, 
where alteration of the primary condition is considered. 

Finally, if it is judged to be safe at S8, S12, or S15, 
overall safety is evaluated. Here, the sum of respective 
reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
(S16). This overall safety coefficient is compared with a 
reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4 , S6 , S7 , S9 , Sll , S13 , and S14 will be explained in detail . 

(i) Secondary Condition Calculation (S2) 

Calculations of the definite factor condition shown 
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in Table 5 and the temporal change condition shown in Table 

6 will be explained. 

According to surface temperature (T d ) of particles, 
properties of particle-constituting solvent (vapor 
pressure p d , molecular weight M d , and specific gravity p 
d ), and indoor environment (air diffusion coefficient D air , 
vapor pressure poo, and room temperature Too), the 
evaporation constant ( a ) of the large and small particles 

2 and 3 is defined as follows: 

wherein R is a gas constant. 

Here, the upper parts of the above-mentioned Tables 

3 and 4 exemplify the case where the residually sprayed 
aerosol is of oil base. In this case, Too (temperature at 
the site far from particles) and T d are set to room 
temperature (T), whereas poo (vapor pressure of oil at the 
site far from particles) is set to zero. 

In the case where the residually sprayed aerosol is 
of water base, it is necessary for P d , M d , p d , T d , and P 
oo in expression (1) to be changed to values based on 
properties of water. Specifically, the above-mentioned 
factor values when the room temperature (T) is 298 [K] (= 
25 °C) and the relative humidity is 60[%RH] are exemplified 
in the lower parts of Tables 3 and 4 . 

Here, T d can be determined, according to the room 
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temperature (T) and relative humidity (H) , from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 
Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
the other hand, P d and poo can be computed from the following 
expressions : 

logioPa = 10.23 - 1750/(T d -38) 
log lo poo = 8.23 - 1750/(T-38) + log lo 0 
wherein ip is an indoor relative humidity. 

From the vapor pressure (P 3 ), water solubility (C 3 ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and from the oil component content (p x ) of the 
floor or wall constituting the spraying site 1, Fugacity 
capacity (ZJ of the spraying site 1 is expressed as: 
Z^p^C'/P 5 ...(2) 

From the vapor pressure (P 3 ) of the compound and the 
surf ace temperature (T d ) of particles, Fugacity capacity (Z LI 
i = 2 , 3 ) of the large and small particles 2 and 3 is expressed 
as : 

Z,--^ ~<3) 

Pf RT d v 7 

Pi = P s exp{6.79(r w IT d -l)} 

Though the vapor pressure of the compound in a liquid state 
(P L 3 ) is computed by use of the melting point T M of the harmful 
material and T d here; in the case where T M cannot be obtained, 
P L S may be set identical to the vapor pressure of the compound 
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in a solid state (P s ). 

From the room temperature (T), Fugacity capacity (Z 4 ) 
of the air 4 is expressed as follows: 

Z 4 =— ...(4) 
RT w 

From the vapor pressure (P s ), water solubility (C s ), 

and octanol/water distribution coefficient (K ow ) of the 

compound, and from the particle-constituting solvent 

component content (p jr j = 5, 6, 7) of the material 

constituting the floor 5, wall 6, and ceiling 7, Fugacity 

capacity ( Zj , j = 5 , 6 , 7 ) of the floor 5 , wall 6 , and ceiling 

7 is expressed as: 

Zj -PjK^C'/I* ...(5) 

Though the ratio ( 7 j) of the large and small particles 

2 and 3 attaching to the floor 5, wall 6, and ceiling 7 may 
be assumed to be such that 100% of them attach to the floor 
5 from viewpoint of safety evaluation, it may be set, for 
example, floor : wall : ceiling = 97:2:1 in view of the 
actually measured values in the past. 

The ratio (7i) of the large and small particles 2 and 

3 falling to the spraying site 1 is determined by the 
above-mentioned ratio (7j), and the area (AJ occupied by 
the spraying site 1 with respect to the size of the room 
(V 4 ). 

The diffusion coefficient (D c ) of the compound in the 
floor 5, wall 6, and ceiling 7 can be computed as follows. 
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Namely, D c can be determined when the diffusion coefficient 

(D ca ) of the compound in the air is multiplied by 10" 9 , whereas 

D oa can be computed from properties of the compound (e.g., 

structural formula, molecular weight, and the like) 

according to Wike and Lee method ("Handbook of Chemical 

Property Estimation Methods," McGraw-Hill Book Company, 

1982), for example. 

Mainly generated in the seven kinds of media are 

degradation reactions caused by light and oxidation. The 

degradation constant (K x ) of the spraying site 1, 

degradation constant (K if i = 2, 3) of the large and small 

particles 2 and 3, degradation constant (K 4 ) of the air 4, 

degradation constant (K s ) of the floor 5, and degradation 

constant ( Kj , j = 6 , 7 ) of the wall 6 and ceiling 7 are defined 

by use of degradation half-life (T 1 ), (T^), (T 4 ), (T 5 ), 

or (Tj) as: 

^ = 0.693/^ 
K t =0.693 jr. 

K 4 = 0.693 /r 4 j, ...(6) 
K 5 = 0.693 It. 
Kj =0.693 fx 

Here, among the degradation half -lives shown in the 
above-mentioned Table 2, it is difficult to actually measure 
the degradation half-life (T 4 ) of the compound in the air, 
whereby it is not usually easy for the degradation constant 
(K 4 ) to be determined. Therefore, when the actually 
measured value is unavailable, it may be determined as being 
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computed by Atmospheric Oxidation program (Atkinson et al., 
1984, Chem. Rev. Vol. 84, pp. 437-470) using the structural 
formula of the compound. Also, the degradation constant in 
the particles and in the interior materials (floor 5, wall 
5 6, and ceiling 7) may be set to zero in view of the safety 

of inhabitants when the actually measured values are 
unavailable. 

Since the volume of particles decreases over time due 
to evaporation of the particle-constituting solvent 

10 component, during the evaporation of particle-constituting 

solvent component, using the diameter (d oi , i = 2, 3) of 
particles immediately after spraying (t = 0) and the 
evaporation constant (a) determined from expression (1), 
the diameter (di, i = 2, 3) of the large and small particles 

15 2 and 3 at time t is represented as shown in the following 

expression ( 7 ) : 
d.^^jdl-lat ...{!) 
4= 3 V*T 4* -(7') 

By contrast, the diameter (d £ ) of particles after 
complete evaporation of the particle-constituting solvent 

20 component is represented as shown in the above expression 

(7' ) by use of the ratio of volume of the compound in the 
product (R a ) and the diameter (d oi ) of particles immediately 
after spraying. 

From the compound content (C a ) and specific gravity 

25 ( p ) of the compound, the ratio of volume of the compound 
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in the product (R a ) is expressed as Cjp. 

Since movement of the particles is under the control 
of gravity and air resistance, according to Stokes ' law using 
the specific gravity ( p d ) of the particle-constituting 
5 solvent and the diameter (dj of particles determined by 

expression (7), the falling speed (Vi, i = 2, 3) of the large 
and small particles 2 and 3 before complete evaporation of 
the particle-constituting solvent component is represented 
by the following expression (8): 

v,-&&-d?-fi(dl-2at) -(8) 

10 S =1 + ^— -{6.32+ 2.01 expf- 8.322 xlOV)} ...(8') 

v t ~{p/ Pd )fidf ...(8") 

wherein g is gravitational acceleration, V is the 
coefficient of viscosity of the air 4, S c is a sliding 
correction coefficient , and /? is a speed coefficient . Here, 
while particles do not conform to the Stokes * law when they 

15 become small, the sliding correction coefficient (S c ) , which 

is a coefficient for correcting this phenomenon, is 
represented by the above expression (8') according to the 
diameter (dj of particles determined by expression (7). 
On the other hand, the falling speed (vj of the 

2 0 particles after the complete evaporation of particle- 

constituting solvent component is represented by the above 
expression (8") using the specific gravity (p) of the 
compound, specific gravity ( p d ) of the particle- 
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constituting solvent, speed coefficient ( /3 ) , and diameter 
(di) of particles determined by expression (7"). 

From transference speed (k 4 ) of the compound in the 
air 4, which will be explained later, floating height of 
5 particles (Hi), size of the room (V 4 = L 4 * W 4 x height), and 

the falling speed (Vi) of particles determined by expression 
( 8 ) or ( 8 " ) , the number of floating large and small particles 
2 and 3 (n L , i = 2, 3) is expressed as: 

— - = -v. -k, ...(9) 

dt H i L 4 W 

10 Here, the floating height (H^ of particles refers to the 

maximum height of the particle distribution from the floor 
5 after the particles are reflected by the spraying site 

1 and then are suspended in the case where the chemical is 
sprayed to the spraying site 1. It is positioned higher as 

15 the particles are finer. 

Using the diffusion coefficient (D e ) of compound in 
the floor 5, wall 6, and ceiling 7, the transference speed 
(k x ) of the compound in the spraying site 1 is expressed as 
(D c /t) 0 - 5 . 

20 From the transference speed (k 4 ) of the compound in 

the air 4 indicated in the following, the transference speed 
( k ir i = 2 , 3 ) of the compound in the large and small particles 

2 and 3 is expressed as k 4 /100. 

From the ventilation rate (G) , the transference speed 
25 (k„ ) of the compound in the air 4 is expressed as GV 4 /A 4 . Here, 
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the above-mentioned A 4 is the cross-sectional area of the 
room with respect to the moving direction of the air 4 
(arrowed direction in Fig. 2). 

Using the diffusion coefficient <D C ) of compound in 
5 the floor 5, wall 6, and ceiling 7 , the transference speed 

(kj, j = 5, 6, 7) of the compound in the floor 5, wall 6, 
and ceiling 7 is expressed as (D c /t) 0-5 . 

From the transference coefficient (k x ) of the compound 
in the spraying site 1, transference speed (k 4 ) of the 

10 compound in the air 4, area (A t ) of the spraying site 1, 

falling speed (v d ) of particles determined by expression (8) 
or (8"), Fugacity capacity (Z 2 ) of the spraying site 1 
determined by expression (2), and Fugacity capacity (Z 4 ) of 
the air 4 determined by expression (4), the transference 

15 coefficient (D 14 ) of the compound between the spraying site 

1 and the air 4 is represented as follows: 

D = 1 ...(10) 

14 l/^A^ZJ + l/^A^ZJ 

From the transference speed (k t ) of the compound in 
the large and small particles 2 and 3, transference speed 

2 0 (k 4 ) of the compound in the air 4, diameter (di) of particles 

determined by expression (7) or (7' ), Fugacity capacity (ZJ 
of particles determined by expression (3), and Fugacity 
capacity (Z 4 ) of the air 4 determined by expression (4), the 
transference coefficient (D i4 , i = 2, 3) between the large 

25 and small particles 2 and 3 and the air 4 is represented 
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D. 4 = —. r -.(11) 

;4 l/(A t .A ; Z ; ) + l/{(/r 4 + v ; )A ; Z 4 } 

wherein K L (= 7Cd*) is the surface area of particles. 

From the transference speed (k 4 ) of the compound in 
5 the air 4, transference speed (kj) of the compound in the 

floor 5, wall 6, and ceiling 7, surface area (A 3 ) of the floor 
5, wall 6, and ceiling 7, Fugacity capacity (Z 4 ) of the air 
4 determined by expression (4), and Fugacity capacity (Z 3 ) 
of the floor 5 , wall 6 , and ceiling 7 determined by expression 
10 (5), the transference coefficient (D 4j , j = 5, 6, 7) of the 

compound between the air 4 and the floor 5 , wall 6 , and 
ceiling 7 becomes : 

£>..= ...(12) 

" 1/(& 4 A.Z 4 ) + 1/(£ ; A.Z.) 

Here, the surface area (A t ) of each of the floor 5, wall 6, 
15 and ceiling 7 is determined by the size of the room (V 4 = 

L 4 x w 4 x height) and the area (A t ) of the spraying site 1. 

The volume (V\) of the spraying site 1 is determined 
by the diffusion coefficient (D c ) in the floor 5, wall 6, 
and ceiling 7. 

20 Using the diameter (d^ of suspended particles 

determined by expression (7) or (7') and the evaporation 
constant (a) determined by expression (1), the volume 
change (dv i /dt) of the large and small particles 2 and 3 is 
represented as: 
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Namely, as the oil component or moisture in the large and 
small particles 2 and 3 evaporates over time, their volume 
(Vi, i = 2, 3) decreases. 
5 Assuming that each of the floor 5, wall 6, and ceiling 

7 before residually spraying the chemical is like a thin 
film and that, as the chemical infiltrates into the film, 
its thickness increases so as to enhance the volume thereof; 
from the surface area (Aj) of the floor 5, wall 6, and ceiling 

10 7 and diffusion coefficient (D c ) of the compound in the floor 

5, wall 6, and ceiling 7, the volume (V 3 , j = 5, 6, 7) of 
floor 5, wall 6, and ceiling 7 is indicated as follows: 
V. = 2^rTtA i ...(14) 

(ii) Fugacity Calculation (S3) 

15 The behavior of the compound in the spraying site 1 

is expressed in the form of differential eguation concerning 
Fugacity (fj as: 

—V l Z l = -jDJtA.ZJt V- change 

dt 

F- change +Yi^ n i V i Z ifi Deposition^) ...(15) 

- D u ( / t - / 4 ) Transfer en ce(V) 

-Ky i ZJ l Degradation 

Here, the terms of V-change, Deposition (i), Transference 
2 0 (1)/ anci Degradation respectively indicate volume change 

{increase over time) of the spraying site 1, attachment 
accompanying the falling of the large and small particles 
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2 and 3, amount of transference of the compound between the 
spraying site 1 and the air 4, and change in amount of 
photodegradation of the compound. 

The behavior of the compound in the large and small 
particles 2 and 3 is expressed in the form of differential 
equation concerning Fugacity (f it i = 2, 3) as: 

Q-V.Z. = -ad.Z.f. V -change 

dt 1 1 2 ' 

F -change -Du{fi~U) Transference^ 

- D i4 ( / ; - f 4 ) Transference^) 

- KyZi f. Degradation 

Here, the terms of V-change, Transference (1), Transference 
(4), and Degradation respectively indicate volume change 
(decrease over time) of the large and small particles 2 and 
3, amount of transference of the compound between the 
spraying site 1 and the air 4, amount of transference of 
the compound between the large and small particles 2 and 

3 and the air 4 , and change in amount of photodegradation 
of the compound. 

The behavior of the compound in the air 4 is expressed 
in the form of differential equation concerning Fugacity 
(f 4 ) as: 

^V A Z, = -GV.Z.f. Ventilation 

dt 4 4 4 4 4 

F - change - D l4 (/ 4 - f ) Transference^ 



-J? n ; D i4 (f 4 -f) Transference^) ...(17) 

- ^ D 4j (f 4 - /•) Transference^]) 

- K 4 V 4 Z 4 f 4 Degradation 
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Here, the terms of Ventilation, Transference (1), 
Transference (i), Transference (j), and Degradation 
respectively indicate amount of discharge of the compound 
outdoors, amount of transference of the compound between 
5 the air 4 and the spraying site 1, amount of transference 
of the compound between the air 4 and the large and small 
particles 2 and 3, amount of transference between the air 
4 and the floor 5, wall 6, and ceiling 7, and change in amount 
of photodegradation of the compound. 
10 The behavior of the compound in the floor 5, wall 6, 
and ceiling 7 is expressed in the form of differential 
equation concerning Fugacity (f jf j = 5, 6, 7) as follows: 
df. . 

fa V j Z i = -^DJtA j Z ! f j V- change 

3 

F- change +Yj^ n yi Z ;f Deposition(i) ...(18) 

- D 4j (fj - / 4 ) Transfer ence(4) 

-KjVjZjfj Degradation 
Here, the terms of V-change, Deposition (i), Transference 
15 (4), and Degradation respectively indicate volume change 

(increase over time) of the floor 5, wall 6, and ceiling 
7 , attachment accompanying the falling of the large and small 
particles 2 and 3, amount of transference of the compound 
between the floor 5, wall 6, ceiling 7 and the air 4, and 
2 0 change in amount of photodegradation of the compound. 

The above-mentioned seven kinds of differential 
equations (15) to (18) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
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(fi to f 7 ). 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
width (dt) which is a minute time unit. Namely, the 
estimation nick time width is used such that solutions of 
the simultaneous differential equations are initially 
determined at a time (t 0 ), and then solutions of the 
simultaneous differential equations are determined at a 
time (t 0 + dt) to which the estimation nick time width is 
added . As solutions are obtained while estimation nick time 
widths are successively added, temporally changing Fugacity 
can be determined. Theoretically, as the set time of the 
estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 
calculation time. By contrast, when the set time is too long, 
solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 
large change. 

Specifically, mass balance is always confirmed such 
that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 
balance starts fluctuating. For example, when the 



53 



SCC 97-06 



fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 

compound input amount/ ( existing amount + degrading 
amount + discharging amount) = 0.95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is the 
compound emission amount (A a ) . Since the temporal amounts 
of compound in the seven kinds of media are determined by 
the simultaneous differential equations, they are summed 
up so as to compute the existing amount as shown in the 
following. Also, the degrading amount and discharging 
amount are as follows : 



With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
explained. 

First, an initial value of estimation nick time width 
( dt ) is inputted ( S2 1 ) . Then , an upper limit set value (e.g., 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 
is the lower limit of difference in mass balance are inputted 




...(19) 
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(S22) . Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
5 is within the range of ±5% or not (S27) . If the fluctuation 

in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 
(or t 0 ) and the mass balance at t = t + dt (or t 0 + dt) is 
at least the upper limit set value or not (S28). If it is 

10 judged to be at least the upper limit set value at S2 8, then 

solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 

15 smaller than the upper limit set value at S28, it is judged 

at S3 0 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 

2 0 the lower limit set value at S30, since solutions are not 

influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 

25 a maximum value (e.g., 0.1 [hour]) of estimation nick time 

width (dt) or not (S32). When the estimation nick time width 
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(dt) is not greater than the maximum value at S3 2, since 
solutions do not diverge, the estimation nick time width 
(dt) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S3 2, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S30, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
are continued without changing the estimation nick time 
width ( dt ) . 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S27 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
balance fluctuation, Fugacity can be computed accurately 
and efficiently. 
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(iii) Computation of Temporal Concentration in Indoor Air 
and Residual Amount (S4) 

The temporal concentration of the compound in the 
indoor air is computed when Fugacity (f 4 ) of the air, 
determined by the above-mentioned item (ii), multiplied by 
Fugacity capacity (Z„) , and Fugacity (f 3 ) of small particles 
multiplied by Fugacity capacity (Z 3 ) are summed up. Here, 
the large particles may be inhaled by a human body depending 
on the kind of chemical. In such a case, calculation is 
effected by use of Fugacity (f 2 ) and (f 3 ) of large and small 
particles . 

The residual amount of the compound is computed when 
Fugacity (f x ) of the spraying site 1 multiplied by Fugacity 
capacity (ZJ, and Fugacity (f 5 ) of the floor 5 multiplied 
by Fugacity capacity (Z 5 ) are summed up. Here, in the case 
where there is substantially no possibility of the spraying 
site 1 coming into contact with the skin, Fugacity (f 5 ) of 
the floor 5 alone may be used for calculation. 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 
indoor air forms a curve shown in Fig. 3A, for example. This 
concentration curve is integrated, an accumulated 
concentration of the compound during a specific period (t x 
to t 2 ) is determined (see Fig. 3B) , and the mean concentration 
in the indoor air is computed from thus determined value. 
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While an arbitrary period is set as the specific period 
depending on the object, an appropriate period is usually 
set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 
concentration in the indoor air, amount of respiration, and 
exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] * amount of 
respiration [m 3 /kg/min] x exposure time [min/day] 
is effected. Here, as the above-mentioned amount of 
respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 

The inhalation safety coefficient is computed from a 
non-influential amount concerning inhalation toxicity 
examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 
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inhalation safety coefficient = 

inhalation non-influential amount [mg/kg/day ] / 
estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned residual amount forms a curve such 
as that shown in Fig. 4A, for example. This residual amount 
curve is integrated, the accumulated residual amount of the 
compound during a specific period (t x to t 2 ) is determined 
(see Fig. 4B) , and the mean residual amount is computed from 
thus determined value. While an arbitrary period is set as 
the specific period depending on the object, an appropriate 
period is usually set in view of the method of use of the 
product and the test period of toxicity data. 

Then, according to the mean residual amount, skin 
attachment ratio, contact area, and body weight, the 
estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean residual amount [mg/m 2 ] x skin attachment ratio [%] 
x contact area [mVday ] ) /body weight [kg] 

is effected. Here, as the contact area, a published value 
(e.g., 4 [mVday]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor 5 where the compound 
exists. As this value, a published value or a value 
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experimentally obtained from a model may be used. 

A model experiment method for the skin attachment 
ratio is as follows. A weight (8 cm x 8 cm x 8 cm; 4.2 kg) 
is placed on a denim cloth (8 cm x 10 cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants . 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 

percutaneous non-influential amount [mg/kg/day]/ 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
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are not many published values . Accordingly, a more accurate 
value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day] /( estimated percutaneous 
exposure amount [mg/kg/day] x percutaneous absorption ratio 
[%] ) 

Here, when the percutaneous absorption ratio is unknown, 
employed is a national guideline (e.g., 10%) which usually 
exists . 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
obtained in the above-mentioned item (v) , hand surface area 
ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
estimated percutaneous exposure amount [mg/kg/day] * hand 
surface area ratio [%] * oral transference ratio [%] 
is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 
transference ratio is a hypothetical value, which is set 
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to 100%, for example. 

In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
amount from tableware or food to mouth be added to the 
estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 
the amount of the harmful material remaining in tableware, 
tableware use area which is the sum of tableware surface 
areas, and oral transference ratio from tableware, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] x tableware use area 
[mVday] x oral transference ratio [%]/body weight [kg] 
is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio ) . As the tableware contamination ratio , 
an actually measured value (e.g., 9% ) or a hypothetical value 
may be used. 

The oral safety coefficient is computed from the 
non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 
as : 
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oral safety coefficient = 

oral non-influential amount [mg/kg/day]/ estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
this embodiment, the environment is initially divided into 
the spraying site 1, large and small particles 2 and 3, whole 
indoor air 4 , floor 5 , wall 6 , and ceiling 7 ; and differential 
equations concerning Fugacity of the compound in them are 
simultaneously formed and are solved, so as to estimate the 
indoor behavior of the compound when a chemical containing 
the compound is residually sprayed. 

Since the environment is thus considered to be seven 
kinds of media, and exchanges between the individual media 
and the like are taken into account, simulation results close 
to the actual behavior of the compound can be obtained. 

The estimation nick time width set when solving the 
simultaneous differential equations is variably set, while 
constantly confirming mass balance of the compound indoors 
after the residual spraying, so that the amount of input 
of the chemical indoors and the resulting solution coincide 
with each other. 

Accordingly, since mass balance of the compound after 
the residual spraying is always monitored such that the 
amount of input of the compound indoors and the resulting 
solution coincide with each other, thereby varying the 
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estimation nick time width; the estimation nick time width 
is set longer when the mass balance fluctuates a little, 
whereas it is set shorter when the mass balance starts 
fluctuating greatly. Namely, when solving simultaneous 
differential equations including a parameter accompanying 
temporal change, the estimation nick time width is 
automatically set in response to the fluctuation in mass 
balance. Consequently, when processed by a computer, an 
accurate solution can be obtained in a short time. 

The method of evaluating safety of a pesticidal 
compound in accordance with this embodiment uses the 
estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is residually sprayed. 

Accordingly, the safety of the compound with respect 
to the human body can be evaluated accurately in a short 
time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 

Though the kind of the floor 5 is assumed to be flooring 
in this embodiment, differential equations similar to those 
mentioned above may also be formed in the case of tatami. 
In this case, however, when determining fugacity capacity 
(Z 5 ) of the floor 5 in equation (5), it is necessary for the 
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particle-constituting solvent component content (see Table 
3) to be changed from that of the flooring to that of tatami. 
Also, when determining the degradation constant (K 5 ) of the 
floor 5 in equation (6), the decomposition half-life is 
required to be changed from that of the flooring to that 
of tatami. 

Though Fugacity is determined by use of Runge- 
Kutta-Gill method in this embodiment, other methods may be 
used for solving differential equations. Runge-Kutta-Gill 
method, however, is preferably used since a program for the 
above-mentioned differential equations can be easily made 
by Basic. Also in the case where differential equations are 
solved by a method other than Runge-Kutta-Gill method, 
similar effects can be obtained when the estimation nick 
time width is set as mentioned above. 

In the following, examples (Examples 1 to 5) of the 
above-mentioned method of estimating an indoor behavior of 
a pesticidal compound and the safety evaluation method using 
the same will be explained with reference to the accompanying 
drawings . 

The evaluation of the method of estimating an indoor 
behavior of a pesticidal compound in these examples is 
effected as the results obtained by estimation and the 
actually measured values obtained by measurement of the 
actual indoor behavior of the compound are compared with 
each other. 
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First, the measurement of the actual indoor behavior 
of the compound will be explained. 

The employed aerosol can (300 [ml]) contains 0.9 [g] 
of d-phenothrin (trademark: Smithrin; C 23 H 26 0 3 : (3- 
phenoxyphenyl ) methyl ( 1R) -cis-trans-2 , 2-dimethyl-3- ( 2- 
methyl-l-prophenyl) cyclopropane carboxylate) and 1.1 [g] 
of d-tetramethrin (trademark: Neo-Pyramin Forte; C 19 H 25 0 4 : 
(1,3,4,5,6, 7-hexahydro-l , 3-dioxo-2H-isoindol-2- 
yl ) methyl ( 1R) -cis-trans-2 , 2-dimethyl-3- ( 2-methyl-l- 
prophenyl ) cyclopropane carboxylate) . 

The indoor environment is assumed to be a six-mat room 
(9.72 [m 2 ]) of a typical apartment in Japan. In conformity 
to the description of the above-mentioned aerosol can, 
spraying was effected at a rate of 60 [sec/ra 2 ] from a distance 
of about 20 [cm] from the floor and wall. The spraying site 
had a width of 10 [cm] or 15 [cm] from the boundary between 
the floor and the wall toward the floor, and a width of 10 
[cm] or 15 [cm] from the boundary toward the wall [see Fig. 
2]. The total spraying period was 2.5 minutes for 2.52 [m 2 ] 
in the region having a width of 10 [cm], and 3.7 minutes 
for 3.78 [m 2 ] in the region having a width of 15 [cm]. 

Five kinds of measurement were effected while the room 
conditions were changed as shown in Table 7. Namely, 
measurement in the case with continuous use (Example 1), 
case with changing ventilation rate (Example 2), case with 
an open window (Example 3), and case with air-conditioner 
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on (Example 4), and measurement for skin attachment ratio 
(Example 5) were performed. 
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In Example 1 , spraying was effected four times at a 
rate of once in two weeks . The air-conditioner was not 
actuated, and the windows of the room were closed. The room 
temperature and humidity were controlled so as to become 
5 25 [ °C] and 60 [%], respectively . The illumination condition 

was such that synthetic light illuminates indoors from 7 
a.m. to 10 a.m. and from 6 p.m. to 11 p.m., while natural 
light was provided through two windows . The ventilation 
rate was set to 1.5 [h _1 ] , which was minimally necessary for 

10 a Japanese woman and an infant to inhabit. The temperature, 

humidity, illuminance, and ventilation rate of the room 
actually measured under such a condition are shown in Table 
7, each being represented by the mean value in measurement. 
The temperature, humidity, and illuminance are measured at 

15 the center of the room. The ventilation rate was confirmed 

before and after the measurement. 

In Example 2, the illumination condition was in the 
dark as natural light was cut. For performing measurement, 
the ventilation rate was changed to 0.5 to 4 [h" 1 ] since the 

20 actual ventilation rate of Japanese apartment with closed 

windows is 0.5 to 3 [h -1 ]. 

In Example 3 , immediately after spraying, a window was 
opened for 5 minutes or 2 hours in conformity to the 
description of the aerosol can. While the illumination 

25 condition was basically in the dark (illuminance 0 [lux] 

in Table 7) , natural light entered when the window was open 
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(mean illuminance 1890 [lux] in Table 7). 

The sampling processing in each of the above-mentioned 
measurements was effected as follows. 

In Examples 1 to 4, the indoor air was sampled after 
5 a specified period from the spraying. Here, as for d- 

phenothrin and d-tetramethrin, at the positions 
corresponding to the respiration regions of an infant and 
an adult , air suction amount was adjusted such that particles 
having a diameter of 30 [/zm] or less can be collected, 
10 whereby not only the compounds in a gas state but also those 

in a suspended particle state were sampled. 

In Example 1, at the centers of floor, wall, and 
ceiling, strips of materials identical thereto each having 
a size of 5 [cm] * 5 [cm] were placed before spraying, whereby 
15 sampling was effected at the floor, wall, and ceiling. As 

for the wall, sampling was effected at three places 
consisting of the inflow and discharge sides of the air and 
a lateral side apart therefrom. 

In the collecting tube that had sampled the air, 
20 flooring strip, wallpaper strip, ceiling strip, hand, and 

denim cloth, d-phenothrin and d-tetramethrin were analyzed 
by the methods shown in Table 8 . In these methods , an organic 
solvent containing BHT { 2 , 6-di-t-butyl-p-cresol ) was used 
to extract these compounds, and cleanup was effected when 
25 necessary, so that analysis was effected by mass 

spectrometer (MSD) attached to a gas chromatography (GC). 
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In each process, in order to prevent the compound from 
vaporizing due to evaporation of the solvent, diethylene 
glycol or octanol is added. Both d-phenothrin (mass: 183 
[m/z]) and d-tetramethrin (mass: 164 [m/z]) were detected 
5 in 9 minutes under the GC condition noted in the following. 

As the apparatus, Shimadzu GC-MS QP-1100EX (Mode EI, 70 eV) 
was used with a DC-60 8 column (inner diameter: 0.53 [mm], 
length: 3 0 [m] , thickness: 0.8 [Am]) programmed to raise 
temperature at a rate of 10[°C/min] from 205 [°C], whereas 
10 the temperature of injection and ion source was 260 [°C], 
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The results of concentrations of d-phenothrin and 
d-tetramethrin in the air and the residual amounts in the 
floor, wall, and ceiling measured by the above-mentioned 
methods will be explained. 

In Example 1 , the concentration in the air was measured 
for 125 days including the time in which spraying was 
effected once (2.5 minutes ) per 2 weeks for 8 weeks . As shown 
in Fig. 7, the concentration in the air exhibited its maximum 
level immediately after spraying, where the average of 
d-phenothrin at the heights of 25 [cm] and 120 [cm] from 
the center of the floor was 752 [^g/m 3 ] and that of d- 
tetramethrin was 1040 [/Zg/m 3 ], and rapidly decreased 
thereafter. Here, since the half-life of each of both 
compounds was short, i.e., about 2 0 minutes, there were 4 
peaks in the 8 weeks of spraying period. The concentration 
in the air after the spraying period was less than 0.10 [ ju 
g/m 3 ] in both of them. Further, there was no substantial 
difference in concentration between the heights of 25 [cm] 
and 120 [cm] from the floor. Also, no increase was seen in 
the concentration in the air due to the increase in the number 
of spraying operations. In Fig. 7, concentrations which 
were less than the detection limit (0.05 [/Zg/m 3 ]) and thus 
could not be measured are plotted at the detection limit. 

The total spraying amount was 250 [g], including 1.11 
[g] of d-phenothrin and 1.3 8 [g] of d-tetramethrin, which 
equal 1.2 aerosol cans of 3 00 [ml] . Since 5 million aerosol 
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cans of 300 [ml] for residual spraying are annually sold 
in Japan, annual consumption per household is calculated 
as 0.13 can. Accordingly, it is seen that the spraying 
amount used for measurement (1.3 cans per room) is 
considerably greater than the actual annual consumption. 

As shown in Table 9, the mean concentration in the air 
during the spraying period and the mean annual concentration 
in the air including 8-week spraying period were calculated 
as 2.35 [yUg/m 3 ]/0.43 [Ag/m 3 ] for d-phenothrin, and 3.25 [ ju 
g/m 3 ]/0.53 [/Zg/m 3 ] f or d-tetramethrin . Consequently, it is 
seen that the amount of existence of d-phenothrin with 
respect to d-tetramethrin in the air corresponds to their 
ratio of existence within the aerosol can. In Table 9, for 
those less than the detection limit value, the mean 
concentration in the air was calculated by use of the 
detection limit value. 
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In Example 2 / concentrations of d-phenothrin and 
d-tetramethrin in the air after 3 . 7 minutes from the spraying 
with changed ventilation ratio were measured . Fig . 8A shows 
only the results of d-phenothrin. This has indicated that 
5 the half -lives in the air at ventilation ratios of 1 . 5 8 [ h" 1 ] , 

0.50 [h -1 ], and 4.14 [h _1 ] are 21 minutes, 58 minutes, and 
11 minutes, respectively. 

In Example 3 , immediately after the chemical was 
sprayed for 2.5 minutes, the window was opened, and the 

10 concentrations of d-phenothrin and d-tetramethrin in the 

air were measured. Fig. 8B shows the results of d- 
phenothrin. This has indicated that the half -lives in the 
cases where the window was opened for 5 minutes , the window 
was opened for 2 hours , and the window was kept closed are 

15 14 minutes, 4 minutes, and 20 minutes, respectively. Also, 

the mean concentration (2.88 [/zg/m 3 ] when opened for 5 
minutes, 1.19 [//g/m 3 ] when opened for 2 hours) was less than 
1/5 that of the closed state (14.8 [/Zg/m 3 ]). 

The measured values of Figs. 8A and 8B are mean values 

20 of the concentrations in the air respectively sampled at 

heights of 25 [cm] and 120 [cm] from the floor. Also, for 
those less than the detection limit and thereby incapable 
of measurement were plotted as the detection limit value 
(0.05 [/zg/m 3 ]). 

25 Table 10 lists the above-mentioned half-lives . These 

results have indicated that, in the residual spraying, the 
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ventilation of room, opening the window immediately after 
spraying in particular, is highly influential in the 
concentration in the air. For reference, Table 10 shows 
half-lives of compounds (d-resmethrin and d-tetramethrin) 
5 in the air when a chemical is sprayed in the space. 
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Meanwhile, in the state where the window is open as 
with Example 3, for example, the half-lives caused by 
photodegradation are estimated according to the above- 
mentioned Atmospheric Oxidation program from the structures 
of compounds to be 2 . 1 hours for d-phenothrin and 1 . 7 hours 
for d-tetramethrin. Since the half -lives in the air are 
shorter than the half -lives caused by photodegradation, the 
decreases in concentration of d-phenothrin and d- 
tetramethrin are considered to be mainly caused by 
ventilation rather than photodegradation. 

In Example 4, the mean concentration when the air- 
conditioner was used (11.6 [jug/m 3 ] for d-phenothrin, 15.5 
[/Zg/m 3 ] for d-tetramethrin) was substantially the same as 
the mean concentration when the air-conditioner was not used 
(14.8 [/Zg/m 3 ] for d-phenothrin, 20.2 [Ag/m 3 ] for d- 
tetramethrin) . This has indicated that the influence of the 
air convection and absorption of the air into the filter 
of the air conditioner upon the concentration in the air 
is substantially negligible. 

Figs. 9A and 9B and the above-mentioned Table 9 show 
the results of measurement of residual amounts at the floor, 
wall, and ceiling for 12 5 days in Example 1. The values of 
the residual amounts at the wall in Fig. 9B and Table 9 are 
means of values on the air inflow and discharge sides and 
their lateral side, whereas the residual amounts of wall 
and ceiling in Fig. 9B are plotted as the mean values of 



79 



SCC 97-06 



residual amounts at the wall and ceiling. 

As the air flows in and out, the residual amount at 
the wall on the inflow side (28 [/Zg/m 2 ] for d-phenothrin, 
35 [/Ug/m 2 ] for d-tetramethrin ) was less than the residual 
5 amount at the wall on the discharge side or lateral side 

(45 [/ig/m 2 ] and 58 [jLLg/m 2 ] on the discharge side, 33 [jU 
g/m 2 ] and 43 [/ig/m 2 ] on the lateral side). Also, while the 
residual amounts of d-phenothrin and d-tetramethrin 
enhanced as the number of spraying operations increased, 

10 they gradually decreased after spraying. Here, half-lives 

of d-phenothrin and d-tetramethrin at the floor were 
respectively 13 days and 20 days, whereas their half-lives 
at the wall /ceiling were respectively 31 to 41 days and 2 4 
to 75 days. These values are longer than those in the air, 

15 since the materials of floor, wall, and ceiling contain a 

large amount of antioxidants. Also, the residual amount at 
the floor, except for the site sprayed with the chemical, 
is greater than the residual amount at the wall or ceiling, 
since the suspended particles attach to the floor. 

20 In Example 5, a hand attachment test in which a hand 

is directly pressed against the floor where the compound 
remains, and a wiping test using a denim cloth instead of 
the hand are performed, so as to measure the skin attachment 
ratios of d-phenothrin and d-tetramethrin at the flooring. 

25 When Student's t-test (P S 0.05) was performed, the skin 

attachment ratios obtained by the wiping test (0.31 [%] and 



80 



SCC 97-06 



0.33 [%] respectively for d-phenothrin and d-tetramethrin 
after 6 hours from spraying, 0.17[%] and 0.09[%] after 24 
hours from spraying) were equal to or higher than the results 
obtained by the hand attachment test (0.22 [%] and 0.32[%] 
respectively for d-phenothrin and d-tetramethrin after 6 
hours from spraying, 0 . 06 [ % ] and 0 . 07 [ % ] after 24 hours from 
spraying). This has indicated that the safety of 
percutaneous exposure can be sufficiently estimated when 
the skin attachment ratio is determined by the wiping test 
alone. 

In the following, simulation results obtained 
according to the method of estimating an indoor behavior 
of a compound in accordance with the present invention will 
be explained. The condition here is the same as that of the 
above-mentioned Example 1. As the primary condition, the 
values described in Tables 1 to 4 of the first embodiment 
were inputted. 

When the behavior of d-phenothrin caused by residual 
spraying during the period after 2.5 minutes to 24 hours 
from the spraying was simulated, the results shown in Table 
11 were obtained. In the table, "suspension period" refers 
to a period in which the particles floating at the position 
lower than 2.4 [m] fall onto the floor or disappear in the 
air, whereas "transference amount" refers to the ratio of 
compound amount transferred within 24 hours with respect 
to the whole spraying amount. 
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According to these results, the main solvent in the 
large particles (diameter: 20 [jUm]) was completely- 
evaporated after 35 seconds, and disappeared from within 
the air after 2 hours. The main solvent in the small 
particles (diameter: 5 [jUm]) was completely evaporated 
after 2.2 seconds, and existed in the air for 30 hours. The 
chemical containing the particles in the air (40[%] of the 
amount of chemical immediately after spraying) was 
distributed to the spraying site by 1 . 0 [ % ] , floor by 7 . 0 [ % ] , 
wall by 0.44[%], and ceiling by 0.10[%], whereas the 
remaining 29[%] went outdoors, and 3.2 [%] was degraded in 
the air. Of the chemical at the spraying site (60[%] of the 
amount of chemical immediately after spraying), 0.16[%] of 
the sprayed amount was evaporated into the air within 24 
hours, 1.7[%] was degraded by light or oxidation, and the 
remaining 58[%] was left at the spraying site. 

When the simulation results concerning the 
concentration in the air and the residual amounts at the 
spraying site and floor /wall/ceiling and the actually 
measured values which is mentioned above were compared with 
each other, relationships shown in Figs. 10A and 10B were 
obtained. In these graphs, solid lines show simulation 
results , whereas plotted points indicate actually measured 
values. These actually measured values were derived from 
the results of Example 1 . 

From these graphs, it can be seen that the estimated 
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results according to the method of estimating an indoor 
behavior of a compound in accordance with the present 
invention are very close to the actually measured values . 
Namely, upon residual spraying, the most part of the chemical 
in the air is discharged outdoors or falls onto the floor. 
The chemical at the spraying site permeated through the floor 
or wall, and was slightly degraded or evaporated. 

In the following, the results obtained according to 
the method of evaluating safety of a compound in accordance 
with the present invention will be explained. Here, indoor 
exposure amounts of d-phenothrin and d-tetramethrin with 
respect to inhabitants including an infant when a chemical 
was residually sprayed at the boundary between the floor 
and wall for four times in eight weeks (sprayed for 2.5 
minutes for each time) were evaluated. 

The inhabitants inhale vapors or particles of d- 
phenothrin and d-tetramethrin existing in the air. Since 
humans inhale aerosol particles having a diameter not 
greater than 8 [ jUm] to the lungs, sampling the particles 
having a size not greater than 30 [/Zm] and vapor chemicals 
is sufficient for evaluating the exposure in inhalation. 

The following Table 12 shows the estimated exposure 
amount in inhalation calculated by use of the mean indoor 
concentration in the air (see "spraying period" in Table 
9), respiration amount, and exposure time as explained in 
item (iv) in the first embodiment. Here, the mean indoor 
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respiration amount of Japanese is 0.220 [L/min/kg] for an 
adult, and 0.235 [L/min/kg] for an infant. Also, the 
exposure time was presumed to be 24 [h/day] , taking account 
of the worst case where inhabitants are confined in the room 
for 24 hours after spraying. 
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In the table, "window closed" refers to the case where 
the window is closed throughout the exposure time, whereas 
"window open" refers to the case where spraying was effected 
in the state with the closed window and then the window was 
opened for 5 minutes (Example 3). 

Also, Table 12 shows the estimated percutaneous 
exposure amount calculated by use of the mean floor residual 
amount (see "spraying period" in Table 9), skin attachment 
ratio, contact area, and body weight as explained in item 
(v) in the first embodiment. Here, it was assumed that the 
contact area for an infant and an adult was 4.0 [mVday], 
the body weight of the infant was 10.2 [kg], and the body 
weight of the adult was 50 [kg]. As the skin attachment 
ratios for d-phenothrin and d-tetramethrin at the flooring, 
values evaluated by integration and extrapolation of the 
measured values in Example 5 were used. Namely, the skin 
attachment ratio was 0.31 [%] for d-phenothrin, and 0.25 [%] 
for d-tetramethrin. These values were used for calculating 
the estimated percutaneous exposure amount. 

Further, Table 12 shows the estimated oral exposure 
amount calculated by use of the estimated percutaneous 
exposure amount, hand surface area ratio, and oral 
transference ratio as explained in item (vi) in the first 
embodiment. Here, while the hand surface area ratio was 
assumed to be about 5[%], and the oral transference ratio 
was presumed to be the worst case of 100[%], only the 
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estimated oral exposure amount of the infant was calculated. 

Table 13 lists the inhalation safety coefficient, 
percutaneous safety coefficient, and oral safety 
coefficient obtained according to the estimated exposure 
amount in inhalation, estimated percutaneous exposure 
amount, and estimated oral exposure amount determined by 
the above-mentioned calculations. At the same time, the 
total safety coefficient determined from the above- 
mentioned three safety coefficients are shown in the table. 
Further shown in the table is the safety coefficient with 
respect to both of d-phenothrin and d-tetramethrin computed 
when the sum of the respective reciprocals of the safety 
coefficient (inhalation, percutaneous, or oral) of d- 
phenothrin and the safety coefficient (inhalation, 
percutaneous, or oral) of d-tetramethrin is determined and 
then is turned into the reciprocal. 



TOTAL 
(WINDOW OPEN) 


80,400 
116,000 


12,200 
15,800 


10,600 
13,900 


TOTAL 
(WINDOW CLOSED) 


21,300 
24,400 


2,780 
3,180 


2,460 
2,820 


ORAL 


923,000 


100,000 


90,300 


PERCUTANEOUS 


369,000 
1,810,000 


465,000 
2,280,000 


206,000 
1,010,000 


INHALATION 
(WINDOW OPEN) 


116,000 
124,000 


14,400 
15,900 


12,800 
14,100 


INHALATION 
(WINDOW CLOSED) 


23,100 
24,700 


2,880 
3,190 


2,560 
2,820 


SUBJECT 
COMPOUND 


d-phenothrin 
INFANT 
ADULT 


d-tetramethrin 
INFANT 
ADULT 


BOTH 

INFANT 
ADULT 
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Here, the non-influential amount of a rat concerning 
d-phenothrin was 18 . 4 [mg/kg/day] (= non- influential in-air 
concentration 210 [mg/m 3 ] * rat respiration amount 0.3 65 
[L/min/kg] * rat exposure time 4 [h/day]; 4-week test with 
exposure of 6 days a week) in inhalation, 1000 [mg/kg/day] 
(21-day test) for percutaneous exposure, and 125 
[mg/kg/day] ( 6-month test ) for oral exposure. On the other 
hand, the non-influential amount of a rat concerning d- 
tetramethrin was 3.22 [mg/kg/day] (= non-influential in-air 
concentration 4 9 [mg/m 3 ] x rat respiration amount 0.365 
[L/min/kg] * rat exposure time 3 [h/day]; 4-week test with 
exposure of 7 days a week) in inhalation, 1000 [mg/kg/day] 
(21-day test) for percutaneous exposure, and 10.75 
[mg/kg/day] (6-month test) for oral exposure. 

As a result, it has been confirmed that the safety 
coefficient of the infant is lower than that of the adult, 
and that the safety coefficients of d-phenothrin were higher 
than the safety coefficients of d-tetramethrin except for 
the percutaneous safety coefficient. Also, since the 
inhalation safety coefficient is the lowest of the three, 
it substantially determines the value of total safety 
coefficient when exposed by all of the three kinds of 
exposure routes . 

(2) Second Embodiment 

With reference to the drawings, the second embodiment 
of the present invention will be explained in the following. 
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The method of simulating an indoor behavior of a pesticidal 
compound in accordance with this embodiment estimates the 
behavior of the compound in the case where the floor in the 
first embodiment is constituted by a rug having ears of fiber. 
For convenience of explanation, members identical to those 
shown in the drawings of the previous embodiment will be 
referred to with numerals or letters identical thereto, 
without their explanations provided. Explained here is a 
case where only the kind of floor at the spraying site is 
a carpet (e.g., the floor spraying region indicated by c 
in Fig. 2 is assumed to be a carpet) , whereas the other part 
of the floor is treated as a mere plane. Accordingly, only 
the portion different from the first embodiment will be 
explained . 

Here, though the kind of floor (whole surface) may be 
a carpet, only the suspended particles fall onto the floor 
other than the spraying site, whereby the behavior of the 
floor does not greatly differ from that in the case of 
flooring. Accordingly, it is often only necessary for 
Fugacity capacity of the floor to be corrected with the 
particle-constituting solvent component content alone. 

First, the carpet is modeled as shown in Figs. 6A, 6B, 
and 6C. Namely, the carpet is divided into a plurality of 
carpet fibers (ears of fiber) planted on a substrate, and 
space portions existing between the fibers . When a chemical 
is residually sprayed in a room in which such a carpet covers 
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the whole floor surface, the chemical infiltrates into the 
carpet fiber at the spraying site, and accumulates in the 
space portion. 

While the method of simulating an indoor behavior of 
5 a pesticidal compound in accordance with this embodiment 

also conforms to the flowchart of Fig. 1, at SI, as the 
primary condition, the carpet-related primary condition 
shown in Table 14 is added to the primary condition of the 
first embodiment (see Tables 1 to 4 ) . Here, the cross- 
10 sectional area (A c ) of the space portion of the carpet refers 

to the area in which the space portion is in contact with 
the indoor air (see Figs. 6A, 6B, and 6C). 
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The examples of values in Table 14 are based on the 
values shown in Figs. 6A, 6B, and6C. Here, one space portion 
is set to a size of 0.1 [mm] x o.l [mm] * 3 [mm], and one 
square of carpet fiber is set to a size of 2.5 [mm] x 2.5 
[mm] x3 [mm]. Also, it is assumed that 16 [pieces ] of carpet 
fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

A c = (0.1 [mm] x o.l [mm] x 16 [pieces] )/l [cm] 
= 1.6 x 10" 3 [mVm 2 floor] 

A cl = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces] )/l 

[cm] 

= 0.192 [mVm 2 floor] 
A x = 1 - A c 

= 1 - 1.6 x 10" 3 [mVm 2 floor] 
A al = A x + A cl 

= 1.19 [m 2 /m 2 floor] 
Here, A x , which is the area occupied by the spraying site 
in the first embodiment, is replaced by the area in which 
the carpet fibers in 1 [m 2 ] are in contact with the indoor 
air. 

As the secondary condition at S2, the carpet-related 
secondary condition shown in Table 15 is added to the 
secondary condition of the first embodiment. Here, k L and 
V ir which are respectively transference speed of the 
compound in the spraying site and the volume of the spraying 
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site in the first embodiment (see Table 6) , are respectively 
replaced by the transference speed of the compound in the 
carpet fiber and the volume of the whole fiber of the carpet . 
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A method of computing the definite factor conditions 
and temporal change conditions shown in the above-mentioned 
Table 15 will be explained hereinafter as exemplified by 
the case using an aqueous aerosol, i.e., case where the 
particle-constituting solvent is water. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T) , and humidity (H), the evaporation 
constant (R d ) of the chemical entered into the space portion 
of the carpet is expressed as: 

C H = 0.24 + 0.46/f I ...(20) 

h =3.06xl0" 4 (7;-7;) I/3 j 

wherein T d is the temperature of the space portion of the 
carpet (corresponding to the wet-bulb temperature), and H m 
is the humidity (corresponding to the saturated humidity 
at the wet-bulb temperature) of the space portion at T d . 
Here, T d and H m can be determined from the above-mentioned 
"mass-based humidity table" according to the room 
temperature (T) and humidity (H). In the case of an 
oil-based product (aerosol), while the vapor constant of 
the chemical can be computed when properties of the water 
and oil are compared with each other and then are multiplied 
by R d , R d = 0 may be assumed in the case where the evaporation 
of the oil is slow. The chemical entered into the space 
portion of the carpet reduces its volume upon evaporation 
and finally disappears, leaving the compound. Thus left 
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compound permeates through the carpet fiber. 

From the diffusion coefficient (D c ) in the 
floor/wall/ceiling and the area (A al ) in which the carpet 
fiber is in contact with the compound, the volume (V\) of 
5 the whole fiber of the carpet is determined as follows : 

V l = ZyjDJtA^A, ...(21) 

The transference speed (k c ) of the compound in the 
product is expressed as k 4 /100 by use of k 4 (= GV 4 /A 4 ) which 
indicates the transference speed of the compound in the 

10 indoor air. 

From the above-mentioned transference speed (k c ) of 
the compound in the product, transference speed (k 4 ) of the 
compound in the indoor air, area (A t ) of the spraying site, 
cross-sectional area (A c ) of the space portion, Fugacity 

15 capacity (Z c ) of the chemical in the space portion determined 

by expression (24) which will be mentioned later, and 
Fugacity capacity (Z 4 ) of the air determined by expression 
(4) of the first embodiment, the transference coefficient 
(D c4 ) of the compound between the chemical in the space 

20 portion of the carpet and the air is expressed as follows: 

D , = ...(22) 

\l{k c AA t Z c ) + 1/(VW1,Z 4 ) 

From the above-mentioned transference speed (k c ) of 
the compound in the product, diffusion coefficient (D c ), 
area (A t ) of the spraying site, area (A cl ) in which the carpet 
2 5 fiber is in contact with the compound in the space portion, 
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Fugacity capacity (Z c ) of the chemical in the space portion, 
and Fugacity capacity (Z 1 ) of the floor determined by 
expression (2) of the first embodiment, the transference 
coefficient (D o1 ) of the compound between the chemical in 
the space portion and the carpet fiber is expressed as: 

D = 1 _ (2 3) 

l/(* e A el AZ e ) + l/(MciA z i) 

wherein the transference speed (k x ) of the compound within 
the carpet fiber is expressed as (D 0 /t) 0 ' 5 . 

The volume (V c ) of the chemical in the space portion 
is determined from the area (A t ) of the spraying site, amount 
of emission of the particle-constituting solvent (O a ), and 
evaporation constant (R d ) of the chemical in the space 
portion determined by expression (20). The volume (V c ) of 
the space portion decreases until the solvent in the chemical 
completely evaporates . 

From the vapor pressure (P 3 ), water solubility (C s ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and property change constant (a) of the product, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is expressed as: 
Z c = (e-"Kow + 1 - e al )C S /P S ...(24) 

Here, the property change constant (a) of the product 
can be computed according to the evaporation constant (R d ) 
of the chemical in the space portion determined by expression 
(20), content (C sd ) of the solvent (Xylene) contained in the 
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product, and dilution (X sol ) of the product. Namely, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is also expressed as: 

Z 0 = [ (ratio of organic solvent in product) K ow + (ratio 
of water in product)] * C s /P s 

By simultaneously forming this expression and the 
above-mentioned expression (24), the property change 
constant (a) of the product is computed. Here, in the case 
of the oil-based product (aerosol), Fugacity capacity (Z c ) 
of the chemical in the space portion is given by expression 
(5) of the first embodiment from the vapor pressure (P s ), 
water solubility (C s ), and octanol/water distribution 
coefficient (K ow ) of the compound - 

Consequently, the behavior of the compound in the 
space portion of the carpet is expressed in the form of 
differential equation concerning Fugacity (f c ) as follows: 



Also, expressions (15) and (17) in the first 
embodiment are respectively replaced by: 





Deposition(i) 
V - change 



-DM* -fx) 
- KV c ZJ' c 



Transference^ j) 
Transference(5) 
Degradation 



...(25) 
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—V 1 Z l = ^DJtA al A t ZJ x V -change 



dt 

F- change + Y n i A 1 AV i Z i f i Deposition{i) 

n ...(26) 

-D e i(fi~ fc) Transference^) 
~ D u(fi-fi) Transference^) 

- K 1 V 1 Z 1 fl Degradation 

Qt-V.Z, = -GV.Z.f. Ventilation 

dt 4 4 4 * U 

F- change n fii*{fA,~fi) Transference(i) 

-^D 4y (/ 4 -/ ; ) Transference^) -( 27 ) 

- K^V 4 Z A f Degradation 

- £> c4 (/ 4 - / c ) Transference(c) 
-DM-fi) Transference® 

In expression (26), terms of V-change, Deposition (i), 
Transference (c), Transference (1), and Degradation 
respectively indicate the change in volume (increase over 
time) of the carpet fiber, attachment accompanying the 
falling of particles, amount of transference of the compound 
between the carpet fiber and the space portion, amount of 
transference of the compound between the carpet fiber and 
the air, and change in amount of photodegradation of the 
compound - 

In expression (27), terms of Ventilation, 
Transference (i), Transference (j), Degradation, 
Transference (c), and Transference (1) respectively 
indicate the amount of discharge of the compound outdoors, 
amount of transference of the compound between the air and 
the particles, amount of transference of the compound 
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between the air and the floor/ceiling/wall, change in amount 
of photodegradation of the compound, amount of transference 
of the compound between the air and the space portion at 
the spraying site, and amount of transference of the compound 
5 between the air and the carpet fiber at the spraying site. 

Accordingly, Fugacity calculation at S3 is performed 
such that eight kinds of differential equations of 
expressions (16) and (18) in the first embodiment and the 
above-mentioned expressions of (25) to (27) are solved by 

10 Runge-Kutta-Gill method. The subsequent steps of S4 to S17 

are similar to those in the first embodiment. 

As explained in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
accordance with this embodiment, in the case where the kind 

15 of the spraying site is a carpet having ears of fiber and 

a space portion between the ears, a differential equation 
for Fugacity in the above-mentioned space portion of the 
carpet is added to seven kinds of simultaneous differential 
equations concerning Fugacity of the compound in the 

20 spraying site, two kinds of suspended particles, indoor air, 

floor, wall, and ceiling, so as to calculate each Fugacity. 

Consequently, even in the case where the chemical is 
accumulated in the space portion, so that the compound 
behaves differently from the case with a flat floor, the 

2 5 behavior of the compound can be estimated accurately, 

whereby various kinds of simulations can be performed. 
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( 3 ) Third Embodiment 

With reference to the drawings, the third embodiment 
of the present invention will be explained in the following. 
5 The method of simulating an indoor behavior of a 

pesticidal compound (hereinafter simply referred to as 
compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment mainly 
estimates the behavior of the compound in the case where 
- 10 the chemical is spatially sprayed in order to terminate 

indoor vermin such as mosquito or fly. Here, the chemical 
encompasses all kinds of insecticides including pyrethroid 
insecticidal compounds and organophosphorus insecticidal 
compounds . 

15 First, as shown in Fig. 11, an environment is 

classified into seven kinds of media consisting of large 
particles ( aerosol particles ) 2, medium particles (aerosol 
particles) 8, and small particles (aerosol particles) 3 
having diameters different from each other, whole indoor 

20 air 4, floor 5, wall 6, and ceiling 7. 

The particles are divided into three kinds since their 
behavior mainly depends on their diameter. Though the 
particles are distributed in the same site immediately after 
spraying, as time passes, the large particles 2 are 

25 distributed into a spraying zone SZ 2 near the floor due to 

their weight, the medium particles 8 into a spraying zone 
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SZ 8 in the middle of the room, and the small particles 3 into 
a spraying zone SZ 3 near the ceiling 7 since they are lighter 
than the other two kinds of particles . 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned seven kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is sprayed in the indoor space, i.e., 
temporal change of compound distribution in the room, is 
estimated and analyzed. 

Here, since the floor 5 mainly comprises three kinds, 
i.e., tatami, flooring, and carpet (rug), the differential 
equation is formed so as to correspond to the kind of floor. 
Explained in this embodiment is the case where the floor 
5 is made of tatami. 

With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 
be mainly divided into a step of dividing an indoor 
environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2 ) ; a step of determining the fugacity 
of the compound in each medium from the differential equation 

(53) ; a step of determining the indoor behavior of the 
compound from the fugacity of the compound in each medium 

(54) ; a step of changing, in response to a fluctuation in 
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mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 
see Fig. 5) ; and a step of evaluating safety of the compound 
with respect to a human body according to the indoor behavior 
5 of the compound (S5 to S17). 

First, at step (hereinafter abridged as S) 1 , a primary 
condition is inputted. The primary condition comprises 
physicochemical properties of the compound (see Table 16) , 
indoor environmental behavior properties of the compound 
10 (see Table 17), indoor environment (see Table 18), product 

properties (see Table 19) , and a spray condition (see Table 
20) . 
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Here, the examples of values in Tables 16 to 20 assume 
those of a six-mat room (9.72 [m 2 ]) of a typical apartment 
in Japan in summer (see Fig. 14). 

Subsequently, using the primary condition, a 
5 secondary condition is determined by calculation (S2 ) . The 

secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 21) determined by the 
primary condition alone, and a temporal change condition 
accompanying a temporal change (see Table 22). Their 
10 calculations will be explained later in detail. 
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TABLE 21 



DEFINITE FACTOR CONDITION 


SYMBOL 


MT TlVTRPl? CW AUCnCOT D A DTIPI EC TA/TA/TCnT A TT3T V ACTED CUD A VTMn 

EVAPORATION CONSTANT OF AEROSOL PARTICLES 
VOLUME RATIO OF COMPOUND IN PRODUCT 
DIFFUSION COEFFICIENT OF COMPOUND IN AIR 
SURFACE AREA OF FLOOR / WALL / CEILING 
DIFFUSION COEFFICIENT OF FLOOR / WALL / CEILING 


nn; (i=2 8 31 
a 
Ra 
D ca 

Aj 0=5,6,7) 
D c 


FUGACITY CAPACITY 

OF AEROSOL PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Zi (i=2,8,3) 
Z 4 

Zj (j=5,6,7) 



TABLE 22 



TEMPORAL CHANGE CONDITION 


SYMBOL 


AEROSOL PARTICLE 

DIAMETER 
FALLING SPEED 
FLOATING NUMBER 


di 0=2,8,3) 
vi 0=2,8,3) 
ni 0=2,8,3) 


AEROSOL ZONE WIDTH 


L zi (i=2,8,3) 


TRANSFERENCE SPEED OF COMPOUND 

IN AEROSOL PARTICLE 
AIR 

FLOOR / WALL / CEILING 


ki 0=2,8,3) 
k 4 

kj (j=5,6J) 


TRANSFERENCE COEFFICIENT OF COMPOUND 
BETWEEN AEROSOL PARTICLE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN AIR AND FLOOR / WALL /CEILING 


D i4 0=2,8,3) 
D 4 j 0=5,6,7) 


VOLUME OF AEROSOL PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


Vi 0=2,8,3) 
Vj 0=5,6,7) 
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Using the secondary condition, seven kinds of Fugacity 
are calculated (S3). Namely, seven kinds of differential 
equations concerning the large, medium, and small particles 
2, 8, and 3, air 4, floor 5, wall 6, and ceiling 7 are 
5 simultaneously formed and are solved by Runge-Kutta-Gill 

method, whereby seven kinds of Fugacity are computed over 
time. Here, an estimation nick time width (minute time 
unit) set when solving the differential equations is 
automatically set so as to be varied in response to a 

10 fluctuation in mass balance of the compound. 

At S4, using thus computed Fugacity (f Lr i = 2, 8, 3) 
of the small, medium, and large particles 2, 8, and 3, and 
Fugacity (f 4 ) of the air 4, temporal concentration of the 
compound in the air is computed; whereas temporal floor 

15 residual amount of the compound on the floor 5 is computed 

by use of Fugacity (f 5 ) of the floor 5. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 

2 0 an estimated exposure amount in inhalation indicating a 

degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 

25 estimated exposure amount in inhalation (S7). At S8, the 

inhalation safety coefficient is compared with a reference 
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value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the floor 5 to which the chemical is attached 
(S9). Thereafter, at S10, it is judged whether to perform 
a percutaneous safety evaluation or not. If the 
percutaneous safety evaluation is to be performed, a 
percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at SI 2, the 
operation returns to SI, where alteration of the primary 
condition is considered. 

If the percutaneous safety evaluation is not selected 
at S10, calculated at S13 according to the estimated 
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percutaneous exposure amount is an estimated oral exposure 
amount indicating the degree of exposure in the case where 
the chemical attached to a hand or the like is taken orally 

(513) . Subsequently, according to the estimated oral 

5 exposure amount, an oral safety coefficient is computed 

(514) . This oral exposure may occur, in particular, when 
an infant puts a chemical-attached hand into the mouth. At 
S15, as with the safety evaluation in inhalation, the oral 
safety coefficient is compared with a reference value 

10 defined in each country, whereby the safety is evaluated. 

If there is a problem in safety at S15, the operation returns 
to SI, where alteration of the primary condition is 
considered. 

Finally, if it is judged to be safe at S8, S12, or S15, 
15 overall safety is evaluated. Here, the sum of respective 

reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
20 (S16). This overall safety coefficient is compared with a 

reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4, S6, S7, S9, Sll, S13, and S14 will be explained in detail. 
25 (i) Secondary Condition Calculation (S2) 

Calculations of the definite factor condition shown 
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in Table 21 and the temporal change condition shown in Table 
22 will be explained. 

The particle number (n oi , i = 2, 8, 3) of the large, 
medium, and small particles 2, 8, and 3 immediately after 
5 spraying is determined from the spraying rate, particle 

diameter (d oi , i = 2, 8, 3) immediately after spraying, 
distributions of particles (three kinds), and spray time. 
The particle number (n oi ) is constant until the bottom of 
the spraying zone (SZ ir i = 2, 8, 3) reaches the floor 5. 

10 According to surface temperature (T d ) of particles, 

properties of particle-constituting solvent (vapor 
pressure P d , molecular weight M d , and specific gravity p d ), 
and indoor environment (air diffusion coefficient D air , vapor 
pressure P°o, and room temperature T°°), the evaporation 

15 constant ( a ) of the large, medium, and small particles 2, 

8, and 3 is defined as: 

wherein R is a gas constant. 

Here, the upper parts of the above-mentioned Tables 

20 18 and 19 exemplify the case where the spatially sprayed 

aerosol is of oil base. In this case, T°° (temperature at 
the site far from aerosol particles) and T d are set to room 
temperature (T) , whereas poo (vapor pressure of oil at the 
site far from aerosol particles) is set to zero. 

25 In the case where the spatially sprayed aerosol is of 
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water base, it is necessary for P d , M df p d , T d , and P°° i n 
expression (28) to be changed to values based on properties 
of water. Specifically, the above-mentioned factor values 
when the room temperature (T) is 298 [K] (= 25 °C) and the 
5 relative humidity is 60[%RH] are exemplified in the lower 

parts of Tables 18 and 19. 

Here, T d can be determined, according to the room 
temperature (T) and relative humidity (H), from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 

10 (Theory and Calculation of Chemical Machines)" (Second 

Edition) (Saburo Kamei ed. , Sangyo Tosho) or the like. On 
the other hand, P d and poo can be computed from the following 
expressions : 

logioPd = 10.23 - 1750/(T d -38) 

15 log 10 Poo = 8.23 - 1750/(T-38) + log lo 0 

wherein 0 is an indoor relative humidity. 

From the vapor pressure (P 3 ) and melting point (T H ) 
of the compound and surface temperature (T d ) of particles, 
Fugacity capacity (Z if i = 2, 8, 3) of the large, medium, 

20 and small particles 2, 8, and 3 is expressed as: 

...(29) 

P/=F 5 exp {6.79(^/^-1)} 

Though the vapor pressure of the compound in a liquid state 
(P L S ) is computed by use of T M and T d here; in the case where 
T M cannot be obtained, P L S may be set identical to the vapor 
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pressure of the compound in a solid state (P s ) . 

From the room temperature (T), Fugacity capacity (Z 4 ) 
of the air 4 is expressed as follows: 

Z 4 = — ...(30) 
RT 

From the vapor pressure (P s ), water solubility (C s ), 
and octanol /water distribution coefficient (K ow ) of the 
compound, and from the oil component content ( p^, j = 5, 
6, 7) of the material constituting the floor 5, wall 6, and 
ceiling 7, Fugacity capacity (Z jA j =5, 6, 7) of the floor 
5, wall 6, and ceiling 7 is expressed as: 
Zj = Pj K ow C s /P s ...(31) 

Since the volume of particles decreases over time due 
to evaporation of the particle-constituting solvent 
component, during the evaporation of particle-constituting 
solvent component, using the diameter (d oi ) of particles 
immediately after spraying (t = 0) and the evaporation 
constant ( a ) determined from expression (28) , the diameter 
(di, i = 2, 8, 3) of the large, medium, and small particles 
2,8, and 3 at time t is represented as shown in the following 
expression ( 3 2 ) : 
d, = 4^ - 2at ...(32) 
4 - ..-(32') 

By contrast, the diameter (dj of particles after 
complete evaporation of the particle-constituting solvent 
component is represented as shown in the above expression 
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(32' ) by use of the ratio of volume of the compound in the 
product (R a ) and the diameter (d oi ) of particles immediately 
after spraying. 

From the compound content (C a ) and specific gravity 
of the compound ( p ) , the ratio of volume of the compound 
in the product (R a ) is expressed as Cjp. 

Since movement of the particles is under the control 
of gravity and air resistance, according to Stokes ' law using 
the specific gravity ( p d ) of the particle-constituting 
solvent and the diameter (d^ of particles determined by 
expression (32), the falling speed (v if i = 2, 8, 3) of the 
large, medium, and small particles 2, 8, and 3 before 
complete evaporation of the particle-constituting solvent 
component is represented by the following expression (33): 

1877 

S c =1 + {6.32 + 2.01exp(-8.322xl0 6 4)} ...(33') 

7.6x10 d i 

v-.-tjplpM ...(33") 

wherein g is gravitational acceleration, 77 is the 
coefficient of viscosity of the air 4, S c is a sliding 
correction coefficient , and j3 is a speed coefficient- Here, 
while particles do not conform to the Stokes' law when they 
become small, the sliding correction coefficient (S c ) , which 
is a coefficient for correcting this phenomenon, is 
represented by the above expression (33' ) according to the 
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diameter (d £ ) of particles determined by expression (32). 

On the other hand, the falling speed (v £ ) of the 
particles after the complete evaporation of particle- 
constituting solvent component is represented by the above 
5 expression (33") using the specific gravity ( p ) of the 

compound, specific gravity ( p d ) of the particle- 
constituting solvent, speed coefficient ( j8 ) , and diameter 
(di) of particles determined by expression (32'). 

From the spray height, aerosol zone width (L 2i ) 

10 determined by expression (35) which will be mentioned later, 

falling speed (vj of particles determined by expression 
(33) or expression (33"), and the above-mentioned particle 
number (n oi ) immediately after spraying, the floating number 
(n if i = 2, 8, 3) of the large, medium, and small particles 

15 2, 8, and 3 after the bottom of the spraying zone (SZi) reaches 

the floor 5 is expressed as: 

7^ -( 34 ) 

wherein t xi is the time necessary for the bottom of the 
spraying zone (SZ L ) to reach the floor 5, and t < t xi . Here, 

2 0 G is the ventilation rate, taking account of the fact that 

the large, medium, and small particles 2, 8, and 3 not only 
disappear as attaching to the floor, wall, and ceiling but 
also are discharged outdoors by ventilation. Also, L zi is 
the dispersion width of suspended particles upon spraying, 

25 which may be either determined by expression (35) or 
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simplified by use of the dispersion width immediately after 
dispersion as it is. 

The width of aerosol zone (L zi/ i = 2, 8, 3) is the 
width (height) of the spraying zone after the bottom of the 
spraying zone (SZ t ) reaches the floor 5. Since the spraying 
zone (S zi ) is absorbed by the floor 5 while changing the 
falling speed (Vi) of particles determined by expression 
(33) or expression (33"), from the spray height and the 
above-mentioned falling speed (Vi), the aerosol zone width 
is defined as: 



wherein t yi is the time necessary for the spraying zone (SZi) 
to be completely absorbed by the floor 5 - 

From the transference speed (Jq, i = 2, 8, 3) of the 
compound in the large, medium, and small particles 2, 8, 
and 3, transference speed (k 4 ) of the compound in the air 
4, diameter (dj determined by expression (32) or (32' ) , and 
Fugacity capacity (Zj_) and Fugacity capacity (Z 4 ) determined 
by expressions (29) and (30), the transference coefficient 
(D i4/ i = 2, 8, 3) between the large, medium, and small 
particles 2,8, and 3 and the air 4 is represented as follows: 

n = 1 _ (36 ) 

l/(£.A ; Z ; ) + l/{(£ 4 + v.)A ; Z 4 } 

wherein Ai (= 7td L 2 ) is the surface area of particles. From 
the ventilation rate (G) , the transference speed (k 4 ) of the 
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compound in the air 4 is expressed as GV 4 /A 4 . From the 
transference speed (k 4 ) of the compound in the air, the 
transference speed (k ± ) of the compound in the particles is 
expressed as k 4 /100. Here, as shown in Fig. 14, the 
5 above-mentioned A 4 is the cross-sectional area of the room 

with respect to the direction of movement of the air 4 
(arrowed direction in the drawing). 

Using the diffusion coefficient (D Q ) of compound in 
the floor 5, wall 6, and ceiling 7, the transference speed 

10 (kj) of the compound in the floor 5, wall 6, and ceiling 7 

is expressed as (D c /t) 0,5 . Here, D c can be computed as follows . 
Namely, D c can be determined when the diffusion coefficient 
(D ca ) in the air is multiplied by 10" 9 , whereas D ca can be 
computed from properties of the compound (e.g., structural 

15 formula, molecular weight, and the like) according to Wike 

and Lee method ( "Handbook of Chemical Property Estimation 
Methods," McGraw-Hill Book Company, 1982), for example. 

From the half-life of transference to air (T 4j , j = 
5, 6, 7) in each of the floor 5, wall 6, and ceiling 7, room 

20 size (V 4 ), volume (Vj) of the floor 5, wall 6, and ceiling 

7 determined by expression (39) mentioned later, and 
Fugacity capacity (Z 4 ) and (Zj) determined by expressions 
(30) and (31), the transference coefficient (D 4j , j= 5, 6, 
7 ) of the compound between the air 4 and the floor 5 , wall 

25 6, and ceiling 7 is represented as the following expression 

(37): 
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D 4J - , , , « -(37) 

° 4i = l/^AjZj+l/ikjAjZ .) "" (37) 

The above-mentioned transference coefficient (D 4j ) 
may also be determined by the expression represented by the 
above-mentioned expression (37') from the speed of the 
5 compound in the air (k 4 ), speed (kj) of the compound in the 

floor 5, wall 6, and ceiling 7 , surface area (A jx j = 5, 6, 
7) of the floor 5, wall 6, and ceiling 7, and Fugacity 
capacity (Z 4 ) and (Z 3 ) determined by expressions (30) and 
(31). The surface area (A 3 ) of each of the floor 5, wall 

10 6, and ceiling 7 is determined by the size of the room (V 4 

= L 4 x w 4 x height ) . 

Here, the half-life of transference to air (T 4j ) in 
expression (37 ) is hard to actually measure and varies under 
various conditions, thus making it difficult to determine 

15 an accurate value thereof. By contrast, all the parameters 

in expression (37 ' ) can be easily determined by calculation 
without necessitating measurement. Accordingly, the 
transference coefficient (D 4j ) is more preferably determined 
by expression (37') than by expression (37). 

20 Using the diameter (d L ) of particles determined by 

expression (32) or expression (32' ), the volume (V if i = 2, 
8, 3) of the large, medium, and small particles 2, 8, and 
3 is expressed by: 
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V,-^df -(38) 
o 

Assuming that each of the floor 5, wall 6, and ceiling 
7 before spraying the chemical is like a thin film and that, 
as the chemical infiltrates into the film, its thickness 
increases so as to enhance the volume thereof; from the 
surface area (Aj) of the floor 5, wall 6, and ceiling 7 and 
diffusion coefficient (D c ) of the compound in the floor 5, 
wall 6, and ceiling 7, the volume (V jA j =5, 6, 7) of floor 
5, wall 6, and ceiling 7 becomes: 
Vj-2y[D~tAj •••(39) 

(ii) Fugacity Calculation (S3) 

The behavior of the compound in the large, medium, and 
small particles 2, 8, and 3 is expressed in the form of 
differential equation concerning Fugacity (f f ) as: 

^-V.Z. = -od.ZJ. V- change 

F - change - D i4 (f ; - f 4 ) Transference (4) 

- KVZ, f. Degradation ■ ■ ■ (40) 

Here, the terms of V-change, Transference (4), and 

Degradation respectively indicate volume change (decrease 

over time) of the large, medium, and small particles 2, 8, 

and 3, amount of transference of the compound between the 

particles and the air 4, and change in amount of 

photodegradation of the compound. Here, among the 

photodegradation constants shown in the above-mentioned 

Table 17, it is difficult to actually measure the 
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photodegradation constant (K 4 ) of the compound in the air, 
whereby it is not usually easy to determine this value. 
Therefore, when the actually measured value is unavailable, 
it may be determined as being computed by Atmospheric 
5 Oxidation program (Atkinson et al., 1984, Chem. Rev. Vol. 

84, pp. 437-470) using the structural formula of the compound . 
Also, the degradation constant in the particles and in the 
interior materials (floor 5, wall 6, and ceiling 7) may be 
set to zero in view of the safety of inhabitants when the 
10 actually measured values are unavailable. Here, the term 

of V-change is effective only while the particle- 
constituting solvent component is evaporating, whereas 
Fugacity (fj is effective only while the spraying zone (SZJ 
exists . 

15 The behavior of the compound in the air 4 is expressed 

in the form of differential equation concerning Fugacity 
(f 4 ) as: 

—V.Z. = -GV A Z.f A Ventilation 
dt 4 4 4 474 

F - change - ^"iA^CA ~ /; ) Transference^) 

- D 4j (/ 4 - fj ) Transference^]) 
-K 4 V 4 Z 4 f 4 Degradation •••(41) 

Here, the terms of Ventilation, Transference (i), 
20 Transference (j), and Degradation respectively indicate 

amount of discharge of the compound outdoors, amount of 
transference of the compound between the air 4 and the 
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particles, amount of transference between the air 4 and the 
floor 5, wall 6, and ceiling 7, and change in amount of 
photodegradation of the compound. 

The behavior of the compound in the floor 5, wall 6, 
5 and ceiling 7 is expressed in the form of differential 

equation concerning Fugacity (f 5 ), (f 6 ), (f 7 ) as: 



•••(42) 



•(43) 



•(44) 

10 Here, the terms of V-change, Transference (4), and 

Degradation in each expression respectively indicate volume 
change (increase over time) of the floor 5, wall 6, and 
ceiling 7, amount of transference of the compound between 
the floor 5, wall 6, ceiling 7 and the air 4, and change 



V 5 Z 5 



dt 

F - change 



-jDJtA 5 ZJ 5 
-K 5 V 5 Z 5 f 5 



V - change 

Deposition(i) 
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Degradation 



4ft 
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in amount of photodegradation of the compound. Also, 
Deposition (i) indicates attachment accompanying the 
falling of the particles. Here, y 5 , y sr and y 7 indicate 
ratios of the large, medium, and small particles 2, 8, 3 
attaching to the floor 5 , wall 6 , and ceiling 7 , respectively, 
and mainly depend on the ventilation rate (e.g., when the 
ventilation rate is 1.58 [1/hour], 7 5 = 0.93, y 6 = 0.06, 
7 7 = 0.01). As 7 5 , 7 6 , and 7 7 are introduced, the behavior 
of the compound can be estimated accurately. When 7 5 , y 
6 , and 7 7 are unavailable, it is preferable to presume the 
worst case in terms of safety evaluation and perform 
calculation with y 5 = 1, y s = 0, and y 7 = 0. 

The above-mentioned seven kinds of differential 
equations (40) to (44) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
(f a to f B ) . 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
width (dt) which is a minute time unit. Namely, the 
estimation nick time width is used such that solutions of 
the simultaneous differential equations are initially 
determined at a time (t 0 ), and then solutions of the 
simultaneous differential equations are determined at a 
time (t 0 + dt) to which the estimation nick time width is 
added . As solutions are obtained while estimation nick time 
widths are successively added, temporally changing Fugacity 
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can be determined. Theoretically, as the set time of the 
estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 
calculation time . By contrast , when the set time is too long , 
5 solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 
large change. 

10 Specifically, mass balance is always confirmed such 

that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 

15 balance starts fluctuating. For example, when the 

fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 

compound input amount/ ( existing amount + degrading 

2 0 amount + discharging amount) = 0.95 to 1.05 

namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is 
determined by the compound content (C a ), spray rate, and 

25 spray time. Since the temporal amounts of compound in the 

seven kinds of media are determined by the simultaneous 
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differential equations, they are summed up so as to compute 
the existing amount as follows : 



existing amount = ^ n i fV i Z i + ^ f k V k Z k 

7 

degrading amount = ^ n i K i f i V i Z i + ^K k f k V k Z k 
discharging amount = GV 4 Z 4 f 4 



■•(45) 



The degrading amount and discharging amount are as 
represented above. 

With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
explained. 

First, an initial value of estimation nick time width 
(dt ) is inputted (S21 ) . Then, an upper limit set value (e.g. , 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 
is the lower limit of difference in mass balance are inputted 
(S22) . Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27) . If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 
(or t 0 ) and the mass balance at t = t + dt (or t, + dt) is 
at least the upper limit set value or not (S28). If it is 
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judged to be at least the upper limit set value at S28, then 
solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S2 8, it is judged 
at S30 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S30, since solutions are not 
influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width (dt) or not (S32) . When the estimation nick time width 
(dt) is not greater than the maximum value at S32, since 
solutions do not diverge, the estimation nick time width 
(dt ) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S3 2, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S30, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
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are continued without changing the estimation nick time 
width (dt). 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
5 process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S27 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
10 case where the fluctuation is within the range of ±5%. 

Accordingly, the calculation is stopped once (S34), the 
lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
15 set without being held constant, while monitoring the mass 

balance fluctuation, Fugacity can be computed accurately 
and efficiently. 

(iii) Computation of Temporal Concentration in Indoor Air 
and Floor Residual Amount (S4) 

20 The temporal concentration of the compound in the 

indoor air is computed when Fugacity (f L ) of the large, medium, 
and small particles 2, 8, and 3, determined by the 
above-mentioned item (ii), multiplied by their Fugacity 
capacity (Z ± ) , and Fugacity (f 4 ) of the air 4 multiplied by 

25 Fugacity capacity (Z 4 ) of the air 4 are summed up. Here, 

the large particles 2 or medium particles 8 may not be inhaled 
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by a human body depending on the kind of chemical. In such 
a case, calculation is effected by use of Fugacity (f 3 ) of 
the small particles 3 alone, for example. 

The floor residual amount of the compound is computed 
5 when Fugacity (f 5 ) of the floor 5 is multiplied by Fugacity 

capacity (Z 5 ) of the floor 5. 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 

10 indoor air forms a curve shown in Fig. 3A, for example. This 

concentration curve is integrated, an accumulated 
concentration of the compound during a specific period (~t 1 
to t 2 ) is determined ( see Fig . 3B ) , and the mean concentration 
in the indoor air is computed from thus determined value. 

15 While an arbitrary period is set as the specific period 

depending on the object, an appropriate period is usually 
set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 

2 0 concentration in the indoor air, amount of respiration, and 

exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] * amount of 

25 respiration [mVkg/min] * exposure time [min/day] 

is effected. Here, as the above-mentioned amount of 



131 



SCC 97-06 



respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 

The inhalation safety coefficient is computed from a 
non- influential amount concerning inhalation toxicity 
examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 

inhalation safety coefficient = 

inhalation non-influential amount [mg/kg/day]/ 
estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned floor residual amount forms a 
curve such as that shown in Fig. 4A, for example. This 
residual amount curve is integrated, the accumulated 
residual amount of the compound during a specific period 
(tj to t 2 ) is determined (see Fig. 4B) , and the mean residual 
amount is computed from thus determined value. While an 
arbitrary period is set as the specific period depending 
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on the object, an appropriate period is usually set in view 
of the method of use of the product and the test period of 
toxicity data. 

Then, according to the mean floor residual amount, 
5 skin attachment ratio, contact area, and body weight, the 

estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean floor residual amount [mg/m 2 ] * skin attachment ratio 

10 [%] x contact area [ m 2 /day ]) /body weight [kg] 

is effected. Here, as the contact area, a published value 
(e.g., 4 [mVday] ) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor 5 where the compound 

15 exists. As this value, a published value or a value 

experimentally obtained from a model may be used. 

A model experiment method for the skin attachment 
ratio is as follows. A weight ( 8 cm x 8 cm x 8 cm; 4.2 kg) 
is placed on a denim cloth (8 cm x 10 cm) with a pressure 

20 similar to that of an infant in contact with a floor, and 

the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 

25 therebetween, the skin attachment ratio is obtained. It has 

been confirmed that the skin attachment ratio obtained by 
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this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants . 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 

percutaneous non-influential amount [mg/kg/day ] / 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
are not many published values . Accordingly, a more accurate 
value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day ]/( estimated percutaneous 
exposure amount [mg/kg/day] * percutaneous absorption ratio 
[%] ) 

Here, when the percutaneous absorption ratio is unknown, 
employed is a national guideline (e.g., 10%) which usually 
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exists . 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
5 obtained in the above-mentioned item (v) , hand surface area 

ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 

10 estimated percutaneous exposure amount [mg/kg/day] * hand 

surface area ratio [%] x oral transference ratio [%] 
is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 

15 transference ratio is a hypothetical value, which is set 

to 100%, for example. 

In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 

2 0 amount from tableware or food to mouth be added to the 

estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 

25 the amount of the harmful material remaining in tableware, 

tableware use area which is the sum of tableware surface 
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areas, and oral transference ratio from tableware, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] x tableware use area 
5 [mVday] x oral transference ratio [%]/body weight [kg] 

is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio ) . As the tableware contamination ratio , 
an actually measured value (e.g., 9%) or a hypothetical value 

10 may be used. 

The oral safety coefficient is computed from the 
non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 

15 as: 

oral safety coefficient = 

oral non-influential amount [mg/kg/day]/ estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 

20 estimating an indoor behavior of a pesticidal compound in 

this embodiment, when a chemical containing the above- 
mentioned compound is spatially sprayed, differential 
equations concerning Fugacity of the large, medium, and 
small particles 2, 8, and 3, air 4, floor 5, wall 6, and 

25 ceiling 7 are simultaneously formed and are solved, and the 

indoor behavior of the compound is estimated according to 
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thus obtained solution. Here, the estimation nick time 
width is variably set, while constantly confirming mass 
balance of the compound indoors after the spraying, so that 
the amount of input of the chemical indoors and the resulting 
5 solution coincide with each other. 

Accordingly, since mass balance of the compound after 
the spraying is always monitored such that the amount of 
input of the compound indoors and the resulting solution 
coincide with each other, thereby varying the estimation 

10 nick time width; the estimation nick time width is set longer 

when the mass balance fluctuates a little, whereas it is 
set shorter when the mass balance starts fluctuating greatly . 
Namely, when solving simultaneous differential equations 
including a parameter accompanying temporal change, the 

15 estimation nick time width is automatically set in response 

to the fluctuation in mass balance. Consequently, when 
processed by a computer, an accurate solution can be obtained 
in a short time. 

The method of evaluating safety of a pesticidal 

2 0 compound in accordance with this embodiment uses the 

estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is spatially sprayed. 

Accordingly, the safety of the compound with respect 

2 5 to the human body can be evaluated accurately in a short 

time. As a consequence, when formulating a chemical such 
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as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 

Though the kind of the floor 5 is assumed to be tatami 
in this embodiment, differential equations similar to those 
mentioned above may also be formed in the case of flooring. 
In this case, however, when determining fugacity capacity 
(Z s ) of the floor 5 in equation (31), it is necessary for 
the particle-constituting solvent component content (see 
Table 18) to be changed from that of tatami to that of 
flooring. 

Though Fugacity is determined by use of Runge- 
Kutta-Gill method in this embodiment, other methods may be 
used for solving differential equations . Runge-Kutta-Gill 
method, however, is preferably used since a program for the 
above-mentioned differential equations can be easily made 
by Basic. Also in the case where differential equations are 
solved by a method other than Runge-Kutta-Gill method, 
similar effects can be obtained when the estimation nick 
time width is set as mentioned above. 

In the following, examples (Examples 6 to 10) of the 
method of simulating an indoor behavior of a pesticidal 
compound in accordance with the third embodiment will be 
explained. 

The evaluation of the method of estimating an indoor 
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behavior of a pesticidal compound in these examples is 
effected, under various experimental conditions (Example 
6 to Example 10), as the results obtained by estimation and 
the actually measured values obtained by measurement of the 
actual indoor behavior of the compound are compared with 
each other. The respective experimental conditions of 
Examples are as follows. 

Example 6 uses, as the spatially sprayed chemical, an 
oil-based insecticide containing pyrethroid. The indoor 
environment is a typical apartment room (3.6 [m] x 3.6 [m] 
x 2.4 [m] ) in Japan . The floor is a flooring in which a wooden 
floor is covered with a polyurethane resin, whereas the wall 
and ceiling are covered with a polychlorovinyl wallpaper. 
The windows are completely closed during and after the 
spatial spraying. The spatial spraying is uniformly 
effected for 10 seconds at a height of 1 . 6 [m] from the floor. 
The other conditions are shown in Tables 16 to 20. 

The conditions of Example 7 are the same as those of 
Example 6 except that its ventilation rate is 0.5 [1/hour] . 

The conditions of Example 8 are the same as those of 
Example 6 except that its ventilation rate is 4 . 1 [1/hour] . 

The conditions of Example 9 are the same as those of 
Example 6 except that the windows are opened for only 5 
minutes after 5 minutes from the spraying . 

The conditions of Example 10 are the same as those of 
Example 6 except that the windows are opened for only 2 hours 



139 



SCC 97-06 



after 5 minutes from the spraying. 

Figs. 12A and 12B show the simulation results and 
actually measured values of the compound concentration in 
the air under the respective conditions of the above- 
5 mentioned Examples. Also, Fig. 12C shows the simulation 

results of the residual amounts in floor and ceiling in 
Example 6. In these charts, solid lines, dotted lines, and 
dashed lines indicate simulation results, whereas plotted 
points indicate actually measured values . 
10 From the above-mentioned results, it can be seen that 

the estimated results according to the method of simulating 
an indoor behavior of the compound in accordance with the 
present invention are very close to the actually measured 
values . 

15 

( 4 ) Fourth Embodiment 

With reference to the drawings, the fourth embodiment 

of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a pesticidal 
2 0 compound in accordance with this embodiment estimates the 

behavior of the compound in the case where the floor in the 

third embodiment is constituted by a rug having ears of fiber. 

For convenience of explanation, members identical to those 

shown in the drawings of the above-mentioned third 
25 embodiment will be referred to with numerals or letters 

identical thereto, without their explanations provided. 
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Explained here is a case where the kind of floor is a carpet. 
Accordingly, only the portion different from the third 
embodiment will be explained. 

First, the carpet is modeled as shown in Figs. 13A, 
5 13B, and 13C. Namely, the carpet is divided into a plurality 

of carpet fibers (ears of fiber) planted on a substrate and 
space portions existing between the fibers . When a chemical 
is spatially sprayed in a room in which such a carpet covers 
the whole floor surface, as particles fall on the carpet, 

10 the chemical infiltrates into the carpet fiber and 

accumulates in the space portion. 

While the method of simulating an indoor behavior of 
a pesticidal compound in accordance with this embodiment 
also conforms to the flowchart of Fig. 1, at SI, as the 

15 primary condition, the carpet-related primary condition 

shown in Table 23 is added to the primary condition of the 
third embodiment {see Tables 16 to 20). Here, the 
cross-sectional area (A c ) of the space portion of the carpet 
refers to the area in which the space portion is in contact 

20 with the indoor air (see Figs. 13A, 13B, and 13C). 
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The examples of values in Table 23 are based on the 
values shown in Figs. 13A, 13B, and 13C. Here, one space 
portion is set to a size of 0.1 [mm] x o.l [mm] * 3 [mm], 
and one square of carpet fiber is set to a size of 2.5 [mm] 
x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 [pieces] 
of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

A c = (0.1 [mm] x o.l [mm] x 16 [pieces] )/l [cm] 
= 1.6 x io -3 [mVm 2 floor] 

A c5 = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces] )/l 

[cm] 

= 0.192 [m 2 /m 2 floor] 
A 5 = 1 - A c 

= 1 - 1.6 x 10" 3 [m 2 /m 2 floor] 
A a5 = A 5 + A cS 

= 1.19 [m 2 /m 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

As the secondary condition at S2 , the carpet-related 
secondary condition shown in Table 24 is added to the 
secondary condition of the third embodiment. Here, k 5 and 
V 5 , which are respectively transference speed of the 
compound in the floor and the volume of the floor in the 
third embodiment (see Table 22), are respectively replaced 
by the transference speed of the compound in the carpet fiber 



143 



SCC 97-06 



and the volume of the whole fiber of the carpet. 
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TABLE 24 



CARPET - RELATED SECONDARY CONDITION 


CVA /fT> AT 


DEFINITE FACTOR CONDITION 




EVAPORATION CONSTANT OF CHEMICAL IN SPACE PORTION 


Rd 


HUMIDITY AT TEMPERATURE (T d ) OF SPACE PORTION 


H m 


VOLUME OF WHOLE FIBER OF CARPET 


V5 


TRANSFERENCE SPEED OF COMPOUND IN PRODUCT 


k c 


PROPERTY CHANGE CONSTANT OF PRODUCT 


a 


TEMPORAL CHANGE CONDITION 




TRANSFERENCE COEFFICIENT OF COMPOUND 




BETWEEN CHEMICAL IN SPACE PORTION AND AIR 


D c4 


TRANSFERENCE COEFFICIENT OF COMPOUND 




BETWEEN CHEMICAL IN SPACE PORTION AND CARPET FIBER 


D c5 


VOLUME OF CHEMICAL IN SPACE PORTION 


V c 


FUGACITY CAPACITY OF CHEMICAL IN SPACE PORTION 


Zc 
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A method of computing the definite factor conditions 
and temporal change conditions shown in the above-mentioned 
Table 24 will be explained hereinafter as exemplified by 
the case using an aqueous aerosol, i.e., case where the 
5 particle-constituting solvent is water. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T), and humidity (H), the evaporation 
constant (R d ) of the chemical entered into the space portion 
of the carpet is expressed as: 

R d =A c h(H m -H)lC H -(46) 
10 C H = 0.24 + 0.46// 

h = 3.06xHr 4 (r-r t J /3 

wherein T d is the temperature of the space portion of the 
carpet (corresponding to the wet-bulb temperature), and H m 
is the humidity {corresponding to the saturated humidity 
at the wet-bulb temperature) of the space portion at T d . 

15 Here, T d and H B can be determined from the above-mentioned 

"mass-based humidity table" according to the room 
temperature (T) and humidity (H). In the case of an 
oil-based product (aerosol), while the vapor constant of 
the chemical can be computed when properties of the water 

2 0 and oil are compared with each other and then are multiplied 

by R d , R d = 0 may be assumed in the case where the evaporation 
of the oil is slow. 

From the diffusion coefficient (D c ) and the area (A a5 ) 
in which the carpet fiber is in contact with the compound, 
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the volume (V 5 ) of the whole fiber of the carpet is determined 
as : 

V 5 =2^DjA a5 •••(47) 

The transference speed (k c ) of the compound in the 
product is expressed as k 4 /100 by use of k 4 (= GV 4 /A 4 ) which 
indicates the transference speed of the compound in the 
indoor air. 

From the above-mentioned transference speed (k c ) of 
the compound in the product, transference speed (k 4 ) of the 
compound in the indoor air, floor size (L 4 * W 4 ), cross- 
sectional area (A c ) of the space portion, Fugacity capacity 
(Z c ) of the chemical in the space portion determined by 
expression (50) which will be mentioned later, and Fugacity 
capacity (Z 4 ) of the air determined by expression (30) of 
the third embodiment, the transference coefficient (D c4 ) of 
the compound between the chemical in the space portion of 
the carpet and the air is expressed as: 

D c4 = —, r -(48) 

From the transference speed (k c ) of the compound in 
the product, diffusion coefficient (D c ), floor size (L 4 x 
W 4 ) , area (A cS ) in which the carpet fiber is in contact with 
the compound in the space portion, Fugacity capacity (Z c ) 
of the chemical in the space portion, and Fugacity capacity 
(Z 5 ) of the floor determined by expression (31) of the third 
embodiment, the transference coefficient (D c5 ) of the 
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compound between the chemical in the space portion and the 
carpet fiber is expressed as: 

D c5 = r^-r r • -(49) 

wherein the transference speed (k 5 ) of the compound within 
5 the carpet fiber is expressed as (D 0 /t) 0 - 5 . 

The volume (V c ) of the chemical in the space portion 
is determined from the floor size (L 4 * W 4 ), spray rate, 
content (C sd ) of the solvent (Xylene) contained in the 
product, dilution (X sol ) of the product, and evaporation 
10 constant (R d ) of the chemical in the space portion determined 

by expression (46). 

From the vapor pressure (P s ), water solubility (C s ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and property change constant (a) of the product, 
15 Fugacity capacity (Z c ) of the chemical in the space portion 

is expressed as: 

Z c = (e- a 'K ow + 1 - e at )C S IP S ■■■ (50) 

Here, the property change constant (a) of the product 
can be computed according to the evaporation constant (R d ) 

20 of the chemical in the space portion determined by expression 

(46), content (C sd ) of the solvent (Xylene) contained in the 
product, and dilution (X sol ) of the product. Namely, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is also expressed as: 

25 Z c = [ (ratio of organic solvent in product) K ow + (ratio 
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of water in product)] * C s /P s 

By simultaneously forming this expression and the 
above-mentioned expression (50), the property change 
constant (a) of the product is computed. Here, in the case 
of the oil-based product (aerosol), Fugacity capacity (Z c ) 
of the chemical in the space portion is given by expression 
(31) of the third embodiment from the vapor pressure (P s ) , 
water solubility (C s ), and octanol/water distribution 
coefficient (K ow ) of the compound. 

Consequently, the behavior of the compound in the 
space portion of the carpet is expressed in the form of 
differential equation concerning Fugacity (f c ) as: 

— v c z c = Ys y n i v A V i Z ifJ L n Deposition^) 

F - change + R d L 4 W 4 Z c f c V- change 

- D c4 (f c - f 4 ) Transference^) 

- D cs (fc - f 5 ) Transference^) 

- KV c Z c f e Degradation •••(51) 

Also, expressions (42) and (41) in the third embodiment are 
respectively replaced by: 

^V 5 Z 5 = 4o~TtA a5 Z 5 f 5 V- change 
dt 

F - change + y 5 (l - A c )^- vV i Z i f- t Deposition(i) 

- D c5 (f 5 - f c ) Transference(c) 

- D 45 (f 5 - f 4 ) Transference^) 

- K 5 V 5 Z 5 f 5 Degradation ■ ■ ■ (52) 
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dt 

F - change 



= -GV 4 Z 4 f 4 



Ventilation 




2 n ; Z) i4 {f 4 - fi) Transference^) 
5) D 4j if* - fj ) Transfer en ce(j) 



- K 4 V 4 Z 4 f 4 Degradation 

- D 4c (f 4 - f c ) Transfer en ce(c) 



-(53) 



In expression (52), terms of V-change, Deposition (i), 
Transference (c) , Transference (4), and Degradation 
respectively indicate the change in volume (increase over 
time) of the carpet fiber, attachment accompanying the 
falling of particles, amount of transference of the compound 
between the carpet fiber and the space portion, amount of 
transference of the compound between the carpet fiber and 
the air, and change in amount of photodegradation of the 
compound . 

In expression (53), terms of Ventilation, 
Transference (i), Transference (j), Degradation, and 
Transference (c) respectively indicate the amount of 
discharge of the compound outdoors, amount of transference 
of the compound between the air and the particles, amount 
of transference of the compound between the air and the floor 
(carpet fiber ) /ceiling/wall, change in amount of 
photodegradation of the compound, and amount of 
transference of the compound between the air and the space 
portion. 

Accordingly, Fugacity calculation at S3 is performed 
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such that eight kinds of differential equations of 
expressions (40), (43), and (44) in the third embodiment 
and the above-mentioned expressions (51) to (53) are solved 
by Runge-Kutta-Gill method. The subsequent steps of S4 to 
5 S17 are similar to those in the third embodiment. 

As explained in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
accordance with this embodiment, in the case where the kind 
of floor is a carpet having ears of fiber and a space portion 

10 between the ears, a differential equation for Fugacity in 

the above-mentioned space portion of the carpet is added 
to seven kinds of simultaneous differential equations 
concerning Fugacity in the three kinds of large, medium, 
and small particles, indoor air, floor (carpet fiber), 

15 ceiling, and wall, so as to calculate each Fugacity. 

Consequently, even in the case where the chemical is 
accumulated in the space portion, so that the compound 
behaves differently from the case with a flat floor, the 
behavior of the compound can be estimated accurately, 

2 0 whereby various kinds of simulations can be performed. 

(5) Fifth Embodiment 

With reference to the drawings, the fifth embodiment 
of the present invention will be explained in the following. 
25 The method of simulating an indoor behavior of a 

pesticidal compound (hereinafter simply referred to as 
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compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment estimates, 
for example, the behavior of the compound in the case where 
the chemical is heated to vaporize by use of an electric 
heating vaporizer such as electronic mosquito repellent in 
order to terminate indoor vermin. Here, the chemical 
encompasses all kinds of insecticides including pyrethroid 
insecticidal compounds and organophosphorus insecticidal 
compounds . 

First, as shown in Fig. 15, an environment is 
classified into 11 kinds of media consisting of 3 kinds of 
condensed particles 9, 10, and 11, high-concentration air 
12, medium-concentration air 13, low-concentration air 14, 
floor 5, wall 6, high-concentration ceiling 15, medium- 
concentration ceiling 16 , and low-concentration ceiling 17 . 

Namely, in the case where the chemical is vaporized, 
though the whole compound is initially evaporated as 
complete vapor (high-concentration air 12), the condensed 
particles 9, 10, and 11 are successively formed when the 
emission rate of the compound (E T ) in the medium- 
concentration air 13 exceeds the saturated amount of the 
compound in the vapor. These condensed particles 9, 10, and 
11 are divided into three kinds depending on the time between 
the occurrence and extinction of particles. The indoor air 
other than the above-mentioned high-concentration air 12 
and medium-concentration air 13, which is mainly respired 
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by humans, is defined as the low-concentration air 14. The 
high-concentration ceiling 15 is a ceiling portion into 
which the condensed particles 9, 10, and 11 are absorbed; 
the medium-concentration ceiling 16 is a ceiling portion 
in contact with the medium-concentration air 13, and the 
low-concentration ceiling 17 is a ceiling portion in contact 
with the low-concentration air 14. 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned 11 kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is heated and vaporized, i.e., temporal 
change in compound distribution in the room, is estimated 
and analyzed. 

Here, since the floor 5 mainly comprises three kinds, 
i.e., tatami, flooring, and carpet (rug), the differential 
equation is formed so as to correspond to the kind of floor. 
Explained in this embodiment is the case where the floor 
5 is made of tatami. 

With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 
be mainly divided into a step of dividing an indoor 
environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2); a step of determining the fugacity 
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of the compound in each medium from the differential equation 

(53) ; a step of determining the indoor behavior of the 
compound from the fugacity of the compound in each medium 

(54) ; a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S3 5, 
see Fig. 5) ; and a step of evaluating safety of the compound 
with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step (hereinafter abridged as S) 1, a primary 
condition is inputted. The primary condition comprises 
physicochemical properties of the compound (see Table 25), 
indoor environmental behavior properties of the compound 
(see Table 26), indoor environment (see Table 27), and 
product properties (see Table 28). 
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INDOOR ENVIRONMENT | 


ROOM SIZE 
FLOOR AREA 
TEMPERATURE 
RELATIVE HUMIDITY 
ABSOLUTE HUMIDITY 

VENTILATION RATE 

INDOOR VAPOR PRESSURE IN PARTICLE - CONSTITUTING SOLUENT 
AIR DIFFUSION CONSTANT 
OIL COMPONENT CONTENT 


WET - BULB TEMPERATURE 
VAPOR PRESSURE OF WATER 


OIL 
BASE 


WATER 
BASE 
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Here, A 12 , A 13 , and A 14 in Table 28 indicate cross- 
sectional areas of the respective airs 12, 13, and 14 with 
respect to the direction of movement of indoor air (arrow 
in Fig. 15) . 

The examples of values in Tables 25 to 28 assume those 
of a six -mat room (9.72 [m 2 ] ) of a typical apartment in Japan 
in summer, using an electric heating vaporizer employing 
an oil-based chemical (see Figs. 16A and 16B) . The electric 
heating vaporizer comprises a core portion of 1 [cm] x 1 
[cm] (at a temperature of 405 [K] ) , a heater portion of 4 
[cm] x 4 [cm] x 0.5 [cm] (at a temperature of 373 [K]), and 
a container in which a chemical is packed (at a temperature 
of 3 03 [K] ) . The example of value of compound content (R a ) 
in Table 28 indicates the volume of compound relative to 
the volume of particles immediately after condensation. 

Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2) . The 
secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 29) determined by the 
primary condition alone and a temporal change condition 
accompanying a temporal change (see Table 30). Their 
calculations will be explained later in detail. 
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TABLE 29 



DEFINITE FACTOR CONDITION 


SYMBOL 


FLOATING NUMBER OF CONDENSED PARTICLES 
EVAPORATION CONSTANT OF CONDENSED PARTICLES 
EMISSION AMOUNT OF CHEMICAL IN GAS STATE 
DIFFUSION COEFFICIENT OF COMPOUND IN AIR 
SURFACE AREA OF FLOOR / WALL / CEILING 
DIFFUSION COEFFICIENT OF FLOOR/ WALL / CEILING 
AIR VELOCITY 


ni (i=9,10,ll) 
a 
E v 
D ca 

Ak(k=5,6,15,16,17) 

D c 

Vf 


FUGACITY CAPACITY 

OF CONDENSED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Zi (i=9,10,H) 
Zj 0=12,13,14) 
Zk(k=5,6,15,16,17) 


DEGRADATION CONSTANT 

OF CONDENSED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Ki (i=9,10,ll) 
Kj 0=12,13,14) 
Kk(k=5,6,15,16,17) 
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TABLE 30 



TEMPORAL CHANGE CONDITION 


SYMBOL 


CONDENSED PARTICLE DIAMETER 


j. /:_q in ill 


TRANSFERENCE SPEED OF COMPOUND 
1JN CUJNDEJNbJbD FAKHULli 
AIR 

FLOOR / WALL / CEILING 


i^. o 1 o 1 1^ 

Kj \1— y,LV, 11) 

kj (j=12,13,14) 
kk(k=5,6, 15, 16,17) 


TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN CONDENSED PARTICLE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN AIR AND FLOOR / WALL / CEILING 


Dii3 0=9,10,11) 

Dj k 0=12,13,14, 
k=5,6,15,16,17) 


VOLUME OF CONDENSED PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


Vi (i=9,10,ll) 

Vk (k=5,6, 15,16,17) 
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Using the secondary condition, 11 kinds of Fugacity 
are calculated (S3). Namely, 11 kinds of differential 
equations concerning the 3 kinds of condensed particles 9, 
10, and 11, 3 kinds of air 12, 13, and 14, floor 5, wall 
6, and 3 kinds of ceilings 15, 16, and 17 are simultaneously 
formed and are solved by Runge-Kutta-Gill method, whereby 
11 kinds of Fugacity are computed over time. Here, an 
estimation nick time width (minute time unit) set when 
solving the differential equations is automatically set so 
as to be varied in response to a fluctuation in mass balance 
of the compound. 

At S4, using thus computed Fugacity (f 14 ) of the 
low-concentration air 14, temporal concentration of the 
compound in the air is computed; whereas temporal residual 
amounts of the compound on the floor 5 is computed by use 
of Fugacity (f 5 ) of the floor 5. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 
an estimated exposure amount in inhalation indicating a 
degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 
estimated exposure amount in inhalation (S7). At S8, the 
inhalation safety coefficient is compared with a reference 
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value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the floor 5 to which the chemical is attached 
(S9). Thereafter, at S10, it is judged whether to perform 
a percutaneous safety evaluation or not. If the 
percutaneous safety evaluation is to be performed, a 
percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 
condition is considered. 

If the percutaneous safety evaluation is not selected 
at S10, calculated at S13 according to the estimated 
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percutaneous exposure amount is an estimated oral exposure 
amount indicating the degree of exposure in the case where 
the chemical attached to a hand or the like is taken orally 

(513) . Subsequently, according to the estimated oral 

5 exposure amount, an oral safety coefficient is computed 

(514) . This oral exposure may occur, in particular, when 
an infant puts a chemical-attached hand into the mouth. At 
S15, as with the safety evaluation in inhalation, the oral 
safety coefficient is compared with a reference value 

10 defined in each country, whereby the safety is evaluated. 

If there is a problem in safety at S15, the operation returns 
to SI, where alteration of the primary condition is 
considered. 

Finally, if it is judged to be safe at S8, S12, or S15, 
15 overall safety is evaluated. Here, the sum of respective 

reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
20 (S16). This overall safety coefficient is compared with a 

reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4, S6, S7, S9, Sll, S13, and S14 will be explained in detail. 
25 (i) Secondary Condition Calculation (S2) 

Calculations of the definite factor condition shown 
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in Table 29 and the temporal change condition shown in Table 
3 0 will be explained. 

The condensed particles 9, 10, and 11 are successively- 
generated during time t 1 and are absorbed by the high- 
concentration ceiling 15. Specifically, the condensed 
particles 9 exist during the period of time t r /3 after 
evaporation, the condensed particles 10 are generated after 
the lapse of time t 1 /3, and the condensed particles 11 are 
generated after the lapse of time 2t 1 /3. Here, the floating 
number (n if i = 9, 10, 11) of the condensed particles is 
constant from their generation till absorption, and is 
expressed, from the specific gravity ( p ) of the compound, 
compound content (RJ , emission rate (E T ) , condensation 
ratio (C r ), and initial diameter (d 0 ) of condensed particles 
9, 10, and 11, as follows: 



wherein G is the ventilation rate. 

According to surface temperature (T d ) of particles, 
properties of particle-constituting solvent (vapor 
pressure P d/ molecular weight M d , and specific gravity p 
d ), and indoor environment (air diffusion coefficient D air , 
vapor pressure P°o, and room temperature too), the 
evaporation constant ( a ) of the condensed particles 9, 10, 
and 11 is defined as follows: 




...(54) 
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a — 



4 D air M c 




...(55) 



wherein R is a gas constant. 

Here, the upper parts of the above-mentioned Tables 
27 and 28 exemplify the case where the heat-vaporized aerosol 
is of oil base. In this case, Too (temperature at the site 
far from particles) and T d are set to room temperature (T) , 
whereas poo (vapor pressure of oil at the site far from 
particles) is set to zero. 

In the case where the heat-vaporized aerosol is of 
water base, it is necessary for P d , M d , p Ar T d , and poo in 
expression (55) to be changed to values based on properties 
of water. Specifically, the above-mentioned factor values 
when the room temperature (T) is 298 [K] (= 25°C) and the 
relative humidity is 60[%RH] are exemplified in the lower 
parts of Tables 27 and 28. 

Here, T d can be determined, according to the room 
temperature (T) and relative humidity (H), from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 
Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
the other hand, P d and poo can be computed from the following 
expressions : 

log 10 P d = 10.23 - 1750/(T d -38) 

log 10 poo = 8.23 - 1750/(T-38) + log lo 0 
wherein 0 is an indoor relative humidity. 
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From the molecular weight (M a ) of the compound, 
emission rate (E T ), and condensation rate (C r ), the amount 
of emission of chemical (E v ) to the high-concentration air 
12 is expressed as follows: 
E v =E T (1-C r )/M a ...(56) 

The surface area (A k , k = 5, 6, 15, 16, 17) of each 
of the floor 5, wall 6, high-concentration ceiling 15, 
medium-concentration ceiling 16, and low-concentration 
ceiling 17 is determined by the room size (V 4 ) and floor area 
(L 4 x w 4 ). 

Using the room size (V 4 ), ventilation rate (G), and 
cross-sectional area (A 14 ) of the room with respect to the 
direction of movement of indoor air, the indoor air velocity 
(v f ) is determined as follows: 

V f =G-^- ...(57) 
A 14 

From the vapor pressure (P 3 ) and melting point (T M ) 
of the compound and surface temperature (T d ) of particles, 
Fugacity capacity (Z if i = 9, 10, 11) of the condensed 
particles 9, 10, and 11 is expressed as: 

7 _ 6x106 

P[RT d -(58) 
PI =P*cxp{6.79(T M /T d -l)} 

Though the vapor pressure of the compound in a liquid state 
(P L S ) is computed by use of T M and T d here; in the case where 
T M cannot be obtained, P L S may be set identical to the vapor 
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pressure of the compound in a solid state (P s ) . 

From the room temperature (T) , Fugacity capacity (Z j7 
j = 12, 13, 14) of the airs 12, 13, and 14 is expressed as 
follows : 

5 Z.=^- ...(59) 

RT V ' 

From the vapor pressure (P s ), water solubility (C s ), 
and octanol /water distribution coefficient (K ow ) of the 
compound, and from the oil component content ( p k , k = 5, 
6, 15, 16, 17) of the material constituting the floor 5, 

10 wall 6, and ceilings 15, 16, and 17, Fugacity capacity (Z k , 

k = 5, 6, 15, 16, 17) of the floor 5, wall 6, and ceilings 
15, 16, and 17 is expressed as: 
z * = P k K ow C s /P s ...(60) 

Mainly generated in the 11 kinds of media are 

15 degradation reactions caused by light and oxidation. The 

degradation constant (K if i = 9, 10, 11) of the condensed 
particles 9, 10, and 11, degradation constant <K jr j = 12, 
13, 14) of the airs 12, 13, and 14, and degradation constant 
(K k , k = 5 , 6, 15, 16, 17) of the floor 5 , wall 6 , and ceilings 

20 15, 16, and 17 are defined by use of the respective half-lives 

of degradation by light/oxidation (T £ ), (T^, and ( T k ) . 
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Ki -°^ (61) 

While T i andr k are often unavailable, K A and K k may be 
nullified in this case from the viewpoint of securing safety 
of inhabitants. WhenTj of the compound in the air is 
5 unavailable, it may be determined by use of Atmospheric 

Oxidation program (Atkinson et al., 1984, Chem. Rev. Vol. 
84, pp. 437-470) from the structural formula of the compound . 

Since the volume of particles decreases over time due 
to evaporation of the particle-constituting solvent 

10 component, during the evaporation of particle-constituting 

solvent component, using the diameter (d 0 ) of particles 
immediately after spraying (t = 0) and the evaporation 
constant ( a) determined from expression (55) , the diameter 
(d if i = 9, 10, 11) of the condensed particles 9, 10, and 

15 11 at time t is represented as: 

d. = p 0 2 -2cd (62) 

From the diameter (dj of condensed particles 9, 10, 
and 11 determined by expression (62), indoor velocity (v f ) 
determined by expression (57), and Fugacity capacity (Z L ) 

2 0 of condensed particles 9,10, and 1 1 determined by expression 

(58), and Fugacity capacity (Z 13 ) of the medium- 
concentration air 13 determined by expression (59), the 
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transference coefficient (D il3 ) between the condensed 
particles 9, 10, and 11 and the medium-concentration air 
13 becomes: 

D il3 = — (63) 

5 wherein A i (= is the surface area of particles, k 13 ( = 

v f + v £ ) is the transference speed of the compound in the 
medium-concentration air 13, and k L (= v f /100) is the 
transference speed of the compound in the particles. Here, 
v L is the floating speed of the condensed particles 9, 10, 
10 and 11, and can be determined by use of hydrodynamics since 

it is under the control of the flow caused by the heat of 
the electric heating vaporizer, gravity, and air 
resistance. 

From the half-life of transference to air ( T jk ) in the 
15 ceilings 15, 16, and 17, volume (V jf j = 12, 13, 14) of the 

airs 12, 13, and 14, volume (V k ) of the floor/wall/ceiling 
determined by expression (66) mentioned later, Fugacity 
capacity (Z^) of airs 12, 13, and 14 determined by expression 
(59), and Fugacity capacity (Z k ) of floor /wall/ceiling 
20 determined by expression ( 60 ) , the transference coefficient 

(D jk ) of the compound between the medium-concentration air 
13 and the ceilings 15 and 16 or between the low- 
concentration air 14 and the floor 5, wall 6, and low- 
concentration ceiling 17 becomes: 
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0.693 

D. t = — i 1 (64) 

D. k = (64*) 

* l/(A ; .A. t Z.) + l/(M ; ,Z A ) 

Here, the total volume V 12 + V 13 + V 14 of the airs 12, 13, and 
14 equals the volume (V 4 ) of the room. 

The above-mentioned transference coefficient (D jk ) 
may also be determined by the expression represented by the 
above-mentioned expression (64 ' ) by use of the transference 
speed (kj) of the compound in the air and the transference 
speed (k k ) of the compound in the interior materials (floor 
5, wall 6, and ceilings 15, 16, and 17). 

Here, the half-life of transference to air (T jk ) in 
expression (64) is hard to actually measure and varies under 
various conditions, thus making it difficult to determine 
an accurate value thereof. By contrast, all the parameters 
in expression (64 ' ) can be easily determined by calculation 
without necessitating measurement. For example, using the 
diffusion coefficient (D c ) of compound, the transference 
speed (k k ) of the compound in the interior material is 
calculated as (D c /t)°- s ; whereas A jk can be determined from 
the contact area between the individual airs 12, 13, and 
14 and the individual interior materials (floor 5, wall 6, 
and ceilings 15, 16, and 17). Here, D c can be determined 
when the diffusion coefficient (D ca ) of the compound in the 
air is multiplied by 10" 9 , whereas D oa can be computed from 
properties of the compound (e.g., structural formula, 
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molecular weight, and the like) according to Wike and Lee 
method ( "Handbook of Chemical Property Estimation Methods , " 
McGraw-Hill Book Company, 1982), for example. Accordingly, 
the transference coefficient (D jk ) is more preferably 
determined by expression (64') than by expression (64). 

Using the diameter {di) of particles determined by 
expression (62) and the evaporation constant (a), change 
in volume (dV i /d J t) of the condensed particles 9, 10, and 11 
is expressed by: 

^—Z-ad, (65) 

Namely, as the particle-constituting solvent component 
evaporates over time, the volume (V i7 i = 9, 10, 11) of 
condensed particles 9, 10, and 11 decreases. 

Assuming that each of the floor 5, wall 6, and ceilings 
15, 16, and 17 before heat-vaporizing the chemical is like 
a thin film and that, as the chemical infiltrates into the 
film, its thickness increases so as to enhance the volume 
thereof; from the surface area (A k ) of the 
floor /wall/ceiling and diffusion coefficient (D c ) of the 
compound in the floor /wall/ceiling , the volume (V k , k = 5, 
6, 15, 16, 17) of floor/wall/ceiling becomes: 
V k = 2^DjA k (66) 

(ii) Fugacity Calculation (S3) 

The behavior of the compound in the condensed 
particles 9, 10, and 11 is expressed in the form of 
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differential equation concerning Fugacity (fj as: 

4Nz ; = ^ad&fi ~ AaC/,-/*) - Wj/i -( 67 ) 
dt 2 

F - change V - change Transference^) Degradation 

Here, the terras of v-change, Transference (13), and 
Degradation respectively indicate the volume change 
5 (decrease over time) of the condensed particles, amount of 

transference of the compound between the condensed 
particles 9, 10, and 11 and the medium-concentration air 
13, and amount of degradation of the compound by light and 
oxidation in the condensed particles 9, 10, and 11 (K^Zi 

10 indicating the degradation speed quantity) . Here, when the 

condensed particles 9, 10, and 11 are absorbed by the 
high-concentration ceiling 15, Fugacity (f t ) immediately 
disappears and is added to Fugacity (f 15 ) explained later. 

The behavior of the compound in the high-concentration 

15 air 12 is expressed in the form of differential equation 

concerning Fugacity (f 12 ) as: 

^V l2 Z^ = E v + v f A i2 (Z l3 f l3 -Z 12 f l2 ) - K l2 V l2 Z l2 f l2 ...(68) 
at 

F - change Emission Transference^) Degradation 

Here, the terms of Emission, Transference (j), and 
Degradation respectively indicate the amount of emission 
20 of the compound as complete vapor, amount of transference 

of the compound with respect to the medium-concentration 
air 13, and amount of degradation of the compound by light 
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and oxidation (K 12 V 12 Z 12 indicating the degradation speed 
quantity) . 

The behavior of the compound in the medium- 
concentration air 13 is expressed in the form of differential 
5 equation concerning Fugacity (f 13 ) as: 

^f-v l3 z l3 = v f K 2 z 12 / 12 - (A 2 + As)ZoA. + AJufJ ■■■(&) 

at 

F - change Transference^) 

- ! nPM» ~ ft) ~ J 5 A 3 ,(/o " ft) 
Transference^) Transference (k) 

-K 13 V B Z 13 f i3 
Degradation 

Here, the terms of Transference ( j ) , Transference (i), 
Transference (k) , and Degradation respectively indicate the 
amount of transference of the compound with respect to the 

10 high-concentration air 12 and low-concentration air 14, 

amount of transference of the compound with respect to the 
condensed particles 9, 10, and 11, amount of transference 
of the compound with respect to the ceilings 15 and 16, and 
amount of degradation of the compound by light and oxidation 

15 (K 13 V 13 Z 13 indicating the degradation speed quantity). 

The behavior of the compound in the low-concentration 
air 14 is expressed in the form of differential equation 
concerning Fugacity (f 14 ) as: 
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.(70) 



F - change Transference^) & Ventilation 



-2D uk (f l4 -f k )-D Ml7 (f l4 -f 17 ) 



Transference(k) Transfer enceQJ) 




Here, the terms of Transference (j) & Ventilation, 
Transference (k), Transference (17), and Degradation 
respectively indicate the amount of transference of the 
compound with respect to the medium-concentration air 13 
and the amount of discharge thereof outdoors, amount of 
transference of the compound with respect to the floor 5 
and wall 6, amount of transference of the compound with 
respect to the low-concentration ceiling 17, and amount of 
degradation of the compound by light and oxidation (K 14 V 14 Z 14 
indicating the degradation speed quantity). 

The behaviors of the compound in the floor 5 and wall 
6 are respectively expressed in the forms of differential 
equations concerning Fugacity (f 5 ) and (f 6 ) as follows: 




•••(71) 



F - change V - change 



Transfer en ce (1 4) 
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—V 6 Z 6 = -jDjiAZ 6 f 6 -D 146 (f 6 - / 14 ) -(72) 
dt 

F - change V - change Transference^^) 

-K 6 V 6 Z 6 f 6 
Degradation 

Here, in each expression, the terms of V-change, 
Transference (14), and Degradation respectively indicate 
the volume change (increase over time) of the floor 5 or 
5 wall 6, amount of transference of the compound with respect 

to the floor 5 or wall 6, and amount of degradation of the 
compound by light and oxidation (K 5 V 5 Z 5 or K 6 V 6 Z 6 indicating 
the degradation speed quantity) . 

The behaviors of the compound in the high- 
10 concentration ceiling 15, medium-concentration ceiling 16, 

and low-concentration ceiling 17 are respectively expressed 
in the forms of differential equations concerning Fugacity 
(fis)/ ( f i 6 )/ and (f 17 ) as follows: 

^rv l5 z l5 = -4b~rtA i5 z l5 f l5 - D a 15 (/ 15 - / 13 ) -(73) 

at 

F - change V - change Transfer enceQS) 

~ -^15^15^15/15 
Degradation 
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^V 16 Z 16 = -jDJtA 16 Z l6 f 16 - D 13 " /13) 



•••(74) 



dt 



F - change V - change TransferenceQ3) 




^V l7 Z l7 = -jDJtA^U -D l4 17 (/ 17 - / 14 ) 
at 

F - change V - change TransferenceQ.A) 



-.(75) 




Here, in each expression, the terms of V-change, 
Transference (13), Transference (14), and Degradation 
respectively indicate the volume change (increase over 
time) of the ceilings 15, 16, and 17, amount of transference 
of the compound between the high-concentration ceiling 15 
or medium-concentration ceiling 16 and the medium- 
concentration air 13 , amount of transference of the compound 
between the low-concentration ceiling 17 and the low- 
concentration air 14, and amount of degradation of the 
compound by light and oxidation (K 15 V 1S Z 1S , K 16 V 16 Z 16 , or K 17 V 17 Z 17 
indicating the degradation speed quantity). 

When the following Fugacity (f a ) is added to the floor 

(k = 5), wall (k = 6), and ceiling (k = 15, 16, 17) by a 
predetermined ratio, accuracy can be improved. 

L=n.y, a ZJ : JV k /Z k --(76) 

Here, V tia and f ia respectively indicate the volume and 
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Fugacity of particles immediately before the condensed 
particles attach to the interior materials. 

The above-mentioned 11 kinds of differential 
equations (67) to (75) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
(f 5/ f 6 , f 9 to f 17 ). 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
width (dt) which is a minute time unit. Namely, the 
estimation nick time width is used such that solutions of 
the simultaneous differential equations are initially 
determined at a time (t 0 ), and then solutions of the 
simultaneous differential equations are determined at a 
time (t 0 + dt) to which the estimation nick time width is 
added. As solutions are obtained while estimation nick time 
widths are successively added, temporally changing Fugacity 
can be determined. Theoretically, as the set time of the 
estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 
calculation time. By contrast, when the set time is too long, 
solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 
large change. 

Specifically, mass balance is always confirmed such 
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that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 
balance starts fluctuating. For example, when the 
fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 

compound input amount/ (existing amount + degrading 
amount + discharging amount) = 0.95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is 
determined by the compound content (R a ) and emission rate 
(E T ) . Since the temporal amounts of compound in the 11 kinds 
of media are determined by the simultaneous differential 
equations, they are summed up so as to compute the existing 
amount as shown in the following: 

existing amount = ^n.VZJ, + £ f k V k Z k + £ f k V k Z k 

i-9 A-S it-12 

degrading amount = £ n.Ky.^J. t + £ KJ k V k Z k + £ K k f k V k Z k I —(77) 

i-9 k-5 A- 12 

discharging amount = v y A u Z l4 / 14 

The degrading amount and discharging amount are as 
represented above. 

With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
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explained. 

First, an initial value of estimation nick time width 
( dt ) is inputted ( S2 1 ) . Then , an upper limit set value (e.g., 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 
is the lower limit of difference in mass balance are inputted 
(S22). Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27) . If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 
(or t 0 ) and the mass balance at t = t + dt (or t 0 + dt) is 
at least the upper limit set value or not (S28). If it is 
judged to be at least the upper limit set value at S2 8, then 
solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S28, it is judged 
at S30 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S30, since solutions are not 
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influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width (dt) or not (S32). When the estimation nick time width 
(dt) is not greater than the maximum value at S32, since 
solutions do not diverge, the estimation nick time width 
(dt ) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S3 2, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S3 0, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
are continued without changing the estimation nick time 
width (dt). 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S27 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
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lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
5 balance fluctuation, Fugacity can be computed accurately 

and efficiently. 

(iii) Computation of Temporal Concentration in Indoor Air 
and Floor Residual Amount (S4) 

The temporal concentration of the compound in the 
10 indoor air is computed when Fugacity (f 14 ) of the low- 

concentration air 14, determined by the above-mentioned 
item(ii) , is multiplied by its Fugacity capacity (Z 14 ) . The 
concentration in the low-concentration air 14 is used 
because the most part in the room is filled with the 
15 low-concentration air 14, and humans are most likely to 

inhale the low-concentration air 14. The concentration in 
the air of other concentrations may also be used when 
necessary . 

The floor residual amount of the compound is computed 
20 when Fugacity (f 5 ) of the floor is multiplied by Fugacity 

capacity ( Z 5 ) . 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 
25 indoor air forms a curve shown in Fig. 3A, for example. This 

concentration curve is integrated, an accumulated 
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concentration of the compound during a specific period (t x 
to t 2 ) is determined (see Fig. 3B), and the mean concentration 
in the indoor air is computed from thus determined value. 
While an arbitrary period is set as the specific period 
depending on the object, an appropriate period is usually 
set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 
concentration in the indoor air, amount of respiration, and 
exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] x amount of 
respiration [m 3 /kg/min] x exposure time [min/day] 
is effected. Here, as the above-mentioned amount of 
respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 

The inhalation safety coefficient is computed from a 
non-influential amount concerning inhalation toxicity 
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examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 

inhalation safety coefficient = 
5 inhalation non-influential amount [mg/kg/day]/ 

estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned floor residual amount forms a 

10 curve such as that shown in Fig. 4A, for example. This 

residual amount curve is integrated, the accumulated 
residual amount of the compound during a specific period 
(tj to t 2 ) is determined (see Fig. 4B) , and the mean residual 
amount is computed from thus determined value. While an 

15 arbitrary period is set as the specific period depending 

on the object, an appropriate period is usually set in view 
of the method of use of the product and the test period of 
toxicity data. 

Then, according to the mean floor residual amount, 

2 0 skin attachment ratio, contact area, and body weight, the 

estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean residual amount [mg/m 2 ] * skin attachment ratio [%] 

25 x contact area [ m 2 /day ]) /body weight [kg] 

is effected. Here, as the contact area, a published value 
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(e.g., 4 [rnVday]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor 5 where the compound 
exists. As this value, a published value or a value 
experimentally obtained from a model may be used. 

A model experiment method for the skin attachment 
ratio is as follows. A weight (8 cm * 8 cm x 8 cm; 4.2 kg) 
is placed on a denim cloth ( 8 cm x 10 cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants. 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 
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percutaneous non-influential amount [mg/kg/day ] / 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
are not many published values . Accordingly, a more accurate 
value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day ]/( estimated percutaneous 
exposure amount [mg/kg/day] x percutaneous absorption ratio 
[%]> 

Here, when the percutaneous absorption ratio is unknown, 
employed is a national guideline {e.g., 10%) which usually 
exists . 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
obtained in the above-mentioned item (v) , hand surface area 
ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
estimated percutaneous exposure amount [mg/kg/day] * hand 
surface area ratio [%] x oral transference ratio [%] 
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is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 
transference ratio is a hypothetical value, which is set 
to 100%, for example. 

In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
amount from tableware or food to mouth be added to the 
estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 
the amount of the harmful material remaining in tableware, 
tableware use area which is the sum of tableware surface 
areas, and oral transference ratio from tableware, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] x tableware use area 
[mVday] x oral transference ratio [%]/body weight [kg] 
is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio) . As the tableware contamination ratio, 
an actually measured value (e.g., 9%) or a hypothetical value 
may be used. 

The oral safety coefficient is computed from the 
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non- influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 
as : 

oral safety coefficient = 

oral non-influential amount [mg/kg/day]/ estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
this embodiment, the environment is divided into the 
condensed particles 9, 10, and 11 that are classified into 
three kinds depending on occurrence and extinction, the 
high-concentration air 12, medium-concentration air 13, and 
low-concentration air 14 that are classified into three 
kinds depending on the compound concentration, the floor 
5, the wall 6, and the high-concentration ceiling 15, 
medium-concentration ceiling 16, and low-concentration 
ceiling 17 that are classified into three kinds depending 
on the compound concentration; and differential equations 
concerning Fugacity in these members are simultaneously 
formed and are solved, whereby the indoor behavior of the 
compound when heat-vaporizing a chemical containing the 
compound is estimated- Here, when solving the simultaneous 
differential equations, the estimation nick time width is 
variably set, while constantly confirming mass balance of 
the compound indoors after the spraying, so that the amount 
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of input of the chemical indoors and the resulting solution 
coincide with each other. 

Accordingly, since mass balance of the compound after 
the spraying is always monitored such that the amount of 
5 input of the compound indoors and the resulting solution 

coincide with each other, thereby varying the estimation 
nick time width; the estimation nick time width is set longer 
when the mass balance fluctuates a little, whereas it is 
set shorter when the mass balance starts fluctuating greatly . 

10 Namely, when solving simultaneous differential equations 

including a parameter accompanying temporal change, the 
estimation nick time width is automatically set in response 
to the fluctuation in mass balance. Consequently, when 
processed by a computer, an accurate solution can be obtained 

15 in a short time. 

The method of evaluating safety of a pesticidal 
compound in accordance with this embodiment uses the 
estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 

20 is heated to vaporize. 

Accordingly, the safety of the compound with respect 
to the human body can be evaluated accurately in a short 
time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 

2 5 easily repeated while changing conditions, thereby making 

it easier to formulate a chemical having a high safety 
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conforming to the aimed object. 

Though the kind of the floor 5 is assumed to be tatami 
in this embodiment, differential equations similar to those 
mentioned above may also be formed in the case of flooring. 
5 In this case, however, when determining fugacity capacity 

(Z s ) of the floor 5 in equation (60), it is necessary for 
the particle-constituting solvent component content (see 
Table 27) to be changed from that of tatami to that of 
flooring. 

10 Though Fugacity is determined by use of Runge- 

Kutta-Gill method in this embodiment, other methods may be 
used for solving differential equations. Runge-Kutta-Gill 
method, however, is preferably used since a program for the 
above-mentioned differential equations can be easily made 

15 by Basic. Also in the case where differential equations are 

solved by a method other than Runge-Kutta-Gill method, 
similar effects can be obtained when the estimation nick 
time width is set as mentioned above. 



20 (6) Sixth Embodiment 

With reference to the drawings, the sixth embodiment 

of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a pesticidal 

compound in accordance with this embodiment estimates the 
25 behavior of the compound in the case where the floor in the 

fifth embodiment is constituted by a rug having ears of fiber. 
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For convenience of explanation, members identical to those 
shown in the drawings of the previous embodiment will be 
referred to with numerals or letters identical thereto, 
without their explanations provided. Explained here is a 
case where the kind of floor is a carpet. Accordingly, only 
the portion different from the fifth embodiment will be 
explained. 

First, the carpet is modeled as shown in Figs. 13A, 
13B, and 13C. Namely, the carpet is divided into a plurality 
of carpet fibers (ears of fiber) planted on a substrate and 
space portions existing between the fibers . When a chemical 
is heated to vaporize in a room in which such a carpet covers 
the whole floor surface, as condensed particles fall on the 
carpet, the chemical infiltrates into the carpet fiber and 
accumulates in the space portion. 

While the method of simulating an indoor behavior of 
a pesticidal compound in accordance with this embodiment 
also conforms to the flowchart of Fig. 1, at SI, as the 
primary condition, the carpet-related primary condition 
shown in Table 31 is added to the primary condition of the 
fifth embodiment (see Tables 25 to 28). Here, the 
cross-sectional area (A c ) of the space portion of the carpet 
refers to the area in which the space portion is in contact 
with the indoor air (see Figs. 13A, 13B, and 13C). 
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The examples of values in Table 31 are based on the 
values shown in Figs. 13A, 13B, and 13C. Here, one space 
portion is set to a size of 0.1 [mm] * 0.1 [mm] x 3 [mm], 
and one square of carpet fiber is set to a size of 2.5 [mm] 
x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 [pieces] 
of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

A c = (0.1 [mm] x o.l [mm] x 16 [pieces] )/l [cm] 
= 1.6 x io -3 [m 2 /m 2 floor] 

A c5 = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces ])/l 

[cm] 

= 0.192 [m 2 /m 2 floor] 
A s = 1 - A c 

= 1 - 1.6 x io -3 [m 2 /m 2 floor] 
A a5 = A 5 + A c5 

= 1.19 [m 2 /m 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

As the secondary condition at S2, the carpet-related 
secondary condition shown in Table 32 is added to the 
secondary condition of the fifth embodiment. Here, k s and 
V 5 , which are respectively transference speed of the 
compound in the floor and the volume of the floor in the 
fifth embodiment (see Table 30), are respectively replaced 
by the transference speed of the compound in the carpet fiber 
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and the volume of the whole fiber of the carpet. 
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A method of computing the definite factor conditions 
and temporal change conditions shown in the above-mentioned 
Table 32 will be explained hereinafter as exemplified by 
the case using an aqueous aerosol, i.e., case where the 
5 particle-constituting solvent is water. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T), and humidity (H), the evaporation 
constant (R d ) of the chemical entered into the space portion 
of the carpet is expressed as: 

R d =A c h{H m -H)IC H ••■(78) 
10 C H = 0.24 + 0.46tf 

A = 3.06xl0" 4 (r-7;j 1/3 

wherein T d is the temperature of the space portion of the 
carpet (corresponding to the wet-bulb temperature), and H m 
is the humidity (corresponding to the saturated humidity 
at the wet-bulb temperature) of the space portion at T d . 

15 Here, T d and H m can be determined from the above-mentioned 

"mass-based humidity table" according to the room 
temperature (T) and humidity (H). In the case of an 
oil-based product (aerosol), while the vapor constant of 
the chemical can be computed when properties of the water 

20 and oil are compared with each other and then are multiplied 

by R d , R d = 0 may be assumed in the case where the evaporation 
of the oil is slow. 

From the diffusion coefficient (D c ) and the area (A aS ) 
in which the carpet fiber is in contact with the compound, 
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the volume { V 5 ) of the whole fiber of the carpet is determined 
as : 

V 5 = 2jDjA a5 -(79) 

The transference speed (k c ) of the compound in the 
product is expressed as k 14 /100 by use of k 4 (= Gv 4 /A 14 ) which 
indicates the transference speed of the compound in the 
indoor air. 

From the above-mentioned transference speed (k c ) of 
the compound in the product, transference speed (k 14 ) of the 
compound in the indoor air, floor area {L 4 x W 4 ), cross- 
sectional area (A c ) of the space portion, Fugacity capacity 
(Z c ) of the chemical in the space portion determined by 
expression (82) which will be mentioned later, and Fugacity 
capacity (Z 14 ) of the low-concentration air determined by 
expression (59) of the fifth embodiment, the transference 
coefficient (D cl4 ) of the compound between the chemical in 
the space portion of the carpet and the low-concentration 
air is expressed as: 

D cl4 = 1 — ■ • • -(80) 

l/(kA c L 4 W 4 Z c ) + l/(k u A e L 4 W 4 Z u ) 

From the above-mentioned transference speed (k c ) of 
the compound in the product, diffusion coefficient (D c ), 
floor area (L 4 x W 4 ), area (A cS ) in which the carpet fiber 
is in contact with the compound in the space portion, 
Fugacity capacity (Z c ) of the chemical in the space portion, 
and Fugacity capacity (Z 5 ) of the floor determined by 
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expression (60) of the fifth embodiment, the transference 
coefficient (D o5 ) of the compound between the chemical in 
the space portion and the carpet fiber is expressed as: 

D = ...(81) 

l/(k e A eS L 4 W 4 Z e ) + l/(k 5 A c5 L 4 W 4 Z 5 ) V ' 

5 wherein k 5 (= (D c /t) 0 ' 5 ) is the transference speed of the 

compound within the carpet fiber. 

The volume (V c ) of the chemical in the space portion 
is determined from the floor area (L 4 x W 4 ), emission rate 
(E T ), condensation rate (C r ), and evaporation constant (R d ) 

10 of the chemical in the space portion determined by expression 

(78) . The volume (V c ) of the space portion decreases until 
the solvent in the chemical completely evaporates . 

From the vapor pressure (P s ) , water solubility (C 3 ), 
and octanol/water distribution coefficient (K ow ) of the 

15 compound, and property change constant (a) of the product, 

Fugacity capacity (Z c ) of the chemical in the space portion 
is expressed as follows: 
Z c =(e- at K ow +l-e"")C s IP* •••(82) 

Here, the property change constant (a) of the product 

2 0 can be computed according to the evaporation constant <R d ) 

of the chemical in the space portion determined by expression 
(78), content (C 3d ) of the solvent (Xylene) contained in the 
product, and dilution (X gol ) of the product. Namely, 
Fugacity capacity (Z c ) of the chemical in the space portion 

25 is also expressed as: 
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Z c = [ (ratio of organic solvent in product) K ow + (ratio 
of water in product)] x C s /P s 

By simultaneously forming this expression and the 
above-mentioned expression (82), the property change 
constant (a) of the product is computed. Here, in the case 
of the oil-based product (aerosol), Fugacity capacity (Z c ) 
of the chemical in the space portion is given by expression 
(57) of the fifth embodiment from the vapor pressure (P s ) , 
water solubility (C 3 ), and octanol/water distribution 
coefficient (K ow ) of the compound. 

Consequently, the behavior of the compound in the 
space portion of the carpet is expressed in the form of 
differential equation concerning Fugacity (f c ) as follows: 



Also, expressions (71) and (70) in the fifth 
embodiment are respectively replaced by the following 
expressions (84) and (85): 





Deposition(i) 



-(83) 



F - change 



+ R d L,W 4 ZJ c 

+ D cl3 (f c -f l3 ) 
-KV c ZJ c 



V - change 

Transference^) 
Transfernce(k) 
Degradation 
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—V 5 Z 5 = JDJtA 5 Z 5 f 5 V -change •••(84) 

dt 



-2* 



YiZJi Deposition^) 



+ D c5 (f s - f c ) Transference^) 
+ D 145 (f 5 - f 14 ) Transfernce(14) 
-K 5 V 5 Z 5 f s Degradation 

%y i4 Z 14 = GV u Z l4 f l4 Ventilation --(85) 

dt 

ii 

- ^ «;A 14 (/ 14 - f : ) Transference^) 

- gA4,-(/i4-/;) Transference(j) 

- K l4 V^Z l4 f lA Degradation 
~ A4 c ( /14 ~ fc ) Transfer en ce(c) 

Namely, as shown in expression (84), the differential 
equation concerning the floor (carpet fiber) includes, in 
addition to expression (71) of the fifth embodiment, the 
terms of Deposition (i) indicating the amount of deposition 
of condensed particles onto the carpet fiber and 
Transference (c) indicating the amount of transference of 
the compound between the carpet fiber and the space portion 

In expression (85), which is differential equation 
concerning the low-concentration air, terms of Ventilation , 
Transference (i), Transference (j), Degradation, and 
Transference (c) respectively indicate the amount of 
discharge of the compound outdoors, amount of transference 
of the compound between the low-concentration air and 
condensed particles, amount of transference of the compound 
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between the low-concentration air and the medium- 
concentration air and floor, amount of degradation of the 
compound by light and oxidation, and amount of transference 
of the compound between the low-concentration air and the 
space portion. 

Accordingly, Fugacity calculation at S3 is performed 
such that 12 kinds of differential equations of expressions 
(67) to (69) and (72) to (75) in the fifth embodiment and 
the above-mentioned expressions (83) to (85) are solved by 
Runge-Kutta-Gill method. The subsequent steps of S4 to S17 
are similar to those in the fifth embodiment. 

As explained in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
accordance with this embodiment, in the case where the kind 
of floor is a carpet having ears of fiber and a space portion 
between the ears, a differential equation for Fugacity in 
the above-mentioned space portion of the carpet is added 
to 11 kinds of simultaneous differential equations 
concerning Fugacity in 3 kinds of condensed particles, air 
divided into 3 kinds depending on the compound concentration, 
floor, wall, and ceiling divided into 3 kinds depending on 
the compound concentration, so as to calculate Fugacity. 

Consequently, even in the case where the chemical is 
accumulated in the space portion, so that the compound 
behaves differently from the case with a flat floor, the 
behavior of the compound can be estimated accurately, 
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whereby various kinds of simulations can be performed. 

(7) Seventh Embodiment 

With reference to the drawings , the seventh embodiment 
of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a 
pesticidal compound (hereinafter simply referred to as 
compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment estimates 
the behavior of the compound in the case where the chemical 
is sprayed over the whole floor surface, and is applied to 
the case where a carpet (rug) made of polymer fiber is spread 
on the whole floor surface. Here, usable as the chemical 
are pyrethroid insecticidal compounds, organophosphorus 
compounds, carbamate compounds, insect growth restrainers 
(IGR), and the like. Its solvent may be water or oil. 
Explained in this embodiment is the case where the solvent 
is water. The vermin to be killed is mainly flea, mite, and 
the like. 

First, the carpet is modeled as shown in Figs. 13A, 
13B, and 13C. Namely, the carpet is divided into a plurality 
of carpet fibers (ears of fiber) planted on a substrate and 
space portions existing between the fibers . When a chemical 
is sprayed over the whole floor area covered with a carpet, 
since the solvent of the chemical is water, the chemical 
infiltrates into the carpet fiber and accumulates in the 
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space portion. The chemical not only attaches to the floor 
but also is partly suspended in the air. 

Accordingly, as shown in Fig. 17, the indoor 
environment covered with the carpet is classified into seven 
5 kinds of media comprising the space portion of the carpet, 

two kinds of suspended large and small particles having 
different diameters, whole indoor air, floor (carpet fiber) , 
wall, and ceiling. 

The suspended particles are divided into two kinds 

10 since the behavior of the particles mainly depend on their 

diameter. Namely, while the individual particles are 
distributed in the same place immediately after spraying, 
as time passes, the large particles are distributed to a 
particle zone DZ 2 close to the floor due to their weight, 

15 and the small particles are distributed to a particle zone 

DZ 3 close to the ceiling since they are lighter than the large 
particles . 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned seven kinds 

2 0 of media are formed and are solved by Runge-Kutta-Gill method, 

whereby the indoor behavior of the compound in the chemical 
when the chemical is sprayed over the whole carpet surface, 
i.e., temporal change in the compound distribution in the 
room and amount of attachment to the interior materials of 

25 the floor, wall, ceiling, and the like, is estimated and 

analyzed. 
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With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 
be mainly divided into a step of dividing an indoor 
environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2 ) ; a step of determining the fugacity 
of the compound in each medium from the differential equation 

(53) ; a step of determining the indoor behavior of the 
compound from the fugacity of the compound in each medium 

(54) ; a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 
see Fig. 5) ; and a step of evaluating safety of the compound 
with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step ( hereinafter abridged as S ) 1, a primary 
condition is inputted. The primary condition comprises 
physicochemical properties of the compound (see Table 33), 
indoor environmental behavior properties of the compound 
(see Table 34) , indoor environment (see Table 35) , a carpet 
condition (see Table 36) , and product properties (see Table 
37). 
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The examples of values in Table 3 5 assume those of a 
six-mat room (9.72 [m 2 ]) of a typical apartment in Japan. 

The cross-sectional area (A c ) of the space portion of 
the carpet in Table 36 is an area in which the space portion 
is in contact with the indoor air (see Figs. 13A, 13B, and 
13C) . Further, the examples of values in Table 3 6 are based 
on the values shown in Figs. 13A, 13B, and 13C. Here, one 
space portion is set to a size of 0.1 [mm] x 0.1 [mm] x 3 
[mm] , and one square of carpet fiber is set to a size of 
2.5 [mm] x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 
[pieces] of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

A c = (0.1 [mm] x o.l [mm] x 16 [pieces] )/l [cm] 
= 1.6 x io~ 3 [m 2 /m 2 floor] 

A c5 = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces] )/l 

[cm] 

= 0.192 [m 2 /m 2 floor] 
A 5 = 1 - A x 

= 1 - 1.6 x 10' 3 [m 2 /m 2 floor] 
Ks = A s + A 1S 

=1.19 [m 2 /m 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2) . The 
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secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 38) determined by the 
primary condition alone, and a temporal change condition 
accompanying a temporal change (see Table 39). Their 
calculations will be explained later in detail. 



208 



TABLE 38 



DEFINITE FACTOR CONDITION 


SYMBOL 


EVAPORATION CONSTANT OF SUSPENDED PARTICLES 
VOLUME RATIO OF COMPOUND IN PRODUCT 
DIFFUSION COEFFICIENT OF COMPOUND IN AIR 
EVAPORATION CONSTANT OF CHEMICAL IN SPACE PORTION 


a 

Ra 
D ca 

Rd 


FUGACITY CAPACITY 

OF SUSPENDED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Zi (i=2,3) 

z 4 

Zj (j=5,6,7) 


DEGRADATION CONSTANT 

OF CARPET SPACE PORTION 
SUSPENDED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


K c 

Ki (i=2,3) 
K4 

Kj 0=5,6,7) 



209 



TABLE 39 



TEMPORAL CHANGE CONDITION 


SYMBOL 


SUSPENDED PARTICLE 
DIAMETER 
FALLING SPEED 
SUSPENDED NUMBER 


di 0=2,3) 
vi 0=2,3) 
ni 0=2,3) 


TRANSFERENCE SPEED OF COMPOUND 

BETWEEN CHEMICAL IN SPACE PORTION AND AIR 
BETWEEN CHEMICAL IN SPACE PORTION AND CARPET FIBER 
BETWEEN SUSPENDED PARTICLE AND AIR 
BETWEEN FLOOR / WALL / CEILING AND AIR 


Dc4 
D C 5 

D i4 0=2,3) 
D 4 j 0=5,6,7) 


VOLUME OF CHEMICAL IN SPACE PORTION 
VOLUME OF SUSPENDED PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


D c 

Di (i=2,3) 
Dj 0=5,6,7) 


FUGACITY CAPACITY OF CHEMICAL IN SPACE PORTION 


Zc 
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Using the secondary condition, seven kinds of Fugacity 
are calculated (S3). Namely, seven kinds of differential 
equations concerning the space portion of the carpet, two 
kinds of large and small suspended particles, indoor air, 
floor (carpet fiber), wall, and ceiling are simultaneously 
formed and are solved by Runge-Kutta-Gill method, whereby 
seven kinds of Fugacity are computed over time. Here, an 
estimation nick time width (minute time unit) set when 
solving the differential equations is automatically set so 
as to be varied in response to a fluctuation in mass balance 
of the compound. 

At S4, using thus computed Fugacity of the small 
particles (taking account of Fugacity of the large particles 
when necessary) and Fugacity of the indoor air, temporal 
concentration of the compound in the air is computed; whereas 
floor residual amount of the compound is computed by use 
of Fugacity of the carpet fiber and carpet space portion. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 
an estimated exposure amount in inhalation indicating a 
degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 
estimated exposure amount in inhalation (S7). At S8, the 
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inhalation safety coefficient is compared with a reference 
value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the floor to which the chemical is attached 
(S9). Thereafter, at S10, it is judged whether to perform 
a percutaneous safety evaluation or not. If the 
percutaneous safety evaluation is to be performed, a 
percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 
condition is considered. 

If the percutaneous safety evaluation is not selected 
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at S10, calculated at S13 according to the estimated 
percutaneous exposure amount is an estimated oral exposure 
amount indicating the degree of exposure in the case where 
the chemical attached to a hand or the like is taken orally 

(513) . Subsequently, according to the estimated oral 
exposure amount, an oral safety coefficient is computed 

(514) . This oral exposure may occur, in particular, when 
an infant puts a chemical-attached hand into the mouth. At 
S15, as with the safety evaluation in inhalation, the oral 
safety coefficient is compared with a reference value 
defined in each country, whereby the safety is evaluated. 
If there is a problem in safety at S15, the operation returns 
to SI, where alteration of the primary condition is 
considered. 

Finally, if it is judged to be safe at S8, S12, or S15, 
overall safety is evaluated. Here, the sum of respective 
reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
(S16). This overall safety coefficient is compared with a 
reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4, S6, S7, S9, Sll, S13, and S14 will be explained in detail. 

(i) Secondary Condition Calculation (S2) 
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Calculations of the definite factor condition shown 
in Table 38 and the temporal change condition shown in Table 
39 will be explained. 

According to surface temperature (T d ) of particles, 
5 properties of particle-constituting solvent (vapor 

pressure P d , molecular weight M d , and specific gravity p 
d ) , and indoor environment (air diffusion coefficient D air , 
vapor pressure P 00 , and room temperature T°o), the 
evaporation constant ( a ) of the suspended particles is 
10 defined as follows: 

a =^L(^.^] „. (86) 

R/° d \T d Tec) 

wherein R is a gas constant, and P d (= 2.2 6 * 10 3 [Pa]), M d 
(= 18 [g/mole]), and p d (= 1 * 10 s [g/m 3 ]) are assumed to 
be the values due to properties of water . Also, in this case, 

15 Too (temperature at the site far from particles) is set to 

room temperature (T), T d to indoor wet-bulb temperature ( = 
292 .5 [K] ) , and P°o (vapor pressure of water at the site far 
from particles) to 1 .89 x 10 3 [Pa]. These values are those 
with room temperature (T) of 298 [K] and relative humidity 

2 0 of 60 [%RH] . 

Here, T d can be determined, according to the room 
temperature (T) and relative humidity (H), from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 

25 Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
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the other hand, P d and P°o can be computed from the following 
expressions : 

logio^d = 10.23 - 1750/(T d -38) 

log 10 Poo = 8.23 - 1750/(T-38) + log lo 0 
wherein i/j is an indoor relative humidity. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T), and humidity (H), the evaporation 
constant (R d ) of the chemical entered into the space portion 
of the carpet is expressed as: 

R d = A c h(H m -H)IC H ...(87) 
C H = 0.24 + 0A6H 
h = 3.06xlO-\T-T d f !3 

wherein H m is a humidity at temperature (T d ) . The chemical 
entered into the space portion of the carpet reduces its 
volume upon evaporation and finally disappears, leaving the 
compound. Thus left compound infiltrates through the 
carpet fiber. 

From the surface temperature (T d ) of particles and the 
vapor pressure (P s ) of the compound, Fugacity capacity (Z if 
i = 2, 3) of the large and small particles is expressed as: 



P[RT d 

P L S =P s cxp{6.79(T M /T d -l)} 

Though the vapor pressure of the compound in a liguid state 
(P L S ) is computed by use of the melting point (T M ) of the 
compound and T d here; in the case where T M cannot be obtained, 
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P L S may be set identical to the vapor pressure of the compound 
in a solid state (P s ) . 

From the room temperature (T) , Fugacity capacity (Z 4 ) 
of the air is expressed as follows: 



5 




Since the carpet fiber and wallpaper covering the 
floor /wall/ceiling are formed by polymers; from the vapor 
pressure (P s ), water solubility (C s ), and octanol/water 
distribution coefficient (K ow ) of the compound, and from the 

10 oil component content ( p t ) of the carpet fiber and wallpaper 

covering the floor /wall/ceiling, Fugacity capacity (Z jr j 
= 5, 6, 7) of the floor/wall/ceiling is expressed as: 
Z j = p } K ow C s IP S ...(90) 

The degradation constants (K c ), (K if i = 2, 3), (K 4 ), 

15 (K jf j = 5, 6, 7) of the carpet space portion, suspended 

particles, air, and floor/wall/ceiling are defined by use 
of their respective half-lives of degradation by 
light/oxidation (T c ), (TJ, (T 4 )/ and ( r j) as follows: 

K e =0.693/t ' 
K. = 0.693 lx. 

I ...(91) 
K 4 = 0.693/t 4 [ V ) 

Kj = 0.693 It. 

20 Mainly generated in the seven kinds of media are 

degradation reactions caused by light and oxidation, and 
the reaction speed quantity at this time is represented by 
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the product of degradation constant, volume, and Fugacity 
capacity. Among the degradation constants shown in 
expression (91), the degradation constant (K 4 ) of the 
compound in the air is hard to actually measure, and it is 
5 not usually easy to determine this value. Accordingly, when 

the actually measured value is unavailable, it may be 
determined as computed by Atmospheric Oxidation program 
(Atkinson et al., 1984, Chem. Rev. Vol. 84, pp. 437-470) 
from the structural formula of the compound. Also, from the 

10 viewpoint of securing safety of inhabitants, the 

degradation constants in the particles and in the interior 
materials (floor/wall/ceiling) may be nullified when their 
actually measured values are unavailable. 

Since the volume of particles decreases over time due 

15 to evaporation, using the evaporation constant (a) 

determined from expression (86) and the diameter (d oi ) of 
suspended particle immediately after spraying (t = 0), the 
diameter (d ir i = 2, 3) of the suspended particles at time 
t is determined as : 

d. = Jdl - 2at ...(92) 
20 ' ' 

d< = y-Mo,- -(92') 

On the other hand, when the particle-constituting 
solvent in the suspended particle is completely evaporated, 
the diameter (d t ) becomes that of the compound alone, thereby 
being held constant. Namely, the particle diameter (d t ) 
2 5 after the complete evaporation of the particle-constituting 
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solvent is represented by the above-mentioned expression 
(92* ) by use of the volume ratio (R a ) of the compound in the 
product and the diameter (d oi ) of the particles immediately 
after spraying. 

From the compound content (C a ) and the specific gravity 
(P) of the compound, the volume ratio (R a ) of the compound 
in the product is expressed as C a /p. 

As shown in Fig. 18, immediately after spraying the 
chemical, a particle zone (DZjJ having a width (W 4 ), length 
(L 4 ), and space width (H (i) ) is formed at a certain height 
from the floor. In conformity to the falling speed (vj of 
suspended particles, the particle zone (DZ ± ) approaches the 
floor. When the bottom of the particle zone (DZi) reaches 
the floor, the compound in the suspended particles transfers 
to the carpet space portion and carpet fiber. 

Since movement of the particles is under the control 
of gravity and air resistance, according to Stokes ' law using 
the diameter (dj of suspended particles determined by 
expression (92), the falling speed (Vi, i = 2, 3) of the 
suspended particles is represented by the following 
expression ( 93 ) : 

v^B&Ldl-Pidl-lat) -(93) 
I877 

S =1 + [6.32 + 2.01exp(-8.322xl0 6 ^-)] ...(93') 

7.6x10 7 4 

v=(p/p,)M 2 -(93") 

wherein g is gravitational acceleration, 77 is the 
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coefficient of viscosity of the air, S c is a sliding 
correction coefficient , and (5 is a speed coefficient . Here, 
while particles do not conform to the Stokes ' law when they 
become small, the sliding correction coefficient (S c ) , which 
5 is a coefficient for correcting this phenomenon, is 

represented by the above expression (93* ) according to the 
diameter (dj of particles determined by expression (92). 

On the other hand, the falling speed (VjJ of the 
particles after the complete evaporation of particle- 

10 constituting solvent component is represented by the above 

expression (93") using the specific gravity (p) of the 
compound, specific gravity ( p d ) of the particle- 
constituting solvent, speed coefficient ( /? ) , and diameter 
(di) of particles determined by expression (92'). 

15 The number of suspended particles (n ir i = 2, 3) also 

decreases upon exchange of the air through a window and the 
like. Here, the small particles float in conformity to the 
indoor air flow, and the large particles, which reduce their 
diameter as their volume decreases, are assumed to conform 

20 to the flow of the indoor air flow. In this case, during 

the period till the particle zone (DZ L ) reaches the floor, 
the number of particles (n L ) would decrease only upon 
ventilation of the room. Accordingly, from the room size 
(V 4 = L„ x w 4 * h 4 ) , ventilation rate (G) , particle space width 

25 ( H a)' i = 2, 3), and particle falling speed (v £ ) determined 

by expression (93) or (93"), the particle number (n ± ) is 
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expressed as follows: 

^ = __A_y;_^ 4 ... (94 ) 
dt H (!) L 4 

wherein k 4 (= GV 4 /A 4 ) is the transference speed of the 
compound in the air. 

From the floor size (L 4 x W 4 ), ventilation rate (G) , 
cross-sectional area (A c ) of the space portion, Fugacity 
capacity (Z c ) of the chemical in the space portion determined 
by expression (101) mentioned later, and Fugacity capacity 
(Z 4 ) of the air determined by expression (89), the 
transference coefficient (D c4 ) of the compound between the 
chemical in the space portion and the air is represented 
as follows : 

D 4 ...(95) 

c l/(k c A c L 4 W 4 Z c )+l/(k 4 A c L 4 W 4 Z 4 ) 

wherein k c (= k 4 /100) is the transference speed of the 
compound on the water surface of the chemical in the space 
portion. 

From the diffusion coefficient (D p ) , floor size (L 4 
x w 4 ), ventilation rate (G), Fugacity capacity (Z c ) of the 
chemical in the space portion, and Fugacity capacity (Z 5 ) 
of the floor determined by expression ( 90 ) , the transference 
coefficient (D c5 ) of the compound between the chemical in 
the space portion and the carpet fiber is represented as: 

D s = ...(96) 

c5 l/(k c A c5 L 4 W 4 Z c ) + l/(k 5 A c5 L 4 W 4 Z 5 ) 
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wherein k 5 (= (D p /t) 0 - 5 ) is the transference speed of the 
compound in the carpet fiber, and A c5 is the contact area 
ratio of the carpet fiber and space portion relative to the 
floor. 

5 The diffusion coefficient (D p ) of the compound in the 

floor/wall/ceiling can be determined when the diffusion 
coefficient (D ca ) of the compound in the air is multiplied 
by 10" 9 , whereas D ca can be computed from properties of the 
compound (e.g., structural formula, molecular weight, and 

10 the like) according to Wike and Lee method ("Handbook of 

Chemical Property Estimation Methods," McGraw-Hill Book 
Company, 1982), for example. 

From the diameter (d A ) of particles determined by 
expression (92) or (92'), transference speed (k 4 ) of the 

15 compound in the air, falling speed (v £ ) determined by 

expression (93) or (93"), Fugacity capacity {Z L ) and (Z 4 ) 
determined by expressions (88) and (89), the transference 
coefficient (D i4 ) of the compound between the suspended 
particles and the air becomes : 

20 A 4 = -(97) 

1/(M-Z,-) +l/{(* 4 + V f )i4jZ 4 )} 

wherein Ai ( = TCd^ ) is the surface area of suspended particles , 
and ki (= k 4 /100) is the transference speed of the compound 
in the particle surface. 

From the diffusion coefficient (D p ) of the 
25 floor /wall/ceiling, room size (L 4 * W 4 x z 4 ) , ventilation rate 
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(G), and Fugacity capacity (Z 4 ) and (Z 3 ) determined by 
expressions (89) and (90), the transference coefficient 
(D 4j ) of the compound between the air and the 
floor /wall/ceiling becomes: 

D 4i = ...(98) 

41 l/ik.AjZ^ + l/ikjAjZj) K 

wherein A 3 is the surface area of the wall and ceiling, and 
kj (= (D p /t) 0,5 ) is the transference speed of the compound 
in the floor/wall/ceiling. 

The volume (V c ) of the chemical in the space portion 
is determined from the floor size (L„ x W 4 ), spray rate, 
content (C sd ) of the solvent (Xylene) contained in the 
product, dilution (X sol ) of the product, and evaporation 
constant (R d ) of the chemical in the space portion determined 
by expression (87). The volume (V c ) of the space portion 
decreases until the particle-constituting solvent in the 
chemical completely evaporates . 

Since the suspended particles are mostly made of 
moisture (particle-constituting solvent), their volume 
decreases as time passes. Since the suspended particles 
reduce their volume till the particle-constituting solvent 
completely evaporates, using diameter (dj of suspended 
particles determined by expression (92) or (92') and the 
above-mentioned evaporation constant ( a ) , the volume 
change (dV i /dt) of suspended particles is expressed as 
follows : 
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Assuming that the floor /wall/ceiling before spraying 
the chemical is like a thin film and that, as the chemical 
infiltrates into the film, its thickness increases so as 
5 to enhance the volume thereof; from the diffusion 

coefficient (D p ) of the compound in the floor/wall/ceiling, 
the volume (Vj) of floor /wall/ceiling is represented as 
follows : 

V } = 2(D p t)°- 5 Aj ...(100) 

10 Fugacity capacity (Z c ) of the chemical in the space 

portion is initially expressed by K ow C s /P s as with the 
above-mentioned Fugacity of the floor /wall/ceiling. When 
the chemical begins evaporating, Fugacity capacity (Z c ) 
gradually approaches C s /P 3 (Fugacity capacity of water). 

15 Accordingly, from the vapor pressure (P 3 ) , water solubility 

(C s ), and octanol/water distribution coefficient (K ow ) of the 
compound, and property change constant (a) of the product, 
Fugacity capacity (Z c ) is expressed as follows: 
Z c = (e- al K ow + 1 - e°')C s I P s ...(101) 

20 Here, the property change constant (a) of the product 

can be computed according to the evaporation constant (R d ) 
of the chemical in the space portion determined by expression 
(87), content (C 3d ) of the solvent (Xylene) contained in the 
product, and dilution (X sol ) of the product- Namely, 

25 Fugacity capacity (Z c ) of the chemical in the space portion 



223 



SCC 97-06 



is also expressed as: 

Z c = [(ratio of organic solvent in diluted liquid) K ow 
+ (ratio of water in diluted liquid)] * C s /P s 
By simultaneously forming this expression and the 
5 above-mentioned expression (101), the property change 

constant (a) of the product is computed, 
(ii) Fugacity Calculation (S3) 

The behavior of the compound in the carpet space 
portion is expressed in the form of differential equation 
10 concerning Fugacity (f c ) as: 

^T V c Z c = Ys^n-MViZifJHu Deposition^) ...(102) 
F - change 

+ R d L 4 W 4 Z c f V- change 

5 

- ^D ck (f c - f k ) Transference^) 
-KV c Z c f c Degradation 

Here, the terms of Deposition (i), V-change, Transference 
(4,5), and Degradation respectively indicate the attachment 
accompanying the falling of the suspended particles to the 

15 space portion, volume change (decrease over time) of the 

chemical in the space portion, amount of transference of 
the compound between the space portion and the air and floor 
(carpet fiber), and amount of degradation of the compound 
in the space portion. Here, Fugacity (f c ) is effective till 

2 0 the chemical in the space portion completely evaporates. 

Also, 7 5 indicates the ratio of suspended particles 
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attaching to the floor. 

The behavior of the compound in the large and small 
particles is expressed in the form of differential equation 
concerning Fugacity (f £ , i = 2, 3) as: 

^■ccd i Z i f i V- change ...(103) 

-Z). 4 (/ £ - / 4 ) Transference(A) 

-Ky.^ i f i Degradation 

Here, the terms of V-change, Transference (4), and 
Degradation respectively indicate the volume change 
(decrease over time) of the suspended particles, amount of 
transference of the compound between the particles and the 

10 air, and amount of degradation of the compound in the 

particles. Here, with the lapse of time, particles are 
absorbed by the carpet space portion and floor ( carpet fiber ) , 
whereby Fugacity (f^ is lost. 

The behavior of the compound in the indoor air is 

15 expressed in the form of differential equation concerning 

Fugacity (f 4 ) as: 

^-V.Z 4 = - GV.Z.f, 
dt 4 4 

F - change ~D c4 (f 4 - f c ) 
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dt ^' 



Ventilation ...(104) 
Transfer en ce(c) 
Transfer en ce{i) 

Transfer en ce(j) 
Degradation 
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Here, the terms of Ventilation, Transference (c), 
Transference (i), Transference (j), and Degradation 
respectively indicate the amount of discharge of the 
compound outdoors, amount of transference of the compound 
5 between the air and the space portion, amount of transference 

of the compound between the air and the suspended particles, 
amount of transference of the compound between the air and 
the floor /wall/ceiling, and amount of degradation of the 
compound in the air. 
10 The behavior of the compound in the floor (carpet 

fiber) is expressed in the form of differential equation 
concerning Fugacity (f 5 ) as: 

^y 5 Z 5 = 7s j?«,-vAsKZ ; Deposition^) ...(105) 
F - change 

- ^D p ItA a5 Z 5 f 5 V - change 
-D c5 (f 5 -f c ) Transference(c) 

- D 45 (/ 5 - f 4 ) Transference^) 

-K 5 V s Z 5 f s Degradation 
Here, the terms of Deposition (i), V-change, Transference 
15 Transference (4), and Degradation respectively 

indicate the attachment of the suspended particles 
accompanying falling to the carpet fiber, volume change 
(increase over time) of the floor, amount of transference 
of the compound between the carpet fiber and the space 
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portion, amount of transference of the compound between the 
carpet fiber and the air, and amount of degradation of the 
compound in the carpet fiber. 

The behavior of the compound in the wall and ceiling 
is expressed in the form of differential equation concerning 
Fugacity (f j7 j = 6, 7) as: 



dt 

F - change 



dt v j z j = ~4 D P ltA i Z jfi V -change ...(106) 



+ r . J? n i v i V i Z i f : IH {i) Deposition 

-D 4j ( fj - f 4 ) Transference 

- KjVjZj fj Degradation 
Here, the terms of V-change, Transference (4), and 
Degradation respectively indicate the volume change 
(increase over time) of the wall/ceiling, amount of 
transference of the compound between the wall/ceiling and 
the air, and amount of degradation of the compound in the 
wall/ceiling. Also, y ^ is the ratio of suspended particles 
attaching to the wall and ceiling. 

The above-mentioned seven kinds of differential 
equations (102) to (106) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
(f 2 to f 7 , and f c ) . 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
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width (dt) which is a minute time unit. Namely, the 
estimation nick time width is used such that solutions of 
the simultaneous differential equations are initially 
determined at a time (t 0 ), and then solutions of the 
5 simultaneous differential equations are determined at a 

time (t 0 + dt) to which the estimation nick time width is 
added. As solutions are obtained while estimation nick time 
widths are successively added, temporally changing 
Fugacity can be determined. Theoretically, as the set time 

10 of the estimation nick time width is shorter, more accurate 

solutions can be obtained, though necessitating a very long 
calculation time. By contrast, when the set time is too long, 
solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 

15 nick time width is set shorter when a very large change occurs 

in a chemical, whereas it is set longer when there is no 
large change. 

Specifically, mass balance is always confirmed such 
that the amount of input of the chemical and the resulting 

2 0 solution coincide with each other, and the estimation nick 

time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 
balance starts fluctuating. For example, when the 
fluctuation of mass balance is set to an accuracy of ±5%, 

2 5 the estimation nick time width is always set so as to 

constantly satisfy the relationship of: 
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compound input amount/ ( existing amount + degrading 
amount + discharging amount) = 0.95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 
5 Here, the above-mentioned compound input amount is 

determined by the floor size (L 4 x W 4 ) , compound content (C a ) , 
and spray rate. Since the temporal amounts of compound in 
the seven kinds of media are determined by the simultaneous 
differential equations , they are summed up so as to compute 
10 the existing amount as shown in the following. Also, the 

degrading amount and discharging amount are as follows : 

3 7 

existing amount = V c ZJ c fi V i Z i + fk V k Z k 

3 7 

degrading amount = K c V c Z c f c + g K.n, f.V.Z., + ^ K k fk V k Z k ' ' ' O- 0 ?) 

discharging amount = GV 4 Z 4 / 4 + ^Gn ; KZ ; / ; 

With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 

15 explained. 

First, an initial value of estimation nick time width 
(dt ) is inputted (S21 ) . Then, an upper limit set value (e.g. , 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 

2 0 is the lower limit of difference in mass balance are inputted 

(S22) . Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
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t = t + dt (or t„ + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27) . If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
5 whether the difference between the mass balance at t = t 

(or t 0 ) and the mass balance at t = t + dt (or t 0 + dt) is 
at least the upper limit set value or not (S28). If it is 
judged to be at least the upper limit set value at S28, then 
solutions become more accurate when the estimation nick time 

10 width (dt) is made shorter. In this case, the estimation 

nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S2 8, it is judged 
at S3 0 whether the difference between the mass balance at 

15 t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 

dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S30, since solutions are not 
influenced by longer estimation nick time width (dt), the 

2 0 estimation nick time width (dt) is doubled so as to change 

its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width (dt) or not (S32). When the estimation nick time width 

25 (dt) is not greater than the maximum value at S3 2, since 

solutions do not diverge, the estimation nick time width 
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(dt) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S3 2, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S30, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
are continued without changing the estimation nick time 
width (dt). 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S2 7 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
balance fluctuation, Fugacity can be computed accurately 
and efficiently. 

(iii) Computation of Temporal Concentration in Indoor Air 
and Floor Residual Amount (S4) 
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The temporal concentration of the compound in the 
indoor air is computed when Fugacity (f 3 ) of the small 
particles, determined by the above-mentioned item (ii), 
multiplied by Fugacity capacity (Z 3 ), and Fugacity (f 4 ) of 
5 the indoor air multiplied by Fugacity capacity (Z 4 ) are 

summed up. Here, in the case where the large particles may 
be inhaled by a human body depending on the kind of chemical, 
calculation is effected while taking account of Fugacity 
(f 2 ) of the large particles as well. 

10 The floor residual amount of the compound is computed 

when Fugacity (f 5 ) of the floor is multiplied by Fugacity 
capacity (Z 5 ). In the case where the solvent is water, in 
particular, it is preferably computed by summing up the 
above-mentioned computed value (f s , Z 5 ) and the value 

15 obtained when Fugacity (f c ) of the space portion is 

multiplied by Fugacity capacity (Z c ). 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 

2 0 indoor air forms a curve shown in Fig. 3A, for example. This 

concentration curve is integrated, an accumulated 
concentration of the compound during a specific period {t 1 
to t 2 ) is determined (see Fig. 3B), and the mean concentration 
in the indoor air is computed from thus determined value. 

25 While an arbitrary period is set as the specific period 

depending on the object, an appropriate period is usually 
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set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 
concentration in the indoor air, amount of respiration, and 
exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] x amount of 
respiration [m 3 /kg/min] * exposure time [min/day] 
is effected. Here, as the above-mentioned amount of 
respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 

The inhalation safety coefficient is computed from a 
non-influential amount concerning inhalation toxicity 
examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 

inhalation safety coefficient = 

inhalation non-influential amount [mg/kg/day]/ 
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estimated exposure amount in inhalation [mg/kg/day] 

(v) Calculation of Estimated Percutaneous Exposure Amount 

and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned floor residual amount forms a 
5 curve such as that shown in Fig. 4A, for example. This 

residual amount curve is integrated, the accumulated 
residual amount of the compound during a specific period 
(t 1 to t 2 ) is determined (see Fig. 4B) , and the mean residual 
amount is computed from thus determined value. While an 

10 arbitrary period is set as the specific period depending 

on the object, an appropriate period is usually set in view 
of the method of use of the product and the test period of 
toxicity data. 

Then, according to the mean floor residual amount, 

15 skin attachment ratio, contact area, and body weight, the 

estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean floor residual amount [rag/m 2 ] * skin attachment ratio 

20 [%] x contact area [m 2 /day ] ) /body weight [kg] 

is effected. Here, as the contact area, a published value 
(e.g., 4 [m 2 /day]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor where the compound exists . 

2 5 As this value, a published value or a value experimentally 

obtained from a model may be used. 
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A model experiment method for the skin attachment 
ratio is as follows. A weight ( 8 cm x 8 cm x 8 cm; 4.2 kg) 
is placed on a denim cloth (8 cm x 10 cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants. 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 

percutaneous non- influential amount [mg/kg/day]/ 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
are not many published values . Accordingly, a more accurate 
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value can be determined from the estimated percutaneous 
exposure amount, and oral non- influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 
5 percutaneous safety coefficient = oral non- 

influential amount [mg/kg/day] / (estimated percutaneous 
exposure amount [mg/kg/day] x percutaneous absorption ratio 
[%]) 

Here, when the percutaneous absorption ratio is unknown, 
10 employed is a national guideline (e.g., 10%) which usually 

exists . 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 

15 obtained in the above-mentioned item (v), hand surface area 

ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 

20 estimated percutaneous exposure amount [mg/kg/day] * hand 

surface area ratio [%] * oral transference ratio [%] 
is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 

25 transference ratio is a hypothetical value, which is set 

to 100%, for example. 
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In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
amount from tableware or food to mouth be added to the 
5 estimated oral exposure amount from hand to mouth to yield 

the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 
the amount of the harmful material remaining in tableware, 
=- 10 tableware use area which is the sum of tableware surface 

areas, and oral transference ratio from tableware, 
1 calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] * tableware use area 
15 [m 2 /day] x oral transference ratio [ % ] /body weight [kg] 

is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount * tableware 
contamination ratio ) . As the tableware contamination ratio , 
an actually measured value (e.g., 9%) or a hypothetical value 
20 may be used. 

The oral safety coefficient is computed from the 
non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 
25 as: 

oral safety coefficient = 
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oral non-influential amount [mg/kg/day]/ estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
this embodiment, when a chemical containing the above- 
mentioned compound is sprayed over the whole floor area, 
differential equations concerning Fugacity of the space 
portion between carpet fibers, two kinds of large and small 
particles, indoor air, whole carpet fiber, wall, and ceiling 
are simultaneously formed and are solved, and the indoor 
behavior of the compound is estimated according to thus 
obtained solution. Here, the estimation nick time width is 
variably set, while constantly confirming mass balance of 
the compound indoors after the spraying, so that the amount 
of input of the chemical indoors and the resulting solution 
coincide with each other. 

Accordingly, since mass balance of the compound after 
the spraying is always monitored such that the amount of 
input of the compound indoors and the resulting solution 
coincide with each other, thereby varying the estimation 
nick time width; the estimation nick time width is set longer 
when the mass balance fluctuates a little, whereas it is 
set shorter when the mass balance starts fluctuating greatly . 
Namely, when solving simultaneous differential equations 
including a parameter accompanying temporal change, the 
estimation nick time width is automatically set in response 
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to the fluctuation in mass balance. Consequently, when 
processed by a computer, an accurate solution can be obtained 
in a short time. 

The method of evaluating safety of a pesticidal 
5 compound in accordance with this embodiment uses the 

estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is sprayed over the whole floor surface. 

Accordingly, the safety of the compound with respect 

10 to the human body can be evaluated accurately in a short 

time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 

15 conforming to the aimed object. 

Though Fugacity is determined by use of Runge- 
Kutta-Gill method in this embodiment, other methods may be 
used for solving differential equations . Runge-Kutta-Gill 
method, however, is preferably used since a program for the 

20 above-mentioned differential equations can be easily made 

by Basic. Also in the case where differential equations are 
solved by a method other than Runge-Kutta-Gill method, 
similar effects can be obtained when the estimation nick 
time width is set as mentioned above. 

25 Finally, a computer program product (storage medium) 

employed in each of the above-mentioned embodiments will 
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be explained. Fig. 19 is a conf igurational view showing a 
storage medium 20 employed in each of the above-mentioned 
embodiments of the present invention. As shown in Fig. 19, 
the storage medium 20 comprises a program area 20a for 
5 storing a program. This program area 20a stores a 

pesticidal compound indoor behavior simulation program 22. 
The pesticidal compound indoor behavior simulation program 
22 comprises a main program 22a for supervising the 
processing; a differential equation forming program 22b for 

10 (SI) inputting a primary condition necessary for forming 

a differential equation concerning a fugacity of the 
compound in each medium in response to an indoor environment 
divided into predetermined media and (S2) calculating a 
secondary condition so as to form the above-mentioned 

15 differential equation; a Fugacity calculating program 22c 

for (S3) determining the fugacity of the compound from the 
above-mentioned differential equation; an indoor behavior 
calculating program 22d for (S4) computing temporal 
concentration in indoor air and residual amount, which 

2 0 constitute an indoor behavior of the compound, from thus 

calculated fugacity; a safety evaluating program 22e for 
(S5 to S17) evaluating the safety of the compound with 
respect to a human body according to the temporal 
concentration in the indoor air and residual amount; and 

25 a minute time unit changing program 22f for (S21 to S35) 

changing, in response to a fluctuation in mass balance of 



240 



SCC 97-06 



the compound, the minute time unit used when solving the 
differential equation. The above-mentioned letters and 
numerals within parentheses refer to the corresponding 
steps in Fig. 1 or 5 . 
5 As the storage medium 20, a disk-shaped storage medium 

such as flexible disk or CD-ROM is used, for example. 
Alternatively, a tape-like storage medium such as magnetic 
tape may be used as well. 

The storage medium 20 may be used in an information 

10 processing apparatus shown in Figs. 20 and 21. Namely, the 

information processing apparatus 30 comprises a medium 
drive unit 32, which can accommodate the storage medium 20. 
As the storage medium 2 0 is thus accommodated, the 
information stored therein can be accessed by the medium 

15 drive unit 32. Consequently, the pesticidal compound 

indoor behavior simulation program 22 stored in the program 
area 20a can be executed by the information processing 
apparatus 30. 

This information processing apparatus is configured 

20 as follows . First, it comprises the above-mentioned medium 

drive unit 32, an image memory 34 for storing image data 
indicating simulation results and the like; a work memory 
(internal memory) 36 in which an operation system (OS) is 
resident; and a display 38 which is a display means. Also, 

2 5 it comprises an input device 44 which is an input means having 

a mouse 40 and a keyboard 42 for receiving input data such 
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as the primary condition; a printer 4 6 for outputting image 
data and the like; and a CPU 48 for controlling execution 
of the pesticidal compound indoor behavior simulation 
program 22 and the like. 

As the medium drive unit, in conformity to the storage 
medium 20, a flexible disk unit, CD-ROM drive unit, magnetic 
tape drive unit, or the like may be used. 

In the following, the processing of the pesticidal 
compound indoor behavior simulation program 22 stored in 
the program area 2 0a of the storage medium 2 0 will be 
explained. This processing is performed by executing the 
pesticidal compound indoor behavior simulation program 22 
read out by the medium drive unit 32. Upon this execution, 
the main program 22a of the pesticidal compound indoor 
behavior simulation program 22 is activated. 

The subsequent processing operations of the 
differential equation forming program 22b, Fugacity 
calculating program 22c, indoor behavior calculating 
program 22d, safety evaluating program 22e, and minute time 
unit changing program 22f correspond to those explained in 
the above-mentioned embodiments and Figs . 1 and 5 . 

As explained in the foregoing, the computer program 
product (storage medium) employed in each of the above- 
mentioned embodiments of the present invention has a program 
for initially dividing the environment into predetermined 
media (constituents), simultaneously forming differential 
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equations concerning Fugacity of the compound in these 
medium, and then solving them, thereby estimating an indoor 
behavior of a compound when a chemical containing the 
compound is residually sprayed. 

Consequently, the environment is assumed to be a 
plurality of media, and exchanges of the compound between 
media and the like are taken into account, whereby simulation 
results close to the actual behavior of the compound can 
be obtained. 

The estimation nick time width set when solving the 
simultaneous differential equations is variably set, while 
constantly confirming mass balance of the compound indoors 
after the residual spraying, so that the amount of input 
of the chemical indoors and the resulting solution coincide 
with each other. 

Accordingly, since mass balance of the compound after 
the residual spraying is always monitored such that the 
amount of input of the compound indoors and the resulting 
solution coincide with each other, thereby varying the 
estimation nick time width; the estimation nick time width 
is set longer when the mass balance fluctuates a little, 
whereas it is set shorter when the mass balance starts 
fluctuating greatly. Namely, when solving simultaneous 
differential equations including a parameter accompanying 
temporal change, the estimation nick time width is 
automatically set in response to the fluctuation in mass 
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balance. Consequently, when processed by a computer, an 
accurate solution can be obtained in a short time. 

Also, the computer program product (storage medium) 
employed in each of the above-mentioned embodiments of the 
present invention has a program for evaluating, by use of 
the above-mentioned estimated results, safety of a 
pesticidal compound with respect to the human body when the 
chemical is residually sprayed. 

Accordingly, the safety of the compound with respect 
to the human body can be evaluated accurately in a short 
time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 
Industrial Applicability 

In the method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention, the indoor environment is divided into 
predetermined media, and exchanges of the compound between 
media and the like are taken into account, whereby simulation 
results similar to the actual behavior of the compound can 
be obtained. Also, when solving simultaneous differential 
equations including a parameter accompanying a temporal 
change, a minute time unit can be automatically set in 
response to the fluctuation in mass balance, whereby an 
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accurate solution can be obtained in a short time when 
processed by a computer. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
5 invention further comprises a step of evaluating, according 

to the indoor behavior of the compound, safety of the 
compound with respect to a human body, whereby the safety 
of the pesticidal compound can be evaluated accurately in 
a short time. Accordingly, when formulating a chemical such 
10 as insecticide containing the compound, simulation can be 

easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 
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CLAIMS 

1. A method of simulating an indoor behavior of a 
pesticidal compound, said simulation method comprising: 

a step of dividing an indoor environment into 
5 predetermined media and forming a differential equation 

concerning a fugacity of said compound in each of said media; 

a step of determining the fugacity of said compound 
in each of said media from said differential equation; 

a step of determining the indoor behavior of said 
10 compound from the fugacity of said compound in each of said 

media; and 

a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving said differential equation. 
15 2. A simulation method according to claim 1, 

further comprising a step of evaluating safety of said 
compound with respect to a human body according to the indoor 
behavior of said compound. 

3. A simulation method according to claim 1 or 2 , 
2 0 wherein said compound is introduced into an indoor space 

as a solution containing said compound is residually 
sprayed; and 

wherein said media are a spraying site, suspended 
particles which are divided into at least one kind according 
25 to size, indoor air, a floor, a wall, and a ceiling. 

4. A simulation method according to claim 3, 
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wherein said differential equation at said spraying site 
is a differential equation stating a relationship among 
temporal change of fugacity of said compound at said spraying 
site, temporal change in volume of said spraying site, amount 
of attachment of said suspended particles to said spraying 
site, amount of transference of said compound between said 
spraying site and another medium, and change in amount of 
degradation of said compound at said spraying site; 

wherein said differential equation in said suspended 
particles is a differential equation stating a relationship 
among temporal change of fugacity of said compound in said 
suspended particles, temporal change in volume of said 
suspended particles, amount of transference of said 
compound between said suspended particles and another 
medium, and change in amount of degradation of said compound 
in said suspended particles; 

wherein said differential equation in said indoor air 
is a differential equation stating a relationship among 
temporal change of fugacity of said compound in said indoor 
air, amount of discharge of said compound outdoors, amount 
of transference of said compound between said indoor air 
and another medium, and change in amount of degradation of 
said compound in said indoor air; 

wherein said differential equation at said floor is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said floor, 
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temporal change in volume of said floor, amount of attachment 
of said suspended particles to said floor, amount of 
transference of said compound between said floor and another 
medium, and change in amount of degradation of said compound 
at said floor; 

wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of attachment of said 
suspended particles to said wall, amount of transference 
of the compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
attachment of said suspended particles to said ceiling, 
amount of transference of said compound between said ceiling 
and another medium, and change in amount of degradation of 
said compound at said ceiling - 

5. A simulation method according to claim 1 or 2 , 
wherein said compound is introduced into an indoor space 
as a solution containing said compound is spatially sprayed; 
and 

wherein said media are suspended particles which are 
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divided into at least one kind according to size, indoor 
air, a floor, a wall, and a ceiling. 

6. A simulation method according to claim 5, 
wherein said differential equation in said suspended 
5 particles is a differential equation stating a relationship 

among temporal change of fugacity of said compound in said 
suspended particles, temporal change in volume of said 
suspended particles, amount of transference of said 
compound between said suspended particles and another 
10 medium, and change in amount of degradation of said compound 

in said suspended particles; 

wherein said differential equation in said indoor air 
is a differential equation stating a relationship among 
temporal change of fugacity of said compound in said indoor 
15 air, amount of discharge of said compound outdoors, amount 

of transference of said compound between said indoor air 
and another medium, and change in amount of degradation of 
said compound in said indoor air; 

wherein said differential equation at said floor is 
2 0 a differential equation stating a relationship among 

temporal change of fugacity of said compound at said floor, 
temporal change in volume of said floor, amount of attachment 
of said suspended particles to said floor, amount of 
transference of said compound between said floor and another 
2 5 medium, and change in amount of degradation of said compound 

at said floor; 
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wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of attachment of said 
suspended particles to said wall, amount of transference 
of said compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
attachment of said suspended particles to said ceiling, 
amount of transference of said compound between said ceiling 
and another medium, and change in amount of degradation of 
said compound at said ceiling. 

7. A simulation method according to claim 1 or 2, 
wherein said compound is introduced into an indoor space 
as a solution containing said compound is heated to vaporize; 
and 

wherein said media are condensed particles which are 
divided into at least one kind according to generation and 
extinction, high-concentration air, medium-concentration 
air, low-concentration air, a floor, a wall, and a ceiling 
which is divided into at least one kind according to compound 
concentration . 
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8. A simulation method according to claim 7 , 
wherein said differential equation in said condensed 
particles is a differential equation stating a relationship 
among temporal change of fugacity of said compound in said 
condensed particles , temporal change in volume of said 
condensed particles, amount of transference of said 
compound between said condensed particles and another 
medium, and change in amount of degradation of said compound 
in said condensed particles; 

wherein said differential equation in said high- 
concentration air is a differential equation stating a 
relationship among temporal change of fugacity of said 
compound in said high-concentration air, amount of 
discharge of said compound, amount of transference of said 
compound between said high-concentration air and another 
medium, and change in amount of degradation of said compound 
in said high-concentration air; 

wherein said differential equation in said medium- 
concentration air is a differential equation stating a 
relationship among temporal change of fugacity of said 
compound in said medium-concentration air, amount of 
transference of said compound between said medium- 
concentration air and another medium, and change in amount 
of degradation of said compound in said medium- 
concentration air; 

wherein said differential equation in said low- 
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concentration air is a differential equation stating a 
relationship among temporal change of fugacity of said 
compound in said low-concentration air, amount of discharge 
of said compound outdoors, amount of transference of said 
compound between said low-concentration air and another 
medium, and change in amount of degradation of said compound 
in said low-concentration air; 

wherein said differential equation at said floor is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said floor, 
temporal change in volume of said floor, amount of 
transference of said compound between said floor and another 
medium, and change in amount of degradation of said compound 
at said floor; 

wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of transference of 
said compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
transference of said compound between said ceiling and 
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another medium, and change in amount of degradation of said 
compound at said ceiling. 

9. A simulation method according to claim 1 or 2, 
wherein said compound is introduced into an indoor space 
as a solution containing said compound is sprayed over the 
whole floor; and 

wherein said media are suspended particles which are 
divided into at least one kind according to size, indoor 
air, a floor, a wall, and a ceiling. 

10. A simulation method according to claim 9, 
wherein said differential equation in said suspended 
particles is a differential equation stating a relationship 
among temporal change of fugacity of said compound in said 
suspended particles, temporal change in volume of said 
suspended particles, amount of transference of said 
compound between said suspended particles and another 
medium, and change in amount of degradation of said compound 
in said suspended particles; 

wherein said differential equation in said indoor air 
is a differential equation stating a relationship among 
temporal change of fugacity of said compound in said indoor 
air, amount of discharge of said compound outdoors, amount 
of transference of said compound between said indoor air 
and another medium, and change in amount of degradation of 
said compound in said indoor air; 

wherein said differential equation at said floor is 
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a differential equation stating a relationship among 
temporal change of fugacity of said compound at said floor, 
temporal change in volume of said floor, amount of attachment 
of said suspended particles to said floor, amount of 
transference of said compound between said floor and another 
medium, and change in amount of degradation of said compound 
at said floor; 

wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of attachment of said 
suspended particles to said wall, amount of transference 
of said compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
attachment of said suspended particles to said ceiling, 
amount of transference of said compound between said ceiling 
and another medium, and change in amount of degradation of 
said compound at said ceiling - 

11- A simulation method according to any of claims 
3 to 10, wherein said floor is constituted by a rug having 
ears of fiber; and 
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wherein a space between said ears is added to said 

media . 

12. A simulation method according to claim 11, 
wherein said differential equation in the space between said 

5 ears is a differential equation stating a relationship among 

temporal change of fugacity of said compound in the space 
between said ears, temporal change in volume of the solution 
containing said compound in the space between said ears, 
amount of attachment of said compound into the space portion 
10 between said ears by falling, amount of transference of said 

compound between the space portion between said ears and 
another medium, and change in amount of degradation of said 
compound in the space portion between said ears . 

13. A computer program product to be used together 
15 with an information processing apparatus comprising input 

means for receiving a data input from outside, display means , 
and readout means for reading out information from a 
computer-usable storage medium; 

said computer program product comprising a 

20 computer-usable storage medium which has a program area for 

storing a program and has a computer-readable program 
materialized in said storage medium for causing, according 
to data input from said input means, said display means to 
display a result of simulation of an indoor behavior of a 

25 pesticidal compound; 

said computer program product comprising: 
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in said program area, 

a program for dividing an indoor environment into 
predetermined media and forming a differential equation 
concerning a fugacity of said compound, 
5 a program for determining the fugacity of said 

compound in each of said media from said differential 
equation, 

a program for determining the indoor behavior of said 
compound from the fugacity of said compound in each of said 
10 media, and 

a program for changing, in response to a fluctuation 
in mass balance of said compound indoors, a minute time unit 
used when solving said differential equation. 

14. A computer program product according to claim 
15 13, further comprising, in said program area, 

a program for evaluating safety of said compound with 
respect to a human body according to the indoor behavior 
of said compound. 

15. A computer program product according to claim 
20 13 or 14, wherein said compound is introduced into an indoor 

space as a solution containing said compound is residually 
sprayed; and 

wherein said media are a spraying site, suspended 
particles which are divided into at least one kind according 
25 to size, indoor air, a floor, a wall, and a ceiling. 

16. A computer program product according to claim 
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15, wherein said differential equation at said spraying site 
is a differential equation stating a relationship among 
temporal change of f ugacity of said compound at said spraying 
site, temporal change in volume of said spraying site, amount 
5 of attachment of said suspended particles to said spraying 

site, amount of transference of said compound between said 
spraying site and another medium, and change in amount of 
degradation of said compound at said spraying site; 

wherein said differential equation in said suspended 

10 particles is a differential equation stating a relationship 

among temporal change of fugacity of said compound in said 
suspended particles, temporal change in volume of said 
suspended particles, amount of transference of said 
compound between said suspended particles and another 

15 medium, and change in amount of degradation of said compound 

in said suspended particles; 

wherein said differential equation in said indoor air 
is a differential equation stating a relationship among 
temporal change of fugacity of said compound in said indoor 

20 air, amount of discharge of said compound outdoors, amount 

of transference of said compound between said indoor air 
and another medium, and change in amount of degradation of 
said compound in said indoor air; 

wherein said differential equation at said floor is 

25 a differential equation stating a relationship among 

temporal change of fugacity of said compound at said floor, 
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temporal change in volume of said floor, amount of attachment 
of said suspended particles to said floor, amount of 
transference of said compound between said floor and another 
medium, and change in amount of degradation of said compound 
at said floor; 

wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of attachment of said 
suspended particles to said wall, amount of transference 
of the compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
attachment of said suspended particles to said ceiling, 
amount of transference of said compound between said ceiling 
and another medium, and change in amount of degradation of 
said compound at said ceiling. 

17. A computer program product according to claim 
13 or 14, wherein said compound is introduced into an indoor 
space as a solution containing said compound is spatially 
sprayed; and 

wherein said media are suspended particles which are 
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divided into at least one kind according to size, indoor 
air, a floor, a wall, and a ceiling. 

18. A computer program product according to claim 
17 wherein said differential equation in said suspended 
5 particles is a differential equation stating a relationship 

among temporal change of fugacity of said compound in said 
suspended particles, temporal change in volume of said 
suspended particles, amount of transference of said 
compound between said suspended particles and another 

10 medium, and change in amount of degradation of said compound 

in said suspended particles; 

wherein said differential equation in said indoor air 
is a differential equation stating a relationship among 
temporal change of fugacity of said compound in said indoor 

15 air, amount of discharge of said compound outdoors, amount 

of transference of said compound between said indoor air 
and another medium, and change in amount of degradation of 
said compound in said indoor air; 

wherein said differential equation at said floor is 

20 a differential equation stating a relationship among 

temporal change of fugacity of said compound at said floor, 
temporal change in volume of said floor, amount of attachment 
of said suspended particles to said floor, amount of 
transference of said compound between said floor and another 

25 medium, and change in amount of degradation of said compound 

at said floor; 
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wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of attachment of said 
suspended particles to said wall, amount of transference 
of said compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
attachment of said suspended particles to said ceiling, 
amount of transference of said compound between said ceiling 
and another medium, and change in amount of degradation of 
said compound at said ceiling. 

19. A computer program product according to claim 
13 or 14, wherein said compound is introduced into an indoor 
space as a solution containing said compound is heated to 
vaporize; and 

wherein said media are condensed particles which are 
divided into at least one kind according to generation and 
extinction, high-concentration air, medium-concentration 
air, low-concentration air, a floor, a wall, and a ceiling 
which is divided into at least one kind according to compound 
concentration . 
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20. A computer program product according to claim 
19, wherein said differential equation in said condensed 
particles is a differential equation stating a relationship 
among temporal change of fugacity of said compound in said 
condensed particles, temporal change in volume of said 
condensed particles, amount of transference of said 
compound between said condensed particles and another 
medium, and change in amount of degradation of said compound 
in said condensed particles; 

wherein said differential equation in said high- 
concentration air is a differential equation stating a 
relationship among temporal change of fugacity of said 
compound in said high-concentration air, amount of 
discharge of said compound, amount of transference of said 
compound between said high-concentration air and another 
medium, and change in amount of degradation of said compound 
in said high-concentration air; 

wherein said differential equation in said medium- 
concentration air is a differential equation stating a 
relationship among temporal change of fugacity of said 
compound in said medium-concentration air, amount of 
transference of said compound between said medium- 
concentration air and another medium, and change in amount 
of degradation of said compound in said medium- 
concentration air; 

wherein said differential equation in said low- 
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concentration air is a differential equation stating a 
relationship among temporal change of fugacity of said 
compound in said low-concentration air, amount of discharge 
of said compound outdoors, amount of transference of said 
compound between said low-concentration air and another 
medium, and change in amount of degradation of said compound 
in said low-concentration air; 

wherein said differential equation at said floor is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said floor, 
temporal change in volume of said floor, amount of 
transference of said compound between said floor and another 
medium, and change in amount of degradation of said compound 
at said floor; 

wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of transference of 
said compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
transference of said compound between said ceiling and 
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another medium, and change in amount of degradation of said 
compound at said ceiling. 

21. A computer program product according to claim 
13 or 14, wherein said compound is introduced into an indoor 
5 space as a solution containing said compound is sprayed over 

the whole floor; and 

wherein said media are suspended particles which are 
divided into at least one kind according to size, indoor 
air, a floor, a wall, and a ceiling. 

10 22. A computer program product according to claim 

21, wherein said differential equation in said suspended 
particles is a differential equation stating a relationship 
among temporal change of fugacity of said compound in said 
suspended particles, temporal change in volume of said 

15 suspended particles, amount of transference of said 

compound between said suspended particles and another 
medium, and change in amount of degradation of said compound 
in said suspended particles; 

wherein said differential equation in said indoor air 

2 0 is a differential equation stating a relationship among 

temporal change of fugacity of said compound in said indoor 
air, amount of discharge of said compound outdoors, amount 
of transference of said compound between said indoor air 
and another medium, and change in amount of degradation of 

25 said compound in said indoor air; 

wherein said differential equation at said floor is 
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a differential equation stating a relationship among 
temporal change of fugacity of said compound at said floor , 
temporal change in volume of said floor, amount of attachment 
of said suspended particles to said floor, amount of 
transference of said compound between said floor and another 
medium, and change in amount of degradation of said compound 
at said floor; 

wherein said differential equation at said wall is a 
differential equation stating a relationship among temporal 
change of fugacity of said compound at said wall, temporal 
change in volume of said wall, amount of attachment of said 
suspended particles to said wall, amount of transference 
of said compound between said wall and another medium, and 
change in amount of degradation of said compound at said 
wall; and 

wherein said differential equation at said ceiling is 
a differential equation stating a relationship among 
temporal change of fugacity of said compound at said ceiling, 
temporal change in volume of said ceiling, amount of 
attachment of said suspended particles to said ceiling, 
amount of transference of said compound between said ceiling 
and another medium, and change in amount of degradation of 
said compound at said ceiling - 

23. A computer program product according to any of 
claims 15 to 22, wherein said floor is constituted by a rug 
having ears of fiber; and 
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wherein a space between said ears is added to said 

media . 

24. A computer program product according to claim 
23, wherein said differential equation in the space between 
said ears is a differential equation stating a relationship 
among temporal change of fugacity of said compound in the 
space between said ears, temporal change in volume of the 
solution containing said compound in the space between said 
ears, amount of attachment of said compound into the space 
portion between said ears by falling, amount of transference 
of said compound between the space portion between said ears 
and another medium, and change in amount of degradation of 
said compound in the space portion between said ears. 
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ABSTRACT 

The present invention relates to a method of 
simulating an indoor behavior of a pesticidal compound, and 
aims to provide a simulation method which can process 
simultaneous differential equations necessary for 
simulation accurately in a short time, and to evaluate safety 
of the compound with respect to a human body. 

In order to achieve the above-mentioned object, the 
simulation method of the present invention comprises a step 
of dividing an indoor environment into predetermined media 
and forming a differential equation concerning a fugacity 
of the compound in each of the media; a step of determining 
the fugacity of the compound in each of the media from the 
differential equation; a step of determining the indoor 
behavior of the compound from the fugacity of the compound 
in each of the media; a step of changing, in response to 
a fluctuation in mass balance of the compound indoors, a 
minute time unit used when solving the differential 
equation; and a step of evaluating, according to the indoor 
behavior of the compound, safety of the compound with respect 
to the human body. 
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